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Abstract. An exact solution is obtained for the diffraction of a dipole field by a 
unidirectionally conducting semi-infinite plane screen. Double Laplace transforms are 
applied to Maxwell’s equations, and the defining conditions of the unidirectionality 
lead to an equation between two complex functions of two complex variables. This 
equation is solved by an extension of the usual function-theoretical method, and we 
can then express the electro-magnetic field in terms of certain complex triple integrals. 
These are transformed into real integrals, so that it is possible to discuss the field behavior 
in the neighborhood of the diffracting edge. The variation of singularity along the edge 
of the screen is given. 

1. Introduction. Diffraction by a unidirectionally conducting body is the subject of 
two recent investigations: Toraldo di Francia [1] has given approximate results, based 
on a physical discussion, for a unidirectional screen of small diameter, while Karp [2] 
has obtained an exact solution for the diffraction of a plane wave by a unidirectional 
semi-infinite screen. Provided only the far field were of interest, well-known reciprocity 
considerations would suffice to extend the plane wave result to the case of dipole inci- 
dence; but a discussion of the near field, and in particular of the physically interesting 
variation of the fields and currents along the edge of the unidirectional screen, requires 
the complete solution of the diffraction problem for a dipole, and this is the problem 
here considered. 

Important applications of the theory of unidirectional screens are: (1) to the measure- 
ment of the angular momentum of electromagnetic radiation, as in Toraldo di Francia’s 
work (see [1, 9, 10]); (2) to microwave problems involving unidirectionally conducting 
components; (3) to problems of propagation over anisotropic media. In any of these 
cases, the problem of this paper plays the role of a canonical problem, in the sense that 
our results as to edge behavior for a semi-infinite screen permit (by way of Keller’s 
geometrical theory of diffraction: see [11], and references given there to earlier work 
by Keller and his collaborators) the deduction of asymptotic results for a large class of 
unidirectional screens of finite size. 

Our analysis is based on a formulation of the diffraction of a dipole field by a uni- 
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directional semi-infinite screen as a Wiener-Hopf problem. Jones [3] observed that 
diffraction problems leading to Wiener-Hopf equations are advantageously treated by 
taking the transform of the differential equation before applying boundary conditions: 
extending Jones’ method to the case of dipole incidence amounts to little more than 
replacing single by double transforms. The known solution [4] for a perfectly conducting 
screen may very easily be derived in this manner. 

For the case of dipole incidence on a unidirectional semi-infinite screen, it is found 
after taking two successive Laplace transforms of the Maxwell equation 


(A,,. + ke = 0, 


that we may express the double transform of the electric field as an unknown vector 
function of the transform variables. Obtaining relations among the components of 
this vector function from the remaining Maxwell equations, and then applying the 
boundary condition and jump conditions derived from the unidirectionality, we show 
in Sec. 2 that the diffraction problem is equivalent to the solution of a single transform 
equation for two unknown complex functions. The transform equation also involves two 
independent complex variables and its solution therefore requires some modification of 
the usual function-theoretical considerations. This solution is derived in Sec. 3, with 
the result that the field components are expressed as Laplace inverses. 

In Sec. 4, it is found that the originals of these inverse transforms are integrals of a 
type introduced by Macdonald [5], plus additional terms which it is possible to transform 
into certain real integrals. These results permit us to give the variation of the near 
field as the diffracting edge is traversed. 

The results are summarized in Sec. 5, in the form of a theorem, and it is verified 
that all conditions of the problem are met. 

2. Derivation of the transform equation. We consider diffraction by a unidirection- 
ally conducting semi-infinite screen: z > 0, — © < y < o,z2 =O, where z, y, z form a 
right-hand rectangular coordinate system. Let e, , hy be the electric and magnetic field 
vectors of an incident dipole field, which we take as an electric dipole with axis perpen- 
dicular to the screen, while remarking that the method which follows is also applicable 
to an arbitrarily oriented dipole. Locating our dipole at (xo , Yo , 2.) with z > 0, we may 


describe the incident field by the Hertz vector (0, 0, II,), where II, = e “*/kR and 


R = [(x — a)? + (y — yo)? + (2 — 2%)’ ]'”. 
The corresponding electromagnetic field components are 


_ ( am, OM, OM, , Kn. ) 
\°  \éz dx ’ dz By’ Oz” ware 
] (1) 
h, = iad 4-0 — a, ‘ 0). 
OY OL 


Time dependence e‘“*‘ is understood. 

Now denote by e, + e, h, + h the total electric and magnetic fields resulting from 
the incidence of e€, , hy upon the given screen. Then the scattered field vectors e, h 
satisfy the time-harmonic Maxwell equations 

(A... + ke = 0, (2) 


V-e=0, (3) 
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—tph = V Xe, (4) 


subject to a set of conditions (boundary condition, two jump conditions and an edge 
condition) of which the first three are intended as a phenomenological description of 
the unidirectional conductivity of the diffracting screen. These conditions are con- 
veniently stated in terms of field components in the direction £ of conductivity and 


the direction » normal to &: 


és (5) 
fhe] = he(z + 0) — he(z — 0) = O, across the screen, (6) 
[e,] = 0, across the screen, (7) 
[h,] = 0, at the edge of the screen. (8) 


We also assume that e, h are integrable at the edge of the screen: this, with (8), yields 


a unique solution. 
Conditions at infinity complete our specification of the scattered field. We impose 


the familiar condition that e, h be exponentially damped solutions of the three-dimensional 

wave equation. The implied behavior at infinity is that of the elementary solution 

e '*“/kR, where we suppose k to have a negative imaginary part: k = k, — tk, (k, > 0). 
Now define the double Laplace transforms E, H of the scattered field vectors e, h: 


a a 


It follows from (2) that E satisfies 
a - 
(2 + k*)E = 0, (9) 
dz 
where K* = k’ + s? + #*. Our condition on e at infinity yields the range of validity 
of (9): we fix K by the choice 
(ke + 8° + €)""|,2100 = +h, 


and note that (9) is meaningful for | Ret | < k, , | Res| < | Im (k? + @)* |. In this 
domain of s, ¢ we write the solution of (9) (notice Im K < 0) as 


A(s, t) exp [—iK(z + 2)], z>0 (10) 
B(s, t) exp [¢K(z — 2)], z2<0. 


E 


We proceed to determine the transform of the scattered electric field by solving for 


A(s, t), B(s, t). 
It is convenient to rotate coordinates from 2, y to &, n. Let a) (0 < a, < 2/2) be the 
angle, measured clockwise, from the z-direction to the direction. Then 


& = X COSA) — YSINa , (11) 


XSIN a + Y COSA .« 


u) 


The corresponding rotated transform variables p, g, with the property pt + qn = sx — ty, 
are given by 
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= § COS a tsin a 
Dp + , 9 (12) 


q = SsiN ay — t COS a. 


We notice that s’ + ¢ = p? + g’ is invariant under the rotation, so that K* = k* + s°+ 


? =k’ + p’ + ¢’. It is clear that (10) is unchanged if A, B are understood to be functions 
of p, q rather than s, ¢. 

We now deduce the basic transform equation of our problem. Observe first that 
(10) permits us to write the transform of (4) as 


—iwuH = (pi, + gi, F iKi,;)-E, (13) 


where i; , i, , is , are unit vectors in the &, 7 and z directions. Similarly, we write the 
2 
>) 


)} as 


transform of ( 
(pi, + gi,  1Ki,;)-E = 0, (14) 
which according to (10) is equivalent to 
pA; + qA, — 1KA, = 0, (15) 
pB; + gB, + 1KB, = 0. 


Since conditions (5), (7) yield A; = B; , A, = B, , we see that B, = — A, , while 


A, = = (pA; + QA,). (16) 


The transform of (6), with [H,] computed from (13), then gives 


qA, = —1iKA,, (17) 
from which 
7 qA: . 
A, = -PtAc,, 418) 
ki +p 
Now calculate [H,] from (13). If we use (17), (18) to eliminate A, , and if we denote 
(twu/2) [H,] = A, the result is our basic transform equation: 
(k? + p)A = k°iK A, exp (—iKz). (19) 
Now simplify (19) by introducing the boundary condition (5). Let 
E,. = | dx | dy e ***‘“e(x, y, 2), 
a = | dx | dy e ****"e(x, y, 2) 
E.(8, t, 2)|,<0 = E.(0) 
and define the component of E. (0) in accordance with (11) as 
FE. (0) = F..(0) cosay — E.,(0) SIN @& . 
Then 
E,.(0) + E_,(0) = A; exp (—7Kz2p), (20) 
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where we see from (5) that 


E,,(0) = -|[ ax [ dy e*** “eo: 


and from (1) that 


TT. 11. | 
€or = (22 COS Ap — = a, sin a.) ° (21) 
. Oz 0: z=0 


v Oz OY 


Since II, differs only by a factor of k/4m from the free-space Green’s function of the 
three-dimensional wave equation, we readily deduce the y-transform (valid for 


te t | < k,) 


“= exp [—ik{(x — 20)? + (y — yo)? + 20}'7 + ty] 
vile “0 § Yo “O05 (22) 


= Te Hy [Kole — 2)" + 20)", 


where 
K, 7 (k? + r, 
and 


(ke + £)" ce = Hk. 


ae 
70 


of the right side of (22) is then obtained from the known integral representation 


The z-transform 


aes re 1 r’°*'? exp [w(z — x9) — (Ko + w’)'20] 
Hy” [Ko{(x — 2)” + 2}'’?] = --— [ So du, (28 

[K U q 5 ] mi Jy, - (Ki + w’)'” ) 
where | w, | < | Im K, |. Recalling that z, > 0, we conclude that the double transform of 
21), valid for Res > w , is 


s— Ww 


- } pee = 73 —-— le ( 4 °)! "e. 
E.,(0) = [ exp [—wre — (Ko + w)” Zo] 4,4, (24) 
Applying (20) and (24), we put (19) in the form 


= @& 


(ke + pA ~kpKel™ [ exp [—war, — (Ko + w’) “Zo] dw + KiKE_(0), (25) 


where | Wp | | Im K, |. Either of the two unknown functions A, Z_,(0) completely 
determines E. 

3. Transforms of the field components. The usual Wiener-Hopf techniques are 
not immediately applicable to (25), since the domains of regularity of the functions 
depend on the two complex variables s and ¢ simultaneously. In this section, we show 
how the difficulty may be overcome by restricting all operations to a suitable range of 
t. For such t, (25) is treated as if s alone were the variable. Certain representation theorems 
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for the Laplace transform are applied to carry out the function-theoretical argument, 
and the resulting knowledge of A yields the transforms of the various field components. 

Now all field components are required to be exponentially damped solutions of the 
three-dimensional wave equation. It follows that A = (iwu/2) [H,] is regular for Re 
(s + 7K,) > 0, while Z_,;(0) is regular for Re (¢K, — s) > 0. The respective half-planes 
of regularity depend on K, = (k’ + ¢°)'”’, but a discussion based on (22) shows that 
the same /-condition suffices for the regularity of both A and E_,(0), namely that ¢ 
be in the strip | Ret| < k.. 

The factorization K = (s + 7K,)'’? (s — iK,)'” then permits us to rewrite (25) as 


awott@ 


kpe'’’(p — ik)~'(s — iK,)'” | exp [—wa, — i(Ko + w’)'2z9](s — w)7' dw (26) 
Ji FA 8) 


= —(p + tk)(s + tKy) ‘A + tk*(p — ik)"'(s — iKo)'?E_,(0), 


where | w, | < | Im K, |, and where the condition | Ret | < k, is met by choosing Re 
t = 0. Then | Im K, | > k, , and we may take | w,) | < k, . Denote the left side of (26) 
by f(s) and the first and second terms on the right by g(s), h(s) respectively. The equation 


f(s) = g(s) + his), (26’) 


where f(s) is regular in the strip w. < Res < k, , and g(s), h(s) are regular in the over- 
lapping half-plane Re s > — k, , Res < k, , may be solved for the unknown functions 
g(s), h(s) by an application of the Wiener-Hopf technique. The procedure appears to 
be equivalent to that given formally by Harrington [6], but the justification may perhaps 
be of interest. 

We reason as follows. The fact that f(s) is regular in a strip suggests that it is there 
represented by a two-sided Laplace transform. We shall assume such a representation: 
the assumption will be justified on the basis of results we obtain for g(s) and h(s). Pro- 
visionally, then, we write 


f(s) = | e “*F(x) dz, (27) 


where 


aforia 


F(x) = = | e*f(¢) dg (| fo | < ke). 
2mt Jr. -ic 
Next observe that A(s) has the form 
h(s) = ik’ sec as '*E_.(0) + 8 *’*h,(s), 
where h, (3) is bounded for Re s < k, . The function 
s-*’*h,(s) = h(s) — ik’ see as '’*E_,(0) (28) 
then meets all conditions of a standard representation theorem [7] for the one-sided 


‘/ E_.(0) is a Laplace transform [namely the 


Laplace transform. Since ik’ sec ca 8 
1/2 > 97 2 ’ . 
and of zk” sec ay E_;,(0)], it follows 


transform of the convolution of the inverses of s 


that A(s) is. 
As to g(s), we notice that A and s"’”” A are transforms. To consider s A, write 


A WHY: * dfh,] 1s 4, Q) 
sA = lim / az ¢ dx. (29) 











1959] DIFFRACTION OF A DIPOLE FIELD 119 


An integration by parts, with an application of our edge condition (8) as e — 0, 
permits us to conclude that sA is a Laplace transform; an argument similar to that 
given for h(s) then shows that g(s) meets the conditions of the representation theorem 
[7] used above. 

Our result (from which the assumed representation (27) follows directly) is that 
g(s), h(s) are Laplace transforms of certain functions G(x), H(x) over the ranges (0, ~), 
(— , 0) respectively. Adding these representations to the now proved (27), we replace 


(26’) by 
x so) 0 
[ e "F(z) dx = [ e “G(x) dx + i] e "H(x) dz. (30) 
The standard uniqueness theorem [8] for two-sided transforms readily yields 
[ e “F(x) dx = g(s), (31a) 
[ e"F@) dx = hi. (31b) 


Our transform equation (26’) is therefore solved for the two unknown functions g(s), 
h(s). We remark that the familiar function-theoretic ingredients (factorization via 
Cauchy’s integral formula, the appeal to analytic continuation and to Liouville’s theorem) 
of the Wiener-Hopf techniques are implicit in the representation and uniqueness theorems 
we have used. At the same time, it appears that the implications and domain of appli- 
cability of the technique may be considerably enlarged by drawing on more general 
representation theorems, which are independent of the classical function-theoretic 
approach of Wiener and Hopf. We develop this point of view, and give applications to 
diffraction theory, in a forthcoming investigation. 

Expressions for the transforms of the various field components follow at once from 
either (31a) or (31b). Let us apply (31a): the substitution of our definitions [see (26), 


(26’)] of g(s), f(s) and 


;, SY utts ts , 
F(z) = | foe’ dé (| fo | < ke) 


in (3la) leads to 


A = ke'’*(s + iKy)'(p + ik)’ -J, (32) 
where 
pve**® (w COS a + tsin ao)(w — iK,)'” : +2 2\1/2 
J = | anne ——_———— exp [—wa — ta(Ko + w)“]dw (33 
(w — s)(w CoS a + tsin a, — tk) 7 - + ) ‘ ) 


and | w)| < | Im K, |. We then apply (19), (18), (17) in succession to find 


exp (tyo)(p — tk) 
ike — 1K)?” (34) 


II 


A; exp (—tkzo) 


pq exp (ty) J. (35) 





A, exp (—thza) = ik(p + ik)(s — ‘x. * 


xp (—ike,) = —B&XPAtyo)(s + 1K)!” = 
A, exp (—tkz,) = camp A i 
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The x — and y — components of A(z, t) follow from (34), (35): 


a exp (tyo) {(k? + p*) cos a + pq sin ao} 
A, exp (—1Kz) = ik(p + ik)\(s — iK,)'” 





J, (37) 


J iveeae es —exp (tyo) {(k° + p*) sin ao — pq Cos ao} = 
A, exp (—iKs,) = —SRARLIE + 2) OR ee — pe oeee) 7 (38) 
1, exp (—7Kz) = om (> + we — Kd 72 


In accordance with (10), the scattered electric field is given by 
l ia ate , — , 
e= 72-73 | | {A(s, t) exp (—7Kz)} exp [sx — i|2z|K — ty] dtds, (39) 
(292) Js J7 . 


where S denotes a vertical contour from s. — 7 © to s% +7 ~, with | | < k, , while 
T denotes a vertical contour from — 7 ~ to + 7 ~. We now substitute (34), (35), (36) 
in (39) to obtain explicit expressions for the field components e , e, , e, . Writing W for 


the vertical contour from w, — 7 © to Ww) + % ©, where | wo | < | Im K, |, and denoting 
(w cos a + ¢ sin ao) by — po , we introduce the complex integrals 
] an i Po 
sles ~ (2xi)?k ; I, I, @ + Wy pe + ik) 
slid de (40) 
exp [sx — wr — 1K |z| — i(Ko + w’)'’ “eo — tly — Yo) 
— sa =. ——2 | dw dt ds 


(s — w)f{(s — iK, )(w rs 1K,)}'"° 


and 


_ * a a a [i a a ; oo ~ .2\1 a om / = ) 
2 | | | SP “J (Ko = ie nae dw dt ds (41) 
(Qn) Is Ja day (s — w)(p + ik) (p + ik) ‘(s — iK,)(w + iK,)} 





in terms of which we express the electric field components as 


e: = I, + G: 9 ( $2) 
l 

e, = G ft: 
2 (43) 

¢ ; G ( { 1) 


where the subscripts of J, G denote partial derivatives, with the — and & derivatives 
defined by 


C a] oO . » 
= CUS Ga ~~ =" BT Ga « (45) 
0& Ox OY 

5] 0 0 . : 
—— = —— COB ds ~ =" SID Ao. (46) 
0&o OX» OYo 


The operational equivalents of 0/d£, 0/d& are accordingly p, po respectively. 
Applying (4), the components of the scattered magnetic field are 


h; = @. (47) 
Bos 

h, = = G; zto (48) 
Tay 

eo (49) 
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The z- and y-components of the field are linear combinations of the &- and 7-components: 


€, = €¢ COS Q + €, SIN a 


€, = —€; SIN a + €, COS a 
and 

h, = h; cosa + h, sin a 

h, = —h; sinay + h, cosa . 


4. Discussion of the field; edge behavior. In this section, we rviate the integrals 
I, G of (40), (41) to an integral 





i a j [ exp [sr—wao—iK | z |—1i(Ko+w’)'"2zo— t(y—yo)] 5 
F=ary I Je dy (s—w) {(s—iK.)(w+iKo)} — oe 
which may be calculated explicitly. It is in fact readily shown that @ = e"“"/kR + F 
is exactly the classical solution of Macdonald [5] for the field of a point source in the 
presence of a semi-infinite screen on which the field vanishes. Macdonald’s result 


® = I, ea - , (51) 
where 
cr ‘iat ; 
In =3 [ H(kR cosh u) du, (52) 
_— ; il H.(kS cosh ) du, (53) 
with 


R = [(x — a)? + (y — yw)” + (2 — &)*J'”, 
S = [(@ — 2)’ + (Y — yo) + & + 2)'), 
ur = sinh {(2/R)(rro)'” cos [(@ — ¢o)/2]}, 
us = sinh ' {(2/S)(rr.)'” cos [(@ + ¢)/2]}, 


and 


x=rcos¢, Lo = To COS do , z=rsing, Zo = To SIN do , 


will therefore facilitate our discussion of the present solution. 
Now consider the definitions (40), (41), (50) of J, G, F. We find 


3a’°G 2 
aE Oke = k°] (54) 
and 
0 ; rs] ‘ y — Be 
(2 + it)( 2 ao it) = k'F, (55) 
so that 
eee LE 2g 
T=F+G+ ik (2+ a, G. (56) 
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Applying (54) to e as given by (42), (43), (44) we have 
oS = (Je. + a ’ Rags ’ | oe 8 


The result F — 


= é 
given by (51), (52), (53), is then combined with (56) to give 


f 


e+e =)%. +4, %. +", 


( 
( 


where 
y 1 ‘6 Y Y 
Vv — G..¢ + G.. + ik (G, .¢¢. G te)» 
Y ] y v 
Vv = Giz + — (G, Eo a G, t), 
tk 
3) Y 1 Y Y 
v - _ ae + 7 Gr cc8. + —_ 3 
ik 


An explicit integral representation for ¥'"” now follows from (41), (59i). 


first that 


oe 
(2 os ik jw = —1 COS a az, 


aS K | 2 | 
———5 exp [st — wx, — 1K | 2z 
= oe P | 


where 


M = 





dw dt ds 


(Vol. XVII, No. 2 


(57) 


~*2/kR) + &, where ® is Macdonald’s fundamental solution as 


(58) 


(591) 


(591i) 


(59111) 


We see 


— (Ki + w)'72, — ty — y)) ———s 773" 
ames Y — ¥o)] 1 GK, (w + iKo| 


But we may integrate over W, 7, S successively to evaluate M. Consider 


M, =5— | (w + iK,y)'” exp [—way — tzo(Ko + w’)'”?] dw. 
“ml . uw 
Let to = To COS do , 20 = 1. SiN @ and apply Cauchy’s theorem to deform the path into 


iK, cos (do + in), with — © < 9» < o; the result is 


Ww = 
V. = l [ exp (—iKory cosh n)Ko sin (@ + tm) dy 
oe, ce (2i1K,)'”* cos [(¢o + in)/2] 
ied: OO eee — . . 
aimee -* cupaas | exp (—7iKoro cosh n) sin [(@ + t)/2] dn 
“ml J —a@ 


= —isin (¢/2)ro'” exp (—iKoro). 


Similarly, we deform S into s = — tk cos (@ + in) to obtain 


1 / , ee 
= — (s — iK,) 
271 J s 


* exp [sx — iK | z |] ds = 


M, 


i sin (@/2)r"'”? exp (—iKor), 
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where x = r cos ¢, z = rsin ¢. We then have 


os [ (M,)(M;) exp [—t(y — y»)] at 


= (rro)~'”? sin (¢/2) sin (¢,/2) il exp [—t(y — yo) — iKo(r + ro)] dé, 
r 


and, applying (23), 
M = —xk(r + ro)(rro)'” sin (6/2) sin (¢o/2)p 'H;” (kp), 
with p = [(y — yo)” + (r + r9)’]'”’. Solving the differential equation 


oM 
(2 + ix he = —12 COS a 


we find 
Vv =r’ sin 2){mik COS a SiN (d_/2) exp (—tkéo) 
rs) ** (r + 10) Hy” (ke’) . 
_ [- “Gy : 7 exp (iki) déo?, (60i) 
where 
z =r cos¢, z=rsing 
rs = [(x)’ + a)”, 15 = £5 COS a + m0 SiN a 
= [(y — ys)” ++ 7r6)’)”, ys, = —£ sin ao + no COS a « 


It is to be noticed that WV’ vanishes on the screen, as does %;,, , so that our boundary 


condition (5) is satisfied. 
Operational considerations now permit us to evaluate ¥, ¥“’. Let G, denote the 


image of G under our two-dimensional transform; the Sa equivalences 
= (—1/k)(Ki + w’)'?(p + po + iWG, , 
Vv? = (—1/k)(Ko + w’)'?q(p + po + ihG, , 
vw? = (i/k){(Ko + w’)(Ko + 8°)}'7(p + po + th)G, , 


follow from (59i), (59ii), (59iii). Then let ¥{" denote the image of ¥'"’; the equivalences 





(2 — q (1) 
ooo ee 
y® ie —i(Ko + + 8)" gan 
” p + ik * 
lead to evaluations 
g (1) | 
(2) asi ov , i 
y? = eit [ exp (ike’) “S| at (60ii) 
é a) 
(3) 1 . , ov , eee 
gy = eit [ exp (ike’) “S| a, (60iii) 
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where ¥'”? is given by (60i). It is clear from (60ii) that Y‘” is an even function of z and 
therefore continuous across the screen: since ®,,, vanishes on the screen, it follows 
that the jump condition (7) is satisfied. 

The behavior of Y‘” in the neighborhood of the edge of the screen follows from 
(601). As r — 0, the behavior is 


Vv? ~r'” sin /2)r(y), (611) 
where 
) oe gee eee ee ) 
vy) = rik cos Q& exp (— tk&,) ae, | (72)! 2 sin ($/, “= L at £ exp (iket?) dts, 
£0 v—@ p r=0 


. ° . . 2 ( 
The leading terms in the expansions of ¥’, ¥“ are: 


Vv? wr’ sin (6/2)v(y) tan ao (61 ii) 


3 —1/2 9),./ nytt: 
v"~-r cos (¢/2)v(y) sec a . (61111) 


For the net edge behavior of the electric field components, we see from (58) that we 
must add terms arising from the differentiation of ® to the ¥‘”, vV'’”, ¥ terms. The 
terms to be added are obtained from the expansion 


@ ~ Ar(rr,)'” sin (¢/2)(sin ¢o/2)(po) Hy” (kpo) 


(po = p|,-.), which follows from our expression for F [see (50)] upon our noticing that 








(2 25 ae \p _  M _ xr + 79) sin (/2) sin (¢o/2) H;” (kp) 
Ox Ox k (rr) ' p 
The explicit results for the leading terms are 
-1/2_° (i 2) , 
eé~r sin (6/2) — — I, (62) 
. ? \k 0&5 
Shas roth \c itana, dv ds 
e,~r sin (@ 2)( ——— ), (63) 
\ k 0&5 
-1/2 - SEC Qo 2) ia 
eé.~r?r COS ( 2)(§ - emg P (64) 
° ? \ ik 0&> 


where v(y) is given above under (61i). 
To complete our description of the field near the edge of the screen, we consider 
(47), (48), (49). The results for the magnetic field components as r — 0 are found 


to be 

he = (0), (65) 

: oe 2 ) Ov 7 

h, ~ (r)'” cos (6/2) -2i(e/p)'”” sec’ ao ap? (60) 

0&5 
staat 2 aim Ov - 
h, ~ (r)'” sin (6/2) -2(€/u)'”” tan ao sec a BE. (67) 
OSo 


It is evident from (66) that our edge condition (8) is satisfied. We remark also that 
the edge behavior given by (62), (63), (64) verifies that obtained by Toraldo di Francia 
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[1] in an approximate treatment of diffraction by a unidirectionally conducting small 
circular disc. The field behavior in the vicinity of a diffracting edge is, as anticipated, 


independent of the shape of the screen. 
5. Summary. The results of the foregoing analysis may be summarized in the form 


of a theorem: 
Theorem. Let a dipole field e, , hy derived from the Hertz vector 


(0,0,1,), with I, = (~“*/kR), 
be incident upon a screen 


z 0, —o <y< o, z=0 


IV 


which has infinite conductivity in the £-direction, where 


& = © COS ay — YSIN a (0 < a < 2/2), 


£ 
and is perfectly insulating in the direction » normal to é: 
n = XSiNa + YY COSa. 


If the resulting scattered field e, h is required to satisfy Maxwell’s equations, to be 
outgoing at infinity, and to meet the conditions (boundary condition, two jump con- 
ditions and an edge condition) of unidirectionality: 


BC: e =0, on the screen 

Jc: [h;] = 0, across the screen 

Pe [e,] = 0, across the screen 

EC: [h,] = 0, at the edge of the screen, 


and also to be integrable at the edge of the screen, then the total (vector) field 
&=e+e, KR=h+h, 


is given uniquely by 


6 = (%, +, %,, +, 6. +I, +¥%}, 


. JOH, oll, cave 
aC = iw}— , — —* + Dd‘, b**’?, 
On 0& 
where the components of the vectors &, 3C are in the &-, 7- and z-directions, and where 


® is Macdonald’s classical solution, vanishing on the screen, and given by 


m = fs — he ; 
with 
I. = : [ H<\?(kR cosh p) d 
ee ee 1 \A OSO wf) Gp, 
I; = s/ H{?(kS cosh yp) dp, 
= —BS 
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and 

R = [(x — 20)? + (y — yo)? + (2 — &)*]'”, 
[(z — 20)* + (y — yo)” + @ + 20)"}”, 
= sinh” (@/Rere)™ cos [(¢ — ¢o)/2]}, 
us = sinh {(2/S)(rro)'”” cos [@ + ¢)/2]}, 


2 
I 


bs 
ea 
| 


x =r cos¢, z=rsing, Lo = To COS Go , Zo = fo SiN do - 


The functions ¥"”, v'”’, Y™ are defined by 
vw" =r” sin @, 2){nik COS ao Sin (¢,/2) exp (—ikép) 


fo (2) , 
2 [Epp HE oxy cts a 
a A\i72 , exp (7k&) dtj/, 


@ (ro) 
é qa) | 
er ag ™ | exp (ikt’) il | dt’, 
-« On |gne? 
é (1) | 
yy? = “iad | exp (iké’) ~— | dt’, 
-« gag’ 
where p = [(y — yo)” + (r + 10)’]'”, and the functions b‘”, b“ are given by 
é 
b’” _ [ (F a + vy) | | dt’ 
o &=t’ 
] 
hb‘? = -[@ = + vy) dt’ 
lgmg? 
with F = — II, + ©. 


Some of the field components are singular as the edge of the screen is approached, 
but these singularities are of the physically admissible 0 (r~’””) type, as may be seen 
from Eqs. (62) through (67) above. It is especially to be noticed, among our results 
on edge behavior, that the y-component of the electric field vanishes along the edge 
of the screen; this follows either from (62), (63) or from (38). We conclude that a suitable 


edge condition for our problem is 


at the edge of the screen, exactly as for a perfectly conducting screen of the same 
geometry. The use of the corresponding condition, namely that the field component 
tangential to the rim of the screen be required to vanish, is therefore suggested for 
further investigations of unidirectional screens. This condition is met by the approximate 
solution of Toraldo di Francia [1]. 

We remark that we have given the highest order terms in r, in discussing edge 
behavior: it is assumed that 1, is finite. If however we let r, — ©, the character of 
the results for field behavior near the edge simplifies further. This type of approximation 
was given by Senior [4] in discussing edge behavior for the incidence of a dipole field 


on a perfectly conducting screen. 
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We conclude by verifying our theorem. Let us write Maxwell’s equations in the form 
(Ag, + k*)3 = 0, 
V-x = 0, 
weE = V X KX. 


Notice now that ¥"” is a wave function, and that ¥“’, ¥ are therefore wave functions. 
It is then clear that b'”, b‘*’ are wave functions. Remarking [see (55), (59)] that 





ov? ov? oy . 
ot + an + dz = k F,, 


we find that (&, 3C) is an electromagnetic field. The boundary condition 
Ce = —log 


is satisfied; and the first jump condition is obviously met, since h; = 0. We have pointed 
out in giving VY ~’ that it is an even function of z; and it is clear that en is continuous 
across the screen, since ®,,, is. The edge condition [h,] = 0 at the edge of the screen 
follows from the edge behavior of hy. 

The author wishes to express his indebtedness and gratitude, for suggesting the 
problem here treated as well as for many valuable discussions, to Professor 8. N. Karp. 
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BOOK REVIEWS 


Switching circuits and logical design. By Samuel H. Caldwell. John Wiley & Sons, Inc., 
New York, 1958. xviii + 686 pp. $14.00. 


This book is needed by every serious student of the logical design of switching circuits. The author 
has performed a much needed service by gathering most of the significant material on this subject into 
a single volume, reworking it into a coherent and well-written textbook. 

The jacket blurb about Professor Caldwell is unnecessarily modest about his qualifications to have 
written this book. The two most important advances in the field of logical design were papers by Shan- 
non and by Huffman, and each of these papers was written as an M.I.T. thesis for which Caldwell was 
a thesis adviser. Understandably he gives much space in his book to exposition of the ideas of these and 
later publications by Shannon and Huffman, and of significant unpublished material by these two men. 
The book also contains unpublished material original with Caldwell, and various items due to other 
M.I.T. men. The author takes so much advantage of the excellent resources of M.I.T. that he errs some- 
what on the side of provincialism, tending to ignore some of the important papers written elsewhere, 
such as Gilbert’s work on frontal switching functions, and Muller’s work on complexity in electronic 
switching circuits. He completely omits all work on logical design done outside the English-speaking 
world. The books by Gavrilov, Higonnet and Grea, and Plechl are never mentioned, and the papers by 
Lunts, Povarov, Shestakov, Cardot, and others are similarly ignored. 

But many useful methods are treated, some of which had not been published or widely known before. 
This book gives a much more full and thorough coverage of the subject of logical design than any of the 
other books on the field, including those mentioned above, and the English language ones by Keister, 
Ritchie, and Washburn, by the staff of the computation laboratory, by Richards (all cited in the book 
under review ), by Phister (Logical design of digital computers, Wiley, 1958), and by Culbertson (Mathe- 
matics and logic for digital devices, Van Nostrand, 1958). 

This book gives an intuitive rather than a formal approach to its subject, and depends on many well 
explained illustrations and many tricky exercises to appeal to the ingenious puzzle-solver type of student, 
rather than attempting to prove theorems or set up formidable mathematics. The book frequently goes 
out to the limits of what is known in a given part of the subject, and explicitly or implicitly states many 
problems which are of current research interest. 

A strong point of this book is the emphasis that it gives to asynchronous circuit design. The other 
recent books on logical design cover this topic very incompletely, despite its practical importance, and 
the many interesting theoretical problems which arise in dealing with it. 

There are two minor difficulties which may stand in the way of the wide use of this book as a text- 
book. The first of these is that most of it (except Chapters 9, 14, and 15) is written about relay circuits, 
but that most electrical engineers who will be teachers or students of courses in this subject are preju- 
diced against relays, as compared to the much faster, more modern, and more glamorous electronic 
switching elements. From both the puzzle solving and the mathematical points of view, relay circuits 
have many interesting properties which are not to be found in the electronic circuits, although much of 
the skill or knowledge gleaned in either area is applicable in the other. Perhaps the advent of the cryo- 
tron (a recent electronic switching element which is functionally equivalent to a relay) will make this 
prejudice obsolete. 

The second of these difficulties is the close knit way in which the book is integrated. This would be 
an advantage when teaching a two-semester course just like the one Caldwell teaches. However, he has 
not only failed to supply asterisks to show which topics may be omitted without loss of continuity, but 
has interwoven these topics closely with the absolutely essential ones. Anyone wishing to teach a shorter 
course or to make room in a two-semester course for other topics not treated in this book, such as com- 
puter programming, system design, or arithmetic units, will have a major editing job in store for himself. 

But neither of these difficulties as a textbook will impair the lasting value of this book as a refer- 
ence book, explaining in clear and heuristic language the main methods of logical design or its current 
value as a source of orientation about the frontiers of the field, with its frank comments on many of 


the kinds of problems we do not know how to solve. 
Epwarp F. Moore 


(Continued on p. 146) 
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GENERAL SOLUTION OF REYNOLDS EQUATION FOR A JOURNAL 
BEARING OF FINITE WIDTH* 


BY 
L. N. TAO 
Illinois Institute of Technology 


1. Introduction. The governing equation of hydrodynamic lubrication, the Reynolds 
equation, was established in 1886 [1]. Since then many works on the solutions of the 
equation have been published. For plane slider bearings, where the film thickness of the 
lubricating fluid is a linear function of the coordinate variable, solutions for both infinite 
and finite widths have been obtained. Michell [2] and Muskat et al. [3] gave their solu- 
tions for finite width bearing in terms of Bessel functions. In journal bearings where 
the film thickness is not a linear function of the variables, Sommerfeld [4] showed that 
for the case of infinitely long bearings the generalized Reynolds equation, a partial 
differential equation, is reduced to an ordinary differential equation. The solution of the 
problem can be immediately derived by integrations. However, for journal bearings 
of finite length, all of these investigators mentioned in their papers that they did not 
succeed in obtaining the general solution for the finite bearing problem. Michell even 
stated that the exact solution of Reynolds equation is only possible if the viscosity is 
constant and the film thickness is a linear function of the variables. Following this line 
of argument, the subsequent works on the solutions of finite journal bearings** are 
approximations by means of a wide variety of mathematical techniques. The aim of this 
paper is to establish an exact and complete solution of the Reynolds equation for finite 
journal bearings. 

In Sec. 2 of this paper the Reynolds equation, by suitable transformations of variables, 
is changed to a comparatively unknown equation, Heun’s equation. In Sec. 3 Heun’s 
equation and its function related to the problems are discussed. The convergency of 
the Heun function is also studied. Then in Sec. 4 the pressure distributions, as well as 
the eigenvalues and their eigenfunctions, are established. Finally, the load carrying 
capacity and the coefficient of friction of the bearing are determine in Secs. 5 and 6 
respectively. 

2. Governing equation. The governing equation of hydrodynamic lubrication [1] is 


9 Ip | ) di 
: G =| = G op] = Gul —, (2.1) 
Ox Our Oz Oz Ox 


where p, h and yu are the pressure, film thickness and viscosity of the lubricating fluid 
respectively, and U’ is the relative velocity of bearing surfaces. In journal bearings x 
and z are taken along the circumferential and axial directions respectively, and the 


film thickness may be written as 


¢ 
h c(l + 7 cos 6) and 7=-, (2.2 
Cc 


**Received March 24, 1958. 
*See bibliography cited in [5] and a recent work [6]. 
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where c is the radial clearance and 7 the attitude, a ratio of eccentricity e to radial 


clearance c. 
Introducing dimensionless quantities 


©, H=-=1+ncos0, P=, (2.3) 


x - Z 
g==, Z=? 
i r 


with r being the radius of the bearing, Equation (2.1) becomes 





an + 7 cos 6)* — Fl +; +=; at + 7 cos 6)" oF | = a a (2.4) 
The boundary conditions of the problem are 
P(—12, Z) = P(x, Z), (2.5) 
ee 4. 3 
een) = ee (2.6) 
P(6, —L/2) = 0, (2.7) 
P(@, L/2) = 0, (2.8) 


where L = I/d and / and d are the length and diameter of the bearing. Conditions (2.5) 
and (2.6) assure the continuity of the pressure and its first derivative along its circumfer- 
ential direction, while the last two conditions, (2.7) and (2.8), state that the pressure 
at free ends is atmospheric. 

In order to remove the non-homogeneous term in Eq. (2.4), we let 


P(6, Z) = ¢(0, Z) + £(8), (2.9) 
where £(@) satisfies the equation of two dimensional journal bearings of infinite length, i.e. 
Q [a + 7 cos 6)° A = a ae (2.10) 
and its boundary conditions 
&(—m) = (x), 
(2.11) 


dé(—7) dé(r 


dé  do- 
The solution of Eq. (2.10) with the given boundary conditions (2.11) is [4] 


ay 6vrn (2 + 7 cos 6) sin 16 
(@¢) ==> —5- ———$—$—— (2.12 
£0) Uc'(2 + n) (1 + 7 cos 6)° ) 


The other function ¢(@, Z) satisfies 


0 
3 la + 1 co 02 | + FL + 9 cos 0) ot) 0 (2.13) 


and the boundary conditions 


t(—7, Z) = ¢(x, Z), (2.14) 
d¢( —7, Z) ee d¢(x, Z) (2.15) 





? 


00 06 
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(0, — L/2) = —&6), (2.16) 
¢(6, L/2) = —é&(@). (2.17) 
Separation of variables 
(0, Z) = $(6)-¥(Z) (2.18) 
gives 
d@ ___3nsin@ dg - 
de 1+ 7 cos 6dé reas (2.19) 
and 
d’y a 
yi ~ \¥ = 0. (2.20) 


From Sturm-Liouville theory [7], it can be shown that Eq. (2.19) with the given boundary 
conditions possesses a sequence of solutions corresponding to a sequence of real, non- 


negative eigenvalues X, . 
The solution of Eq. (2.20) is then 


4, ~ 0, yy, = A, cosh (A,)'?Z + B; sinh (A,)'?Z (2.21) 
X=0,  =AZ+B. 
Furthermore by the substitution of 


(2.22) 


Eq. (2.19) becomes 


a= fo "Et [we (ergo Lt] 
~~ ’ oa (2.23) 


This is known as Heun’s equation [8]. 
3. Heun’s equation. Heun’s equation 
2iz2 — 1)\(z — ay’ + [atBt+ D2’ — fatBt+1—é6+aly + Dj}z (3.1) 
+ ayly’ + oB(z — gy = 0, 


with a, g, a, 8, y and 6 being arbitrary constants, is closely related to the algebraic form 
of Lame’s equation (if y = 6 = 3, a + 8 = }). Clearly witha = g = lora=q=0, 
Eq. (3.1) is reduced to the hypergeometric equation. 

With regular singularities at 0, 1, a and ~, the Riemann’s P-function of Eq. (3.1) is 


0 l a © 





Ps 0 0 0 « ef (3.2) 
lli—-y 1-8 v¥+85—-@+ A) B J) 


Therefore the solutions of Eq. (3.1) for y being non-integer are 


y=1+ Die” (3.3) 
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and 
y=2 [1 + > ax | (3.4) 
The first solution is 
y=1+ > e,2" = F(a, q; a, B, y, 632), (3.5) 
1 
where F(a, 9; a, 8, y, 6; 2) is known as the Heun function with 
ayc, = aBq 
Qa*y(y + leo = (aBg)? + {((a+B+1 — 8 + (y + dajaBq — ayaB 
and 
a(n + Ivy t+ neu = [(atB+tn— n+ any + 6+2n —- 1) (3.6) 
0.0 
+ aBgle, — [(n — DfatB+n— 1) + aBle,-, . 


The other solution is 


~ 


y = gf + az | =z "F(a,g,;a-y+1,8—7+1,2 —7, 62) (3.7) 


with 


a oe ae ) — sae ap owen 
q, = \g + — l(a + B) + da 1) + (1 a 





Now to establish the convergence of the function (3.5) we consider the series (the 
Heun function with argument z = 1), 


} Be (3.8) 
n 0 
and the recurrence relation 


and ; ae | (3.9) 





where 
Qo — afbq, ' 
P, = : (iy +n — 1), 
n 
i \ (3.10) 
Q,=—a3alatBt+n— bdnt+any+6+n-—-1) + aBq], | 
n 
Rk, = 4 [(n — l(a+B+n — 1) + af]. | 
n J 
Let 
: 3 : 
| Fe (3.11) 
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and 
P = im?,, Q = lim Q, , R = lim R, 
then the recurrence relation (3.9) becomes 


, —— 
N. = 5p WW (3.12) 


By repeating substitutions, Eq. (3.12) may be written as an infinite continued fraction, 





r R, Pav Rasr _PrsoRnss | 
N, =- 0. ~ :.. ~~ coe, (3.13) 


It can be shown [9] that this continued fraction is convergent when the roots (p, , p2) 


of the equation 
Pi? —Qe+R=0 (3.14) 
are | p, | # | p» |, or p; = p, and it is divergent when p, # p, but | p, | = | p2 |. Also in 
the convergent case with | p, | < | p. | or p, = p, we have 
lim NV, — p, . (3.15) 


In the present series P = a, Q = 1 + a, and R = 1, the roots of Eq. (3.14) are 1 
and 1/a. This, in turn, shows that the Heun function (3.5) is convergent for all | z| < | a] 
if|a|> 1orforall|z| < 1lif|a|<1lora=1. 

4. Pressure distribution. By comparison of Eqs. (2.23) and (3.1) we get 


1 + l 1 2 
a=, gat! « f= 5820+ 0)" 
" sa * (4.1) 
l 1 
Ss: 6 = —~* 
_s 2 
Since a > 1 when n < 1, the Heun function of the problem is convergent for all 6. And 
a = 1 when 7 = 1, the Heun function (or the hypergeometric function) is divergent. 
However, 7 = | implies a contact of bearing surfaces. This is no longer a case of hydro- 


dynamic lubrication. Hence the complete solution of Eq. (2.13) is 


£ = bobo + D> [A; cosh (A,)'7Z + B; sinh (A,)'?Z](C,F; + D-H), (4.2) 


where 
i r(! tn Ltn. 3+(9+4)'? 3-94+4)'7 1 1, 
| an’ fy’ 2 , 2 79°97") 
- (2 +n 1+(0+4,)'" 4-(94+ 4," 3 1 u) 
“ur q; ; scone es eg ow OE. 
é' (2a *™? 2 ?, 2 19°92 
7 13 1+ (3 - )] (4.3) 
a (; s.) E Qn \4 ri} |) 
Wo A,Z + B,, 


do Co | (1 + n cos 0) * dé+ D, 


and A, , B; , C; and D, are constants. 
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Conditions (2.16) and (2.17) indicate that ¢ is an even function of Z, so 


F, and H, are, therefore, even and odd functions of @ respectively. This gives C. = 0. 


From conditions (2.16), (2.17), (2.7) and (2.8) we then get 
C, = Dy, = 0. 
Also conditions (2.14) and (2.15) require 


DH f(r) = 0, 
Since D; ~ 0, 


Now the pressure distribution with A, to be determined is 
P(6, Z) = H*(0) + >> A; cosh (A,)'?Z-H(A,; ; 8), 


where 





6vrn(2 + 7 cos 8) sin 6 
* = => —— <« a ~ eet ae > 
H™(0) = &(8) Uc'(2 + n°)(1 + 7 cos 6)” 


and the eigenvalues \; are the roots of 


H(\; 7) = 0. 
H,(@) are orthogonal with respect to the weight function 
w(6) = (1 + 7 cos 6)" 


over the interval (— 7, 7); that is 


ar 


| w(6)H,(0)H,(0) d@=0 for i¥j 
hence 


[ w(0)H*(6)H,(6) dé 


® 





A; =-— 


cosh [(A,)'L /2] / w(0)H(0) do 


(4.4) 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


(4.9) 


5. Load carrying capacity. Relating the pressure distribution to the load carrying 
capacity W of the bearing, W is resolved into two components along and perpendicular 


to the line of centers, 
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L/2 , 
W, = pU’r’ / / P cos 0 dé dZ, 
L/2 -* (5.1) 


L ea 
a pur’ | / Psin 6 do dZ. 
—L/2 -* 


Due to the fact that P is an odd function of 6, we bave 


W, =0 (5.2) 
and 
W, L2avnl. r 1/2 
Us? = Ut pl — 4 + La 2 sinh [()'7L/2) } “HA. jw du. (53) 


This result indicates that the line of centers is perpendicular to the load line. The same 
result has been obtained for the infinitely long journal bearing. With A; always negative 
in the range 0 < @ < z, the last term of Eq. (5.3) is therefore always negative. This 
implies that the load capacity of finite length bearing is less than that of infinite length 
bearing. This result has also been found experimentally by electric analogy [10] and 
theoretically by approximation methods [3, 6]. 

6. Coefficient of friction. The frictional force acting on a journal bearing may be 
expressed by 


+0 24 , 
=f" i i Free penn ted (1+ ncos@) 5, |deaz. (6.1) 


By Eq. (5.3) we obtain 


fa sn 4 Bw, . (6.2) 


The coefficient of friction is defined as the ratio of fractional force to the total load W. 
And W, = 0, we have 


J 2rul ae F cn 
* W cel— n) ?W, + 2r (6.3) 


BIBLIOGRAPHY 


1, O. Reynolds, ‘“‘On the theory of lubrication,’ Phil. Trans. Roy. Soc. London, Ser. A, 177, 157-234 
(1886) 

A. G. M. Michell, “‘The lubrication of plane surface,’’ Z. Math. u. Phys. 52, 123-137 (1905) 
. (a) M. Muskat, F. Morgan and M. W. Meres, ‘‘7'he lubrication of plane sliders of finite width,” J. 
Appl. Phy. 11, 208-219 (1940); (b) M. Muskat and F. Morgan, ‘‘The theory of the thick film lubrication 
of a complete journal bearing of finite length,”’ J. Appl. Phy. 9, 393-409 (1938) 
4. A. Sommerfeld, ‘‘Zur Hydrodynamischen Theorie der Schmiermittel Reibung,’’ Z. Math. u. Phys. 50, 
97-153 (1904) 

5. M. C. Shaw and E. F. Macks, “Analysis and lubrication of bearings,’’ McGraw-Hill Book Co., New 
York, 1949, p. 242 

6. M. J. Jacobson, A. Charnes and E. Saibel, ‘‘The complete journal bearing with circumferential oil 
inlet,’’ Trans. ASME 77, 1179-1183 (1955) 

7. E. L. Ince, ‘‘Ordinary differential equations,’’ Dover Publ., New York, 1944, Chap. 10 


to 


Ww 








136 L. N. TAO [Vol. XVII, No. 2 


8. (a) E. T. Whittaker and G. N. Watson, ‘A course of modern analysis,’’ 4th ed., Cambridge Univ. 
Press, London, 1952, pp. 576-577; (b) A. Erdelyi (Ed.), ‘““Higher transcendental functions,’’ vol. III, 
McGraw-Hill Book Co., New York, 1955, pp. 57-62 

9. O. Perron, ‘‘Die Lehre von den Kettenbruchen,’’ Chelsea Publ. Co., New York, 1950, Kap. 7 

10. A. Kingsbury, “On problems on the theory of fluid film lubrication with an experimental method of 
solution,’ Trans. ASME 53, 59-75 (1931) 














137 


SHAPE-PRESERVING SOLUTIONS OF THE TIME-DEPENDENT 
DIFFUSION EQUATION* 


BY 
FRANK 8S. HAM 
General Electric Research Laboratory 
Schenectady, New York 


Abstract. Exact solutions to the time-dependent diffusion equation are exhibited 
which correspond to the diffusion-limited growth of ellipsoidal precipitate particles 
with constant shape and dimensions proportional to the square root of the time. The 
asymmetry of the diffusion field in these solutions is consistent with the preservation 
of the particle’s shape during growth even if the diffusivity is anisotropic. Limiting 
cases for simpler geometries are derived and shown to be in agreement with previously 
known results for radially symmetric particles and isotropic diffusion. Similar solutions 
for hyperboloidal surfaces are exhibited and generalizations are considered analogous 
to those discussed by Danckwerts for one-dimensional diffusion. 

I. Introduction. Exact solutions of the time-dependent diffusion equation are known 
[1-5] which correspond to the growth of radially symmetric particles in one, two, and 
three dimensions with radius proportional to the square root of the time. The initial 
density of diffusing material is uniform and the density at the surface of the particle 
is specified. The rate of growth is directly proportional to the concentration gradient 
at the surface. It is the purpose of this paper to show that these solutions are special 
eases of a more general class of shape-preserving solutions appropriate to ellipsoidal 
and hyperboloidal surfaces. The generalized boundary condition for the rate of growth 
is that the velocity of motion of the surface along its normal is proportional to the normal 
component of the concentration gradient. These more general solutions are of interest 
in connection with studies of diffusion-limited phase transformations for which the 
interface between the phases is not radially symmetric or for which the diffusivity is 
anisotropic, or both. 

II. General solutions. We shall examine only those solutions of the isotropic time- 
dependent diffusion equation 


2 ns DG, (1) 
ot 


which are functions of the reduced variables 


u = xt" A 
~ yt”, 
w = 2t'” (2) 


where x, y, z are Cartesian coordinates of position and ¢ is the time. In (1) pis the density 
of diffusing material and D an appropriate diffusivity, assumed for the present to be a 
scalar. Introducing the vector s = (u, v, w) and the gradient operator V, = (0/du 
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0/dv, 0/dw), we find that (1) becomes 


—}s-V.p = DVip. (3) 
We now introduce general ellipsoidal coordinates according to the relations [6] 
u = ((f — a y(& — a’ \(& — a’)/a'(a’ — b’)]'” 
v= [(& — BYE — bY — b°)/b°(b* — a’)]'” 
w = £itE,/ab (4) 


with £, >a >£& >b> & > 0. We find that in these coordinates (3) becomes 


1 mt... oO (Asbyle 2) , 
Z bs h? of ~ hyhehg > 0&; h: 0&; ’ (o) 


where 


h; = (& — ud) -— 8.) /(E -— aye — b*)]'”, (6) 


the variables £, , & , &; being placed in cyclic order. 


We shall be concerned with those solutions of (5) which are functions of one of the 
’s only, £; . Because h,h.h3/h; equals 
fé) = ((E — aE — b°)]'” (7) 
times a function of £;,, and &;,. , under these restrictions (5) becomes 
dp d dp 
—= Ef é,.) = = D= | fe) = |- (8) 
2 ENE) oe de, LI? ae, 
The general solution of (8) is 
r** exp [—f/4D] 
é) = c| Et dt, (9) 
he : f(t) 


where c and d are arbitrary constants. 

The surface —£, = & > a > b > O corresponds to the growing ellipsoid in the original 
xz, y, 2 coordinate system with principal axes of length 2¢'”* (¢ — a*)'”*, 2t'”? (& — b?)'”?, 
2t'’"¢, along the z, y, z axes respectively. Hence the solution of (1) which equals p, on 


this surface and (p. + p,) at infinity is from (9) 


p(é:) = p. + pof{l — [FE)/F(E)]}, (10) 
where 
vp) . [ exp l—#/4D] at 
Fé) = [ (0) (11) 


with f(t) defined in (7). At t = 0 this corresponds to the uniform density p = p, + p, . 
The rate at which the ellipsoid grows is determined by & , its shape by the factors 
[1 —(a/é,)*]'”, [1 —(b/&)7]'”*. To determine & we shall require that the instantaneous 
velocity v in the z, y, z coordinate system of a point (uo , % , Wo) on the surface §; = & 
shall have a component along the normal é, to the surface which is proportional to the 
normal component of the gradient of p: 

(p. — p,(v-&:) = D(E;-Vp). (12) 


But from (2) we find 
r = (zr, y,2) = t's, (13) 
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v= i°"s, (14) 
Ve = t'*V.». (15) 
Moreover, if —; is a unit vector we can show from (4) that 
(s-&;) = &/h; . (16) 
Hence (12) becomes 
3(p. — p.d& = D(dp/dé;) (17) 
to be satisfied for —; = & = & . From (10), therefore, & must satisfy 
2(Pe — Ps&ol (Eo) = — DpoF(é), (18) 


where F’(é) = (d/d&) F(é). This equation has a solution if pp < (p, — p,). 

Physically (12) corresponds to the requirement that the flow of diffusing material 
to the surface be equal to the rate at which the amount of the new phase must be added 
locally if the surface is to advance with a velocity v. In (12) p, represents the density 
of the diffusing material in the precipitate phase. In exhibiting (10) and (18), therefore, 
we have demonstrated the existence of a solution to (1) which satisfies this condition, 
with v given by (14), at every point of an ellipsoidal surface. The shape of the diffusion 
field for a density that is uniform at ¢ = 0 and constant at the surface during growth 
is therefore consistent with the growth of an ellipsoid of arbitrary but constant shape 
and dimensions proportional to t'’*.* 

Another interesting solution is obtained from (9) and (17) if we use é; = & . Surfaces 
of constant & (0 < & <b < a) are hyperboloids of two sheets, with &, = 0 being the 


u, v plane. Therefore 


plEs) = po + (p. — po) [GEs)/G(Eo)] os 
with 
ae) = [SP - (20) 


is a solution of (1) which equals p) on the xz, y plane and p, on the hyperboloid 


= £ <b<a, or 


wr 


9 9 


_t > ieee ae 

e-a°b-a a~~* @1) 
This hyperboloid has a trace on the x, z plane which is a hyperbola with foci at y = 0, 
z = + al” and as asymptotes the lines y = 0, x = + [(a° — &§)/&]'’z. On the y, z 


« 


*This analysis is not appropriate for the description of the growth of finite plates or rods of zero 
thickness (the limiting ellipsoidal forms), for the velocity of growth of the finite dimensions of the par- 
ticle is found from (18) to be infinite in this limit. However, our ellipsoidal solutions reduce in this limit 
to those for planes or cylinders given in Sec. III. Thus the length or breadth of the particle becomes 
infinite in any finite time, and the thickness grows at a finite rate given by the simpler solutions for one 
or two dimensions. Such infinite growth velocities obviously make the solutions inadequate representa- 
tions of any physical situation involving a finite three-dimensional particle, but this is not surprising 
since we do not expect the boundary conditions used in our analysis to be reasonable when any dimension 
of the particle is smaller than a few atomic diameters. 
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plane the foci are at x = 0, z = + bt'’”’ and the asymptotes x = 0, y = + [(b® — &)/ 


eats 

If we impose the boundary condition (12) or (17) on the surface , = & , we find that 

£ must satisfy the equation 

2(p. — p.)f&oG(éo) = D(p, — po)G’ (Eo), (22) 
where G’ = (dG/dé). Therefore (19) and (22) correspond to the diffusion-limited transfer 
of material between one branch of a dissolving hyperboloidal surface of two sheets and 
the plane containing the conjugate axes of the hyperboloid. Thus if the density of 
diffusing material on this plane is maintained at a constant value p, by an external agency, 
the hyperboloid retains its shape as it dissolves, preserving its asymptotes and eccen- 
tricity while the distance between its foci and the plane increases in proportion to the 
square “oot of the time. 

A second type of solution of (1) involving hyperboloidal surfaces corresponds to 
precipitation on such a surface instead of its dissolution as above. In this case we replace 
t~'”? in (2) by (r — t)~’””. The sign of the left-hand side of (3), (5), and (8) is corre- 
spondingly changed, as is that of the exponent in (9). The integrals corresponding to 
(10) and (11) diverge because of the positive sign of the exponential, so that we can 
not use a solution of this type to describe the dissolution of an ellipsoidal particle, at 
least not with the boundary condition of a given density of diffusing material at infinity. 
However, there is no such difficulty with (19) and (20), so that with a positive sign 
in the exponential this solution corresponds to transfer of material between the zx, y 
plane and the hyperboloid given by (21) with ¢ replaced by (r — t). The shape of the 
hyperboloid is preserved, but at time ¢ = 7 the surface has become a cone with its apex 
touching the x, y plane. The solution is valid only for t < r. 

III. Limiting forms of the general solutions. Morse and Feshbach [6] have listed 
prescriptions whereby one can pass from ellipsoidal coordinates to most of the more 
familiar coordinate systems by appropriate stretching, compressing, and translating. 
If these or similar prescriptions are applied to the solutions (10) and (19) and the 
boundary conditions (18) and (22) as well as the coordinates, solutions to (1) appropriate 
to these simpler geometries are obtained. These agree for radially symmetric particles 
with the results of Rieck [1], Frank [2], and Zener [3]. 

As an example, if we set & = a’ + wu”, & = b° + v”, and &, = w’ in (4), place b = 
a sin g, and then let a become infinite, we obtain u = wu’, v = v’, w = w’, or rectangular 
coordinates. Making these substitutions in (10), setting £ = a’ + wo’ , and taking the 

h 


limit, we obtain an expression of the same form but with F(é,) replaced by 


a « 


F(us) = | exp [—u°/4D] du. (23) 


This replacement in (18) also gives us the correct boundary condition. This is the usual 
solution [5] in one-dimensional diffusion for precipitation from a semi-infinite medium 


> 


onto the plane x = w/t’ 
Similar manipulations of (19) and (22) yield a one-dimensional result for the diffusion 


of matter between the plane z = w/t'’” and the plane z = 0. Finally, if we start from the 
I I y; 


second type of solution for hyperboloidal surfaces discussed above, we get a similar 


result for the planes z = wi(r — t)'”” and z = 0, all exponentials in the expressions 


corresponding to (19) and (22) having positive arguments in this case. 
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The results obtained for the various coordinate systems and surfaces as limiting 
forms of the general ellipsoid can all be put in the form of (10) and (18) if the function 
f(t) appearing in F'(£) is defined appropriately. The proper definitions are listed in Table 
I, along with values of a, b, and c which when substituted into the equation 

2 2 2 
x y Z 
 momeaaaes. ~ oe ae ee. § 24 
g-F*F-F i-¢ (24) 
define the surface over which £, and p(,) are constant. We restrict ~, to values greater 
than a, b and c when these numbers are finite. An infinite value for a, b, or c simply 
deletes the corresponding term from (24). 


TABLE I. 








Ellipsoidal solutions in special coordinate systems. 





Coordinate System Surface Parameters* f(t)t 

a b c 
Ellipsoidal a b 0 [(t? — a?)(t? — b?)]}*/2 
Rectangular 0 co co 1 
Oblate spheroidal a 0 0 t[t? — a®ji/2 
Prolate spheroidal a a 0 [t? — a?] 
Elliptic cylinder a 0 oo [t? — a®}t/2 
Circular cylinder 0 0 oo t 
Spherical 0 oO 0 {2 





*See Eq. (24) of text. 
tSee Eqs. (10), (11), and (18) of text. 


TABLE IT. 


Hyperboloidal solutions in special coordinate systems. 








Coordinate System Surface Parameters* f(t)t 

a b 
Ellipsoidal a b [((a? — #2)(b? — ¢2)]¥2 
tectangular ro) ro) 1 
Prolate spheroidal a a {a? — ¢?] 
Elliptic cylinder oc b [b? — ¢2]1/2 





*See Eq. (21) of text. 
+See Eq. (19), (20), and (22) of text. 
1 


The hyperboloidal solutions (19) have interesting and distinct limiting forms only 
for the rectangular, prolate spheroidal, and elliptic cylinder coordinates. These solutions 
and the boundary condition all have the form of (19) and (22) with f(¢) in (20) appropri- 
ately defined. The proper definitions are listed in, Table II together with values of a 
and b which when substituted into (21) (with & replacing &) define the surface over 
which £, and p(£,) are constant. We restrict £, to values between zero and the lesser of 
a and b. 

If the argument of the exponential in (20) is made positive and ¢ replaced by (7 — ¢) 
in (21), then (19), (21), (22) and Table II give the hyperboloidal solutions for which 
the scale of the geometry is shrinking instead of expanding. 
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Other types of coordinate systems may be obtained from the ellipsoidal coordinates, 
but either the limiting surface of interest is identical with one we have already listed 
(for conical coordinates, for example, the interesting surfaces are spheres) or the limit 
of (10) or (19) is indeterminate (for paraboloidal and similar systems). In the latter case 
we can show directly that for these systems (5) does not separate in the sense that if p 
is a function of only one coordinate the other coordinates disappear from the equation. 

Evaluation of the integrals appearing in our solutions can not be done simply, 
although those for rectangular and spherical coordinates can be expressed in terms of 
tabulated error functions. Frank [2] has given tables for the evaluation of (11) for 
spherical and circular cylinder coordinates. Since ¢t? > a’ > 6b’ in integrals of type (11), 
the more complicated integrals can be evaluated by expanding the factor multiplying 
the exponential in a power series in (1/t) and tabulating integrals of the form f%t" 
exp (—?°)dt. In integrals of type (20), t? < b” < a’, so that use of a power series in ¢ is 
appropriate. Also, if £/D* is sufficiently small, the exponential factor in the integrand 
can be replaced by unity. The resulting integrals are then integrable in terms of simple 
functions, the ellipsoidal result becoming an elliptic integral. This approximate pro- 
cedure is not satisfactory for cylindrical and rectangular coordinates, since (11) diverges 
in these cases if the exponential is omitted. 

In numerical evaluations of these solutions it is convenient to note that the shape of 
the surface on which the boundary condition is satisfied is fixed by the parameters 
a’/t? and b’/£ . With these chosen, the boundary condition (18) or (22) is then an 
equation for & , which determines the rate of growth of the surface. 

IV. Generalizations. It is possible to generalize our solutions to considerably more 
complicated situations analogous to those for a plane interface discussed by Danckwerts 
[4] and Crank [5]. Thus if the diffusion constant has a discontinuity at a concentration 
p. , We can use one expression of the form (9) in the region p < p, , another for p > p, . 
Conditions of continuity, the shape of the precipitate surface, a growth condition like 
(17) at the surface, and prescribed values of p at infinity, the discontinuity, and the 
surface then may be used to determine the parameters and integration constants in the 
solutions in the two regions. Alternatively, we can specify the normal velocity of two 
confocal surfaces in terms of the rate at which the diffusing material is transported 
across each surface, in analogy with Danckwerts’ “Class A” problems [4, 5]. If we do not 
relate the surface velocity to the rate of transport but instead determine our solutions 
by specifying densities at the surfaces, we have the generalization of ‘‘Class B” problems. 
Finally, we can treat the problem in which the growth of the new phase depends on the 
diffusion of two or more substances. This was solved by Frank for radial symmetry [2]. 

V. Anisotropic diffusivity. We shall now show that the preceding solutions are 
also appropriate when the diffusivity is anisotropic. Choosing coordinate axes such that 
the diffusivity tensor is in diagonal form, we have in place of (1). 


9 a2 Ps 9? 
f= D,—§+ D,—5 + D, - (25) 
t Ox oy Oz 
The transformation 
2’ = 2/Di”, (26) 
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then brings (25) to the form of (1) in terms of the primed coordinates, with D = 1. 
An ellipsoidal or hyperboloidal surface centered at the origin is described by the 
equation 
Ax’? + By’? + C2? + Day + Eyz + Fer — G = 0, (27) 
with suitable values of the coefficients. The principal axes of this surface are in general 
different from those of the diffusivity tensor. The boundary condition for the diffusion- 
limited growth of this surface is now , 
(p. — p,.(v-n) = n-DVp (28) 


in place of (12), v being as usual the velocity of a point on the surface, D the diffusivity 


tensor 
\D, O 0} 

>= | 0 D, O | (29) 
| 0 O D,| 


and n the vector 
n = (2Axr + Dy + Fz, 2By + Dz + Ez, 2Cz + Ey + Fz) (30) 
which is normal to (27) at the point r = (z, y, 2). The length of n is of no concern since 
(28) is linear in the components of n. 
If the surface preserves its shape during growth, v is proportional to the position 
vector Pr: 
v = Kf(or. (31) 
Substituting (29), (30), and (31) into (28) and introducing the change of variable (26), 
we obtain the boundary condition appropriate to shape-conserving growth expressed in 


the primed coordinates, 
(p. — p.)(v’-n’) = (n’-V"p). (32) 
Here n’ is the vector 
n’ = [2AD,2’ + D(D,D,)'’y’ + F(D,D;)'”2’, 
2BD.y’ + D(D,D,)'?x’ + E(D,D,)'’’z’, (33) 
2CDsz’ + E(D,D;)'?y’ + F(D,D;)'2'), 
which is normal to the surface into which (27) is transformed by (26) in the primed 
coordinates. Also 
v’ = Kf(ir’. (34) 


This boundary condition (32) is identical with (12) (with D = 1), however, if v’ is 
identified as the velocity of a point on the surface in the primed coordinates.* Since in 
this coordinate system the diffusivity is isotropic, the solution of (1) and (32) for constant 
initial density is given immediately by our earlier results. These determine K and f(t) 


*Note that v’ is not the vector obtained by substituting from (26) into the components of v in (31). 
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(= ¢-*) in (34), and these in turn determine v in (31). Thus on transforming back to the 
original coordinates we have the solution to the problem of the growing surface with 
anisotropic diffusion, and we see that it indeed corresponds to the preservation of the 
shape of the surface during growth.* Surfaces of constant concentration are not confocal 
with this boundary surface in the original coordinates, although they are in the co- 
ordinates for which the diffusivity is isotropic. 

VI. Discussion. Of the solutions of the time-dependent diffusion equation which we 
have exhibited, those pertaining to the growth of ellipsoidal particles are of particular 
interest in the theory of diffusion-limited precipitation from supersaturated solution. 
As the author has pointed out elsewhere in a discussion of various aspects of precipitation 
theory [7], the existence under appropriate conditions of solutions corresponding to the 
growth of ellipsoidal particles of constant shape indicates that in these situations the 
diffusion field is not responsible for any change that may actually occur in the particles’ 
shape. A change should be attributed to some other cause, therefore. This result contrasts 
with views that have been widely held in the absence of accurate solutions of the diffusion 
equation for asymmetric particles or an anisotropic diffusivity. Of course, one can not 
use the diffusion equation or simple boundary conditions to describe the initial stages of 
formation of an actual precipitate particle, and our solutions are exact for a uniform 
initial solute distribution only under the physically artificial restriction that the particles 
have zero dimensions at ¢ = 0. Our solutions are, however, useful in approximating real 
diffusion-limited growth when the particles have reached a size that permits the neglect 
of atomic fluctuations and other complications of the nucleation process and when 
transient effects associated with finite initial dimensions of the particles have sub- 
stantially vanished. 

The usefulness of these solutions for asymmetric surfaces is broadened by their 
applicability to problems for which the diffusivity is anisotropic. Thus the growth of a 
radially symmetric particle with anisotropic diffusion is mathematically equivalent to 
growth of an asymmetric particle with isotropic diffusion. Indeed, the general problem 
of an arbitrarily oriented ellipsoidal particle of zero initial dimensions with anisotropic 
diffusion can be solved exactly by means of this transformation, as we have shown. 
In all cases we find that the surface preserves its shape as it grows, the rate of transport 
of material across each point of the surface sufficing to give a local surface growth 
velocity appropriate to this result. 

The solutions for hyperboloidal surfaces have no immediate practical application 

*This result may appear at first sight to be incorrect in the limiting case in which one or more 
components of the diffusivity tensor are identically zero. Thus if one interprets the vanishing of D, and 
Dz in (29) as requiring that no transport occur except in the z-direction, one must conclude that a grow- 
ing particle extends its dimensions only parallel to the z-axis, so that its shape is not preserved. How- 
ever, a discontinuity in the solute density (or equivalently an infinite concentration gradient) then 
exists on the cylinder which is parallel to the z-axis and tangent to the particle. On the other hand, the 
solution of the text predicts in the limit (D, , D2 — 0) that, for example, a sphere of zero initial radius 
grows at a finite rate and remains spherical. Surfaces of constant concentration are spheroids of revolu- 
tion which are tangent to the sphere and have their major axes along the z-direction. An infinite con- 
centration gradient therefore exists at the circle of tangency, but in this limiting case the product of a 
vanishing component of diffusivity and an infinite gradient leads to a finite current density. The dif- 
ferent predictions for the limiting case therefore depend upon how one deals with an infinite concen- 
tration gradient. The result of the text is valid if we associate with an infinite gradient a current density 


obtained by approaching the limit from finite diffusivities. 
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so far as the author is aware. They do correspond to the diffusion-limited dissolution 
of a conical or wedge-shaped surface placed initially with its apex touching a plane on 
which the density is maintained constant by an external agency. However, our solutions 
require that the density vary appropriately in the region between the plane and the 
cone even at ¢ = 0, and this condition would be difficult to achieve in practice. 
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BOOK REVIEWS 


(Continued from p. 128) 


The gyroscope—theory and applications. By James B. Scarborough. Interscience Pub- 
lishers, New York, London, 1958. xii + 257 pp. $6.50. 


This is an introduction to the mechanics of the gyroscope, written on an elementary level and cover- 
ing a fair amount of theory, but with emphasis on applications. Part I (theory) contains a chapter on 
vector calculus, one on the fundamental principles of mechanics and two chapters on the theory of the 
gyroscope, treating the free gyro, forced precession and the heavy gyro. Part II (applications) is de- 
voted to all kinds of gyros on vehicles: gyro-compass, spherical pendulum, stabilizers of ships and of 
monorail cars. The last chapter is devoted to astronomical applications. 

One of the features of this book is its simplicity of reasoning. With very few exceptions everything 
follows from two simple sets of equations, (16.4) and (25.2). The treatment of the gyro-compass is 
particularly elegant, also the representation of the pseudo-regular precession by means of a Taylor 
expansion. The part of the book dealing with applications is very complete. Some numerical examples 
are helpful for the appreciation of the importance of various factors. Frequent references to modern 
developments are particularly valuable. 

In some instances, especially in the part of the book devoted to theory, the simplicity mentioned 
above is carried too far. The definitions of the concepts “gyroscope”’ (p. 37), “heavy gyro”’ (p. 65), 
‘spherical pendulum” and “gyroscopic pendulum” (p. 161) are rather vague. The pseudo-regular pre- 
cession of the heavy gyro is discussed under very restrictive initial conditions (p. 70). A gyro having 
two degrees of freedom should not be referred to as a gyro of one degree of freedom (p. 143). The sta- 
bility investigations on pp. 55 and 61 start from incorrect assumptions; the results are either incorrect 
or mere first approximations. Equations as (22.2) through (22.4) or (27.2) through (27.4) where one 
side is a scalar, the other one a vector, should be avoided. 

Though this list might be continued, the value of the book as a whole is beyond question. 

A suggestion concerning terminology: the term “‘precession”’ for a certain partial motion of a gyro 
stems from the verb ‘‘to precede’’ (lat.: praecedere); one should not derive a new verb ‘‘to vrecess’’ 
from this noun but rather say that a gyro precedes if it moves in a precession. 


H. ZIeGLER. 


A course in multivariate analysis. By M. G. Kendall. Hafner Publishing Co., New York, 
1957. 185 pp. $4.50. 


Until recently multivariate analysis has not been well covered by the text books in theoretical 
statistics. At almost the same time three monographs appear now in this field: the one by M. G. Ken- 
dall, another called Introduction to multivariate statistical analysis by T. W. Anderson and a third by 
S. N. Roy: Some aspects of multivariate analysis. The two latter ones are published by John Wiley 
& Sons. 

Dr. Kendall’s book is meant to be an introduction to the subject, perhaps especially for a reader 
interested more in applying these statistical methods than in their mathematical derivation. Indeed, 
one attractive feature of the book is the emphasis given to model construction and the interpretation 
of numerical results. The validity of the available techniques is critically examined, e.g. the centroid 
method in the first chapter, and alternative methods are examined and compared. A great number of 
detailed numerical examples illustrate the theoretical discussion and add a good deal to the value and 
usefulness of the book. The exposition is clear and the book reads very well. 

The book begins with a brief introduction describing what multivariate analysis is and to what fields 
it has been applied. The second chapter is devoted to component analysis and contains some elementary 
but basic results, which are used throughout the rest of the book. A brief discussion of factor analysis 
follows, and the fourth chapter is concerned with functional relationships with stochastic elements. 
This topic is often treated in an obscure way in the literature and often causes difficulties to the student 


(Continued on p. 164) 
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EIGENOSCILLATIONS OF AN ELASTIC CABLE* 


BY 
PENTTI LAASONEN** 
University of California, Los Angeles 


Summary. The eigenoscillations of a cable, supported at its end points, in a homo- 
geneous gravitational field is investigated under the assumption that the bending 
stiffness is negligible but the behavior with respect to tension stresses is perfectly elastic. 
The problem involves two coupled second order equations and one independent second 
order equation; it is shown to be definitely self-adjoint and an iterative method for its 
solution is suggested. 

In the particular case of a shallow cable, that is with negligible sag, the asymptotic 
eigenvalues are obtained. It turns out that the gravitational field has such a “‘stiffening”’ 
effect that the eigenvalues related to some oscillation modes may be substantially 
greater (the lowest of them eightfold) than those given in the classical theory on vibrating 


strings. 


Notations. 

A = cross sectional area 

L = total length 

s = arc length parameter of the cable when unstrained; 


II 


mass density 

pA = mass per unit length 

= Young’s modulus; 

zr, y, 2 rectilinear coordinates of a typical cable point; the gravity acting in the 
direction of negative y-axis, the end points of the cable having coordinates 
(+ 41, 0, 0); 

) displacement components from equilibrium position; 


ty F > 
II 


u,v, U 
= tension force acting on the cable; 
f = increment of the force due to the motion; 
a = length constant determined by Equation (3); 
e = aug/EA = apg/E equilibrium strain at the lowest point; 
a = L/2a = characteristic parameter; 
w = circular frequency of eigenoscillations; 
\* = w'a/g eigenvalue parameters; 
h = = w'L’/4a9 
p = s/a,q = 2s/L dimensionless variables; 


A, B, M,N, S, T matrices. 
Derivation of the equations. From Hooke’s law follows 
a? + yy? + 2? = (1 + F/EA)’. (1) 


*Received April 28, 1958. The preparation of this work was sponsored in part by the Office of Naval 


Research. 
**On leave of absence from the Finland Institute of Technology, Helsinki. 
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Hence the components of the tension force are 
Fz! Fy, Fz 


[+F/EA’ 1+F/BA’ 1+F/BA 


The influence of the flexural rigidity being neglected, Newton’s law gives 


(==;) = ” 
[+Frau,” 


Fi A ; ” 
(ra), = #00 + 9 " 
| 8 


Fz; ‘ 
a — a 
Fea) = ot. 
If time derivatives are taken to be identically zero, then Eqs. (1) and (2) give rise 
to the stationary equilibrium solution. Using boundary conditions and choosing the 
origin of s at the lowest point one obtains 


(a2, = a log [s/a + (1 + 8?/a’)'””] + 6s, 
| yo = (a? + 8”)'”? + 6°/2a + K, 


Zz = 0 
Fo = ug(a’ + 8°)’ 


The integration constants a and K in these expressions are determined by the conditions 





log [L/2a + (1 + L’/4a’)'”?] + Lpg/2E = 1/2a, (3) 
(1 + L’/4a’)'” + L’pg/8aE + K/a = 0. 
After introducing the small oscillation displacements by 
x(s, t) = 2(s8) + Us, bt) 
y(s, t) = ys) + v(s, t) 
2(s, t) = 2,(s) + w(s, t) 
[F(s, t) = Fo(s) + f(s, 


and substituting these expressions in (1) and (2) the linearization is performed: 








ut s/a_ _ ff 

(1+ s8°/a’)'’” EA | 
(. __ f+ FAQ + s'/a)uy _ y = ult | 
(i+ 3 a’)'* [1 4 e(] + 8", ao)"} : Mute y 
(sia + BA + 8/0 ay 
a+ 8°, a’) *[1 Sal + s*/a*)*”] MUit y 
oe Eo 
i+ +dit+ehy'y. * "° 


When f is eliminated, the dimensionless variable p = s/a introduced, and the variables 
separated by u(s, t) = u(s) - e’**, etc., then the following system of ordinary differential 
equations is obtained: 
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(2a tt 4 ns 28, 


(1 + p)[1 + 1 + p)'”] 
} 


(ee a toll + PY + ne 


(4) 
0, 


Il 


(1 +p )[l + (1 + p)'”] 


aA eee | iia “ 
( + py? + -) + A*u 0. (5) 


\ 


Two first equations apparently determine the modes and frequencies of the eigen- 
oscillations in the zy-plane, whereas the third equation applies to the perpendicular 
oscillations, completely independent of the first ones. In the following the main atten- 
tion will be given to the first oscillations. 

The self-adjointness of the problem. Since the eigenvalue problem, consisting of 
the two second order differential equations (4) and boundary conditions 

u(ta) = v(ta) = 0, (a = L/2a) 

may be regarded as arising from a regular variational problem, it is self-adjoint. The 
exact definition of a definitely self-adjoint problem, due to Bliss [1, 2], requires that if 


the problem is written in the form 
d 
|- u = Au + ABu, 
dp 


| Mu(a) + Nu(b) = 0, 


where u is an n-vector, A and B are (n X n) matrices and M and N constant (n X n) 
matrices with (M, N) being of rank n, then there exists a non-singular matrix T such 
that the following conditions are fulfilled: 


d ae 
[Ar+m+aT =o, (6) 
L TB + B’T = 0, 
MT '(a)M’ = NT'(b)N’, (7) 


(here the prime stands for the transposition of the matrix in question) and that the 
matrix 


S = TB 


is symmetric and definite or semi-definite. 
In the present case, solve the Eqs. (4) with respect to u’’ and v’’ and regard u, v, u’, 
and v’ as four elements of u. One finds then that the coefficient matrices are 


l 


A= +p)" + +p) 
0 0 0 0 
|0 0 0 0 | 
00 ptep2—pylt+p)'? p— el — 2p)\l+p)’ 2] 
0 0 -1— «1 — 2p)\(1 +p)” —p — 3ep(1 + p*)'” 
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| 0 0 0 0 

4 | 0 0 0 0 
(1 TP) | ag? — (1 +p)” p 0 0} 

p —-1—el+p)” 0 0 





| 
(M,N) = 
}0 000100 0; 
lo 000010 0) 
The solution for (6), (7) is now 
ey om 
(l+p)(1+e1+p)'”] 
[ 0 0 1+ai +)" +p 
| 0 0 +p p+el1+p)” 
}—-1 —el +p)” —p 0 0 
L —p —p —el+p)” 0 0 
and the matrix S is semi definite: 
(1 0 0 0 
lo 100 
Ss = 
00 0 0) 
00 0 0) 


The definitely self-adjoint character of the problem thus being proved, a number of 
properties associated with such problems may be applied for the present problem: for 
example, the existence of a countably infinite number of eigensolutions, the expansion 
theorem, etc. The orthogonality of two independent eigensolutions (u, , v;), (uz , v2) 
now takes the form 
[ (UjU2 + U,v2) dp = 0. 


a 


Furthermore, an iterative method for the determination of the eigensolutions is applicable 
and proved to be convergent [3]. The sequence of consecutive iterates is obtained directly 


from (4): 


oy 2 ) 
a | on eT er 
mee J ( i¢e’ *jj™ cp + (l+p)” J ‘. apf “P, 


ee ae a See 
— | (1 7am | m ip +( its c) f ved “D; 


the integration constants must be chosen to make the functions satisfy the boundary 


conditions. 








1959] EIGENOSCILLATIONS OF AN ELASTIC CABLE 151 


Finally, it may be seen directly from the form of the Eq. (5) that the problem for the 
determination of the lateral oscillations is definitely self-adjoint also and hence subject 
to similar properties and procedures as those mentioned above. 

Case of a shallow cable. The general solution of the problem above depends es- 
sentially on two dimensionless parameters only, the relative cable length a = L/2a 
and the minimum strain e. The latter is always smaller than the ratio between the yield 
stress and Young’s modulus and hence small compared with 1. In the case of a shallow 
cable the first parameter also is small compared with 1, since a is approximately equal 
to the smallest radius of curvature. For a closer study of this particular case we replace 


the variable p by the new variable 
q = 2s/L = p/a 
in order to have fixed boundary values gq = + 1. The Eqs. (4) and (5) thus obtain the 


form 


( u’ + aq’ + 1 +.°q’)*7u’ ) ie 0 
. + —5-1757 u = 0, 
e(1 4 aq ){l ++ e(1 + a g)'””] 7 


(8) 
agqu’ + aqv’ + 1 + Oar) 
| = > 9 > 9.472 dv = 0, 
(22 + a°q){1 + (1 + a°q')'”] i 
( — ona y + dw = 0 (9) 
(1 +aq) ™ +e ; 


where the prime now indicates differentiation with respect to g. The pertinent boundary 
conditions are 

u(t1) = vo(+1) = w(+)]) = 0. (10) 

After noticing that the last equation apparently has the asymptotic eigenvalues 

 ~ (nx/2)’, n=1,2,---, 
for a and ¢ tending to zero, we will leave this equation aside and study the problem 
(8), (10) for small values of a and e. By letting, separately, a — 0 one arrives at the 
equations 


( wu’? + edu 0, 


j 
ly + (1 + &)dv 


II 


0, 


which are not coupled and give the well-known eigenfrequencies of longitudinal and 
transversal oscillations. Since there is, for a small ¢, a wide gap between the smallest 
eigenvalues 7/4 and ~ x*/4 of both groups, it is of interest to find out what are actually 
the smallest eigenvalues, say when eA may be regarded to be still essentially smaller 
than 1. To this end observe, that the Eqs. (8) are equivalent with the following two 
equations, where C, and C, are arbitrary integration constants: 


u’ + aq’ + al | udq + aq | v in| + €', + €.aq = 0, (11) 
agqu’ ona yp’ . : Y 
: > 92 > / — ) anal , => 0. 1 
e+ (1 tag) ++ \E [ u dq [ 1 a | + Cag C; (12) 
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“ 


Then, considering first the solutions of group I, that is, those for which u is odd and v 
is even, one immediately finds that C, = 0. Assume that \ remains bounded for a and 
¢ tending to zero and, since both integral terms of Eq. (11) have ¢ as a factor, neglect 
these terms besides two first terms: 


u’ + aq’ + CC, = 0. (13) 
Further, substitute u from this equation into the Eq. (12) and neglect all terms which 
are not of the lowest order in a or e: 
ea 
vy +X | » dq — Cieq = 0. 
All even solutions of this equation which satisfy the condition v(1) = 0, are, up to a 
constant factor, of the form 
v = COS Kg — COSk, 
whereby C, and X are 
C, = —«’ cosk/a, A=K«k. 


The substitution into (13) gives odd solutions 


u = alsin xq/x — q cosxqg + eq cosk/a’]. 
From the condition u(1) = 0, « is determined as a solution of the characteristic equation 
tan k/k = 1 — x’e/a’. 
Hence, the roots of this equation squared give the eigenvalues. 
Quite similarly one finds that the asymptotic solutions of group II, whose wu is even 
and v odd, are 


RN =n, 
u = —(a/nr)(cos nx + cosnrg + nrq sin n7q), 
(vy = sin nrg. 


Hence in group I the solutions, the eigenvalues as well as the modes of eigenfunctions, 
depend on the ratio e/a’, whereas in group II the solutions are independent of this 
ratio. The dependence is illustrated in Fig 1, where the solid curves give some lowest 
eigenvalues from the group I, the broken lines eigenvalues from group IT. In particular, 
for a’/e = n’x’ there exist double eigenvalues n’z’. 

Finally, it may be indicated what the physical significance of the ratio ¢/a° is. Ex- 
pressing ¢ as the ratio ¢/E, where a is the equilibrium stress at the lowest point, and 
a by 

a = L/2a = Lpg/2E = Lpg/2c, 
we obtain 

e/a? = 40°/L’p’g’E = 4-(L,/L)’, 
where the characteristic length 


L, = (0/E)"-(a/pg) 











EIGENOSCILLATIONS OF AN ELASTIC CABLE 153 
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depends on the material in question and its equilibrium stress. For most materials this 
length, when computed for stresses reasonably close to their tensile proportionality 
limit, amounts to several hundreds of feet. For instance, for copper, « = 30,000 psi, 
L, = 330 feet, and for steel, c = 40,000 psi, ZL, = 430 feet. Therefore, the eigenfrequencies 
of a tightly stressed cable or wire do essentially differ from those corresponding to the 
well-known eigenvalues 4 = (nz/2)° only at spans of length comparable with this 
measure. 
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ASYMPTOTIC DEVELOPMENTS FOR A BOUNDARY VALUE PROBLEM 
CONTAINING A PARAMETER* 


BY 
RICHARD C. MacCAMY 
Carnegie Institute of Technology 


1. Introduction. It is the purpose of this paper to point out a curious type of singu- 
larity which can arise in the perturbation of the solutions of boundary value problems 
containing a parameter. We shall consider two such problems each for functions of z 
and y continuous in y < 0, (2, y) bounded, and harmonic in y < 0. The two functions 
satisfy the following sets of boundary condition, K being a positive constant and g(x) 
a given function: 

Problem I (A,) u,(z,0) = 0 for {|x|>1 


(By u,(x,0) + Ku(x,0) = g(x) for |x| <1, 


” 


Problem II (Ay) u,(x, 0) — Ku(zx, 0) = 0 for |x|>1 
(By) u,(x, 0) = g(x) for |x| <1. 


We admit the possibility of the function g(x) depending on K provided it be analytic 
in K for K sufficiently small. The two problems clearly have much the same character, 
the one deriving from the other essentially by interchanging the roles of the intervals 

z|<1land|z| > 1. In fact it can be shown that they are equivalent, that is solution 
of the one yields the solution of the other [1]. The two problems indicate strikingly 
the need for caution in the study of perturbations for despite their apparent similarity 
we shall find entirely different behavior of the solutions for small K. 

Physically (I) is a problem in heat conduction while (II) governs the diffraction of 
surface water waves by a rigid dock of finite width [2]. Differences in the problem become 
more apparent when we specify the behavior of the solutions of (I) and (II) for large 
x + y°. It is easily seen in fact that the following behaviors hold for (I) and (II) re- 
spectively, [2]. 


(C;) u(x, y) — ¢ log (2? + y*) = 0 + y)' as 2? +y’' — &, © constant. 


(Cy) ulz,y — Te*e'** = (2) as x—+o, y bounded 
r 


vy ik l 
u(x, y) — ReX*e'™* = of as x—-—o, y bounded, 
r 


for some constants T and R. 
If one sets K = 0 in (I) and (II) the problems become formally the same namely 


(A,) u,(z,0) = 0 for |x|>1 
*Received May 9, 1958. This work was supported by the Office of Naval Research, Contract Nonr- 


222(25). 
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(Bo) u,(z,0) = g(x) for |x| <1. 


This last problem is trivial having as solution the function*, 
Fy , 2 
U(x, y) = --| g(t) log [(a — )° + y'J dt. (1.1) 
v 1 


One may ask now whether the solutions of problems (I) and (II) tend to u,(z, y) as 
K — 0. A glance at conditions (C;) and C;,;) leads one to suspect that such is the case 
for the solution of (I) while for the solution of (II) difficulties might arise, a suspicion 
which proves correct. We shall see that the solution of (I) is in fact an analytic function, 
u(x, y; K) for complex K of sufficiently small absolute value, with u(a, y; 0) = wo(a, y). 
The solution of (II) is also analytic in K but for a region 0 < | K | < p. It has in fact 
the rather curious asymptotic expansion, f 


Ku(a, y; K) ~ Le >. Am (x, y)K"(K log K)" Ao = 0. (1.2) 


m=O nn0 
Thus the solution of (II) does not possess a limit as K — 0. The apparent physical 
paradox introduced in the finite dock problem can be removed by more careful attention 
to physical units but the odd development (1.2) for the solution of (II) remains. 
2. Developments for Problem (I). The development of the solution of Problem (I) 
for small K is elementary. To solve (I) we try the function, 


u(x, y; K) = -- | f(t; K) log [(x — #)? + y’] dt. (2.1) 
Ww J-1 


This function clearly is harmonic in y < 0 and satisfies (A,) and (C,). By a well known 
formula, it will satisfy (B,;) provided f(t; K) is a solution of the integral equation, 
ae & uh 3 
f(z; K) + | f(t; K) log | x — t! dt = g(x) lz| <1. (2.2) 
w Ju; 
The usual method of successive approximations establishes that (2.2) possesses a solution 
f(x; K) for suitably small K which is analytic in K. Further we have f(z, 0) = g(z). 
Substitution in (2.1) yields a solution u(x, y; K) of (1) analytic for K sufficiently small 
and with u(z, y;0) = u)(2, y). 
3. Developments for Problem (II). An integral representation for Problem (II) can 
be obtained in the usual manner by introducing the Green’s function G(z, y, t; K) defined 


by, 
. r\ 2 21 1 K Ky a Kn 2 2 
G(x, y, t; K) = —1/2 log [((x — O° + y'] + 9 ¢ | e log [((x — t)” + n°] dn (3.1) 
4. wien %e'* z-t ; 
We have, 
y vies y 
G = ( = «= ne © 
° KG (a-—t+y 


Thus setting, 


*For this we need Hélder continuity of g(z). 
tIt is to be noted that asymptotic developments of the form (1.2) have been observed before in 


Problem (II) but in an entirely different connection [3}. The meaning of the asymptotic development is 





2 


made clear in Sec. 3. 
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+1 
u(r, y;K) = —= f(t; KG, y, t; K) dt (3.2) 
we find, in the same manner as (2.2) was obtained from (2.1), 
u(x, 0; K) = f(x; K) + */ f(t; K)G(z, 0, t; K) dt on |x| <1 


while u(x, y; K) satisfies (A,;,) for any choice of f. It is further easy to show that the 
function defined by (3.2) is harmonic in y < 0, continuous in y < 0 and satisfies (C;,). 
Thus (3.2) yields a solution of Problem (II) if f(x; K) is a solution of the integral equation, 


g(x) = f(x; K) + & / f(t; K)G(z, 0, t; K) dt on |x| <1. (3.3) 


Note that the function f(x) is related to the solution u by, 
f(x; K) = u,(x, 0; K) — Ku(z, 0; K). (3.4) 


In order to express u(x, y; K) as a function of the parameter K we must so express 
G(x, y, t; K) Observe first that, 


ate” %e**\*-* _ ‘2 A,,(x = t, y)K", 
m=O 


0 


Fes , Kn og I: = 2 2 
K log ((x — 0)? + 9] dy (3.5) 


“y¥ 


-K™> at [a log ((@ — 0 + wl dn = D Bale, y, OK", 


m=} . Jy 


the series converging absolutely and uniformly for bounded z, y, t, and | K | The term, 
fim Ke | e-*" log [(x — 0)? + n°] dn (3.6) 
0 me 


is more complicated. Note first that, 


I = e* [ e” log E — f+ aa dr 


= —2e"" log K + ef e log [((K |x —t|)? + 7] dr, 


hence, 


lim (I + 2 leg K) = 2/ c*hetdo <S, (3.7) 
K-—0 0 
where y is Euler’s constant. 

To obtain a more detailed description of J we study it as a function of complex K. 
The integral (3.6) defining J converges, and yields an analytic function as long as 
0<|K|,0 <argK < x/2. I may be continued analytically to a larger sector by shifting 
the path of integration. In order to continue J(K) into arg K > 2/2 we shift the path 
into the sector arg 7 < 0 so that it keeps the points + 7 | x — ¢ | always to the right 
and runs to ~ along the ray arg » = — arg K, thus keeping Kn positive for large | 7 |. 
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In particular if 7° denotes the value of J after one counter-clockwise circuit of the origin 
we have, 


~ 


I* = Ke*” | e *" log [(a — 0)? + 7°] dn + Ke*” | e*" (log [(x — 0)? + »’)] dn, 


70 


where c denotes a closed path surrounding + 7 | x — ¢| in the negative sense and (log 
{(z — t)? + n’])~ denotes the value of the logarithm after 7 circles the points + 7| xz — ¢ | 
in the negative sense. Hence, 

r—t nilr—t 


-f[= oK eK, ni [ ad e *"dn + 7 | 


\ /0 70 


e*" dn — 2mi | eo in} 


= —4rie“” cos K | x — t}. (3.8) 


Equations (3.7) and (3.8) together imply that J(x, y, t) has the form, 


I = {-2e™"cosK |x —1t]|} log K + >. ets,» OK”, (3.9) 


m=O 


) ) = — 2y. Combining (3.5) and (3.9) we have finally, 


where by (3.7 


G(x, y, t; K) = {-e*” cosK | x —t!} log K + bo» Gils. 9, Dk” 
sisi (3.10) 


= —log K > H,(z, y, 0K" + >> G, (x, y, OK” 
the series converging for suitably small K. For reference we write down the first few terms: 
H, = -!1 H,=-y H, = —3|2—1%|? + iy’, 
Gy = (xi — y) — log [(x — ° + y’|' . 


For the case y = 0 all of the terms in the development may be obtained from formulae 
given in [4]. 
With the developments (3.10) at hand we are ready to discuss u(x, y; A). Suppose 
g(x) = Be 0 Im(x)K”. Then substituting (3.10) in (3.3) we have, 
f(z; K) - z. K"(K log K) | f(OH,,(x, 0, t) dt 
da a (3.11) 


a + 20 


+= > K” I 


Now we remark that the power products K”(K log)" can be ordered according to in- 


{()G,(x,0, t) dt = >> g,(2)K" on |a| <1. 


m=O 


creasing degree of vanishing as K — 0, i.e. 


. K™(K log K)’ - ‘ us 
lim >; B ——- = () if m’ +n’ >m+n or m+n’ =m+n and n>n’. 
xo A”(K log K) 
We proceed to develop the solution of (3.11) in the set of these power products. Passing 


to limit K = 0 in (3.11) we find, 


foo(x) = lim f(a; K) = g(x). 


K-—0 


Then (3.11) becomes, 
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f(x; K) — foo(x) — : K log K [ f(t; HH (x, 0, t) dt = of K log K) on |2| <1. 


Dividing by K log K and passing to limit K = 0, 
(®) = lim | — fol | a f wil in 
for(a im | K log K at 2 foo(t) dt = es go(t) dt. 


Reentering in (3.11), 
SF 

f(a; K) — foolx) — fo(a)K log K + [ f(t; G(x, 0, t) 
Ww J-1 


= g(x)K + of K) on |2z{ <1. 


Dividing by K and passing to limit K = 0, 
+1 


fio(z) = lim [f(z; K) — fool) — for(z)K log K|/K = g,(z) -} | foo(H\Go(x, 0, t) dt. 
J-1 


K-0 
We have indicated here the first three steps of a process which can be continued 
indefinitely. The process leads ultimately to, 


f(a; K) ~ > D foa(z)K"(K log K)". (3.12) 


m=0 n=0 


That the successive f,,,(2) can be determined by recursion one sees as follows. We say 
f(x; K) has an estimate of degree (m, n) if, 


f(x; K) = P(K, K log K) + o(K"(K log K)"), 
where P is a polynomial (with coefficients depending on x) of degree (m, n). Since, 
KG = —K log K + o(K log kK) 


we see that the product of KG with a polynomial of degree (m, n) is a polynomial of 


degree (m, n + 1) plus terms o(K”"(K log K)"*'). Suppose then that we have shown 
f(x; K) to have an estimate of degree (m, n), that is have computed foo , for , «++ 5 fmn iN 


(3.12). Substituting this estimate in the integral in (3.11) we obtain a polynomial of 
degree (m, n + 1) plus terms of o([K”"(K log K)"*"], with coefficients determined by the 
known quantities, foo, --- , fmm » Right hand side of (3.11) has estimates of degree (m, 0) 
for all m, hence (3.11) yields for f(x; K) an estimate of degree (m,n + 1). 

Substituting the series (3.12) and (3.10) in the integral representation (3.2) it is 


readily seen that u(x, y; K) has the form, 


Ku(r,y; K) ~ > > u,(2, yK"(K log K)". (3.13) 
mee acd 

4. Convergence of the asymptotic developments.* The development (3.13) is asymp- 
totic and the question arises as to whether the series actually converges to yield the 
function u(x, y; K). For the special case g(x) = — e’**, which governs the diffraction of 
waves by a rigid dock, an alternative procedure exists for obtaining the development 
(3.13). This second technique answers the convergence question while at the same time 

pointing up the relations which exist between the coefficients of the development. 


*The idea of this section was suggested to the author by Professor Hans Lewy. 
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For simplicity we consider only the even part of the solution that is the solution 
of Problem (II) with g(x) = — cos Kx. The odd part could be treated in the same 
manner. The solution has the representation (3.2) with f(z; K) a solution of the equation, 

. i ih a ‘ F , 
f(x; K) + f(t; K)G(a, 0, t; K) dt = —cos Kz on zx <1. (4.1) 
WT J-1 
G(x, 0, t; K) is defined as a function of complex K by (3.10), that is, 


G(x, 0, t; K) = —{cosK | x — t|} log K + >> G,(z, 0; t)K”. (4.2) 
m=O 
The usual method of successive approximations applied to (4.1) shows that the solution, 
as a uniform limit of functions analytic in K, is analytic for | K | ¥ 0 sufficiently small. 
Let K*, K~ denote the image points of K on a logarithmic Riemann surface of K after 
positive and negative circuits respectively of the origin. By (4.2), 


G(x, 0, t; K*) — G(z, 0, t; K) = —2xi cos K | x — t | 
hence if we form (4.1) for K” we find, 


f(z, K") + E | f(t, K°)G(a, 0, t; K) dt (4.3) 


+1 
— 2K [ f(t, K*) cos K | x — t| dt = —cos Kz. 


J=1 


Subtracting (4.1) from (4.3) and noting that f is even with respect to ¢ 


fiz, K*) — fiz, K) +>] f(t, K*) — f(t, K)G(a, 0, t; K) dt 
~ a+] 
= 2K | f(t, K°) cosK |x —t| dt = 21K cos Kx | f(t, K’) cos Kt dt 
e 1 « 1 
hence the difference, f(z, K”) — f(x, K) satisfies, except for a multiplicative constant 
the same integral equation (4.1), as f(x, K) itself. Therefore, 
+1 


f(x, K*) — f(x, K) = 2iKf(z, K) | f(t, K*) eos Kidt 


1 


or, 
atl 


f(x, K) — f(x, K-) = 2Kf(z; A) | f(t, K) cos Kt dt. (4.4) 


; l 
The function f(z, K) thus satisfies a non-linear difference equation in K. If we let, 
l 


MK) = | f(t, K) cos Kt dt 
we find by multiplying (4.4) by cos Az and integrating, 
MK) — NK) = 2KNAK)MK). (4.5) 
This difference equation for \(K) is still non-linear but we can form from it a linear 


equation. As in the derivation of (3.12) we have from (4.1), 
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f(t, K) = -—1+ 011) as K-00. 


Hence, 


MK) = —2+0(1) as K-0 (4.6) 


so that \(K) [and also \(K_)] is different from 0 for K sufficiently small. Thus for small 
K, (4.5) can be divided by \(K) \(K~) to yield, 
1 1 ; 
NK) ~ NK) > 2K. 

It follows that, 


ee log K + S(K) or XK) = — “= , (4.7) 
* S(K) - = log K 





where S(K) is single-valued, and being bounded by (4.6), is a power series, i.e., 








S(K) = 1 (i roe » SK°), (4.8) 
n=1 
Substituting (4.7) in (4.4) we find, 
f(z, K) -]1 + 2k f(z, K-) = 0. (4.9) 
S(K) — * log K 
Tv 
Now observe that, 
7(K~) = S(K) - K log K + 2iK = 7(K)} 1+ =k " 
i S(K) — = log K 


Also 7(0) ¥ 0, hence for sufficiently small K (4.9) can be written 
7(K)f(z, K) — 7(K")f(z, K-) = 0 


i.e., 7(K)f(z, K) is a single-valued function, p(z; K), which is continuous at K = 0, 


and thus, 





f(x, K) = ek p(t, K) = > p,(2)K". (4.10) 
S(K) — = log K -_ 


From (4.10) two interesting facts may be observed. First we can write, since S(0) = 3, 
1 
a 
f(x, K) 2p(x; K) 1 am Z(K) ’ 
where 


ZK) = 2™ log K — > S,K". 


n=1 


Since | Z(K) | < 1 for K sufficiently small, say | K | < Ky, we have, 
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f(a, K) = 2p(x; K) >> (Z(K))” (4.11) 
which it is readily seen leads to a convergent double power series in K and (K log K). 
Second we see by (4.10) and (4.11) that f(z, K) is completely determined by the two 
analytic functions p(x, K) and S(K). That is to say in the double series development 


f(a, K) = >> > a,,.(z)K"(K log K)” (4.12) 
m=0 n=0 

for f(z, K) there are at most two simple infinities of the coefficients a,,, which are sufficient 
to determine the others. In reality the analytic functions p(z; K) and S(K) are related 
so that only a single infinity of the coefficients need be determined. One sees this by 
using (4.11) to express the S, and p,(x) as functionals of f(z: K), i.e., of the a,,;(x). If 
the operation in (4.11) is carried out and the series (4.12) is substituted on the left, we 
obtain a series of relations, the first few of which we write down here: 


2po(x) = Aoo(x), 


4po(x) 
TT 


— Qo,(x), 
2p,(x) — 2po(x)S, = aro(2), 


8 ; 
“3 PolX) = Aor , 
T 





8S, pol 
4pi(x) _ 8Sipolz) _ aed 
T T 
We see from these relations first that the p,(x) are determined by the a;o(x) and the 
S, and then that the a;;(x) for 7 > 1 are determined by the a,,(x). The p,(x) and {S,} 
and hence f(z, K) are completely specified if one knows the values of the a;.(x) ¢ = 0, 
1,2 +++. 

The intermediate step of computing the a,;(x) may be eliminated as follows. Sub- 
stituting (4.10) and (3.10) in (4.1) yields, 


+1 


(c,K) +2 p(t, K{—K log K (cosK | 2 —t}) + 2G, (x, 0, OK} dt 
P T 


=] m=O 
i eee K e 
= —COos Kal SK) —-— log K | on |z| <1. 
T 
We can equate here terms which involve log K and those which do not, obtaining 
4 f—) ? lo} 


-:/ p(t, K) coosK|x—t|dt= “ cos Kx on |2x| <1 (4.13) 
Tr J-1 


ry ] ps ry - f \ u7m+ Y ry r 
p(x, K) + =| p(t, K) >> G(x, 0, )K"*' dt = —S(K) cosKr. (4.14) 
wT J-1 


m=0 


Since p(z, K) is an even function of x (4.13) may be rewritten, 


+1 
| p(t, K) cos Kt dt = —1. (4.15) 
-1 
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Set p*(x, K) = p(x, K)/S(K). Recalling that S(O) ~ 0 p*(z, K) is an analytic function 
of K for K small, i.e., 


p*(z, K) = 3 p*(x)K". 


n=0 


Substituting in (4.14) and equating coefficients of K’, r = 0, 1, 2, --+ , yields, p§ = — 1, 
r-1 +1 
pea) +2 Ff plOG,-s-s(x, 0, 0) at 
TW n=-0 J-1 


(—1)"*" ie 


= ——— 2x” ifr is even, r>0 (4.16) 


2r! 
= 0 if r is odd. 
Equation (4.16) determines the p,(x) recursively. For example: 
p(x) = 
n(x) = if p*()Go(x, 0, t) dt = 2 Gi ~ 


—} [1 + 2) log (1 + 2) + (1 — 2) log (1 — 2) — 22}. 


Once the p*(x) are found Eq. (4.15) can be used to determined S(K). Substituting 
p(x, K) = S(K)p*(z; K) in (4.15) and equating coefficients of K" yields, 


1 ‘Se (—7)” f* . 2 
= > i on" dt * 
2 p> 2m! I PP-an()t™ dt 


(4.17) 





+1 
a > = 8," | Pinot" dt=0, i=1,2,3,- 
: -1 


Equation (4.17) determines the S, recursively since the coefficient of S, is, 


+1 
[ pxo at = -2 #0. 
1 


For example: 
, oo: a 
S, =i['> ()a=i-2% 5 (log 2 1). 
Once S(K) and p*(z; K) are found, p(z; K) is obtained as their product. 
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BOOK REVIEWS 


(Continued from p. 146) 


of statistics. Canonical analysis is discussed and illustrated by examples in Chapter 5, while Chapters 
6 and 8 deal with distribution problems and tests of multivariate hypotheses. Among other topics 
treated in this book are tests of homogeneity, discriminant functions and some history of multivariate 
analysis. An extensive bibliography and a set of exercises end the book. 

Ur GRENANDER 


The analysis of multiple time-series. By M. H. Quenouille. Hafner Publishing Co., New 
York, 1957. 105 pp. $4.75. 


The last decades have witnessed rapid progress in statistical time series analysis, as far as a single 
series is concerned. This was made possible by the results in probability theory describing the proba- 
bility structure of time series of various types. For multiple (vector valued) time series the progress 
has been slower, one of the reasons perhaps being that we know less about the basic properties of vector 
valued stochastic processes. It is typical that the linear prediction problem solved in the early forties 
for scalar processes, is still not quite settled in the vector case, although recent work of Wiener and 
Masani may change this. 

In this book the author collects available methods and develops new ones for analyzing multiple 
time series. His starting point is the assumption that the underlying process is of the finite parameter 
type: a moving average, an autoregressive scheme or a combination of both. This restriction may be 
reasonable in the present preliminary stage of development of the theory. 

After a brief discussion of these schemes and of their covariance properties and identification prob- 
lems, the author studies the numerical behavior of five artificially generated series. The next chapter, 
Practical complications, contains a discussion of how various factors, may influence the anaylsis. Such 
situations occur e.g. when 1) the errors are serially correlated 2) the sample size is too small, so that 
large sample approximations are inadequate or 3) there are unknown trends or seasonal variation. 
After two chapters on estimation and testing the author investigates an econometric example in detail 
trying to apply the methods described in the book. 

The emphasis of the book is on the formal development of statistical techniques, and one would 
have appreciated a more complete critical examination of the validity of the methods used. However 
it is clearly too early to ask for definitive and convincing general methods of analysis of multiple time 
series, and the need for further work is obvious. At present we will have to use the tools of analysis 
available to us, even if they leave a good deal to be desired. When applied with some caution they 
should enable us to extract some information from the data. 

This book is the first in the new series ‘‘Griffin’s statistical monographs and courses’’ in which has 
also appeared the book on multivariate analysis reviewed elsewhere in this journal. 

Ur GRENANDER 


Applied mathematics for engineers and physicists. By Louis A. Pipes. McGraw-Hill Book 
Co., New York, Toronto, London, 1958. xi + 723 pp. $8.75. 


This second edition of a book first published in 1946 remains essentially the same as the first edition 
as regards subject matter, style and purpose. New material has been added to most chapters making 
the book somewhat larger than its predecessor. Particularly there is a greater emphasis on methods of 
numerical analysis some of which have been developed during the last decade as a result of the stimulus 
furnished by the high speed computers. 

The sections on theory and applications of matrices and of Laplace transforms have been revised 
and contain much material not usually found in a text of this type. The book continues, however, to give 


recipes for solving specific problems with considerable disdain for sound mathematics. 
G. F. NEWELL 


(Continued on p. 172) 
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MATRIX DIFFERENTIAL SYSTEMS WITH A PARAMETER IN THE BOUNDARY 
CONDITIONS AND RELATED VIBRATION PROBLEMS* 


BY 
JULIAN ADEM 
Instituto de Geofisica, Universidad Nacional de M éxico 
and 


Instituto Nacional de la Investigacién Cientifica 
México, D.F., (MEXICO) 


In a series of papers by Moshinsky [1, 3] and Adem and Moshinsky [2] the descrip- 
tion of some physical processes has been given by a column matrix or vector instead: 
of a single function, achieving a simplification in the mathematical solution of such 
processes and making it possible to solve a variety of problems. In the examples con- 
sidered in [2, 3], for example, the above formulation leads to a matrix eigenvalue problem 
which is a generalization of the Sturm-Liouville problem and that can be treated along 
similar lines. 

In this paper, by describing a vibratory system in the aforementioned way, we are 
led to a more general eigenvalue problem that contains a parameter both in the differ- 
ential equations and in the boundary conditions. This problem can be considered as a 
generalization of the corresponding scalar case, examples of which have been considered 
recently by Bauer [4] and Morgan [5]. 

In the first part of this paper, we give examples of vibratory systems that give rise 
to this new type of eigenvalue problem. 

In the second part of this paper, we formulate and discuss the general matrix differ- 
ential system and obtain some properties of the eigenvaJues and eigen column matrices. 

I. Vibrating string with a concentrated mass at a point. Consider a string with a 
concentrated mass at a point which separates two regions of different density and length 
of which the string is composed. The string is fixed at the ends. We wish to find the 
transverse displacements of the string if it is displaced initially into a position and 
released from rest at this position with no external forces acting. 

Let the length of the string be L, + ZL, , and M the concentrated mass at a distance 
L, from the left end point of the string. We shall designate I and II as the regions to 
the left and right of the concentrated mass respectively, and Y, and Y, as the corre- 
sponding displacements. Region I has a density p, , and a length LZ, , while II has a 
density p, and a length L, . The space-independent variable shall be called x, in I and 
t_ in II; x, increases from left to right, while x. increases from right to left, so that at 
the end points we have z, = 0, x, = 0 and at the point with the concentrated mass, 
z= L,,%=1,. 

Considering that in both regions the equation of a vibrating string applies and 
considering boundary and initial conditions we obtain the following differential equations 


and boundary conditions: 


*Received Tune 11, 1958. Numbers in square brackets refer to the bibliography at the end of the 


paper. 
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(2 Y 
OX, 


Y, 


where Ci = T/p,; (i = 


initial displacements prescribed at regions I and II respectively. 


Ys) .: 
‘) T ( O2%e/s,=Ls ™ 
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ay, 


‘ 


OX; 7 
Y,(0, t) = 
Y(L, ’ t) = 


ay, 


1 
7 at” 





(x, , 0) = F,(x,); 


(2%) ~ (2%) - 
wn Veta; 








| 

fy 

- (=e) | 
hf 


Y.(x, , 0) = F,(22), 
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(1) 


) 
| 


| 


(2) 


(3) 
(4) 


1, 2); T is the tension on the string; F,(x,) and F,(x.) are the 


Applying the change of independent variables 


2, = L,2} (¢ = 1, 2), 


we have that now the range of x{ is the same as that of xj and we can therefore put 


% = 2; 
Problem (1), (2), (3), 


where 


and 


Let 


= Z. 


1 ( 21) 1 ( ous) in (=H) | 
Fi \ C2) ens + | T hae Fo 


(4), becomes: 


] 9 - 
Sl “Se —- 0 a - 1 
ax! Ci ar * a<l, 
y,(0, t) = y.(0, 2) = 0 
y.(1, 2) = y(1, 


yi(x, 0) = f(z); 0) = f,(z), 


(Oy; /Ob) pao = (OY2/dt), «0 , 


y2(x, 


y(x, t) = Y,(z,; , d) 
f(x) = F,(x,). 


yz, t) = yx) cos wt 


then problem (1’), (2’) becomes: 


s 


( 


\L, dx 


1 d’y, 


w” , 

t = = 9 
L dx’ + cc Yi 0 (2 i, “); 
y,(0) = y.(0) = 0, 
yi(1) = y2(1), 


a dy; 





L, da 


l 2s _ Fey) = 0 
z=1 


(1’) 


) 


| 
fy 
| 
/ 


(3’) 
(4’) 


(5) 


(6) 
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The solution of Eq. (5) is 
y(ax) = A, cos @L,;z/C,;) + B, sin @L,;x/C,), 
where A, , B; and w are going to be determined from the boundary conditions. 


From conditions (6) it follows that A; = 0. 
From conditions (7) we get the frequency equation: 


(1/C,) sin @L,/C;) cos wL,/C,) + (1/C;) cos @L,/C,) sin @WL,/C,) (8) 
— (Mw/T) sin @L,/C,) sin @L,/C,) = 0 


whose roots are real and form an infinite denumerable set. Let the positive roots of 
Eq. (8) be w; , w., *** , @,, °** . They are the eigenvalues of the problem (5), (6), (7). 
To each eigenvalue there corresponds a solution: 


Y;;(x) - sin (w;L,2/C,) 
i , (9) 
Yo(x) = oe Sata(Ce sin (w,;L,2/C,) 


sin (w;L./C2) J 

where we have used the notation y,;(x), y2;(x) instead of y,(x), y2(x) to distinguish 
that it is a solution corresponding to the w; eigenvalue. y,;(x) and y2;(x) can be con- 
sidered as the components of a column matrix (vector) which will be denoted by y;(z) 


y (zx) _ wn 
Yo;(2). 


and called eigen column matrix: 


Proceeding as usual, we assume as solution of our problem an infinite series of the 
type: 


y(x, t) = ° Byy,(z) cos wt. 


i=1 

This series is a solution of (1’), (2’), (4’) provided it is convergent and admits two 

successive term by term differentiations with respect to ¢ and zx. Assuming that this 

condition is satisfied, there remains, then, to determine the constants B; , so that (3’) 
is satisfied: 

f(2) = y oa = 8B, eae = > By((2). (10) 

fo(x) ~? Yo;(2) = 
Assuming the validity of expansion (10) the coefficients B; can easily be determined 
and for it we shall establish first the orthogonality relation of the eigen column matrices 


corresponding to the problem (5), (6), (7). 
Since y,(x) and y,(z) are solutions of (5) corresponding to «, and w, respectively, 


the following relation is satisfied: 


=~ d’y; d’y; 2. L, 
pn  ® (vi a = a vs») - > Ce (we = ws) YirY ip e 








t=1 


Therefore 








i= dys |' yf yk: 
) | Vir iy ys | = @ - 09 [ Do G2 Yinis de. 


a~ 7T, | Y* ae dx |, 
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From boundary conditions (6) and (7) it follows that 


2 1 
> [ Bya@y.@ az + yy.) = 0 (11) 
¢=1 70 Ci 7 


or with matrix notation: 


[ yi@)Hy,(2) de + yi(1)Jy,(1) = 


where 


od 0 . y= [07 0 
0 L,/Cc?J 0 O 


and y{(x) denotes the transposed of y,(x). 

This is an orthogonality condition for the eigen column matrices of problem (5), (6), (7). 
Assuming that the series (10) can be integrated term by term after being multiplied 

by y/(x)H, we can now, using (11), determine the coefficients B; in (10), which are 


given by 


| ((L,/C?)f,(x) y(t) + (L2/C2)f2(x)yo;(2)] dx + (M/T)f,(1)y,,() 











B; = : al (12) 
| ((Li/Ciyii(a) + (Le/CDy2(a)] de + (M/T)y?,(2) 
or more briefly with matrix notation 
1 
[ yi@Ht@) ax + #OI¥.) 
B; = e é 
[ vi@ry.a) ae + yidy) 
0 
Substituting (9) into (12) we obtain 
] n? @ids |} de SM 
2 E +s Ce l= sin” (w;L,/C,) + ri 
— Lyf 4 cay cin WHA ae 4. Ea Bin (wsdn/C,) wy La 
= 7a | f(x) sn —— C, dx + C? sin (w, mle f.(x) sin C, dx (13) 


+ 7 f(1) sin 2 


and the formal solution to the problem (1), (2), (3), (4), is 


orn 
| ¥i(z, , t)| 


YZ , i) 


© ai ; Ti 
> B, om ores/Cy) COS w,t, 
sin (w,L,/C2) 


sin (w,x,/C2) 
4 


where w, are the positive roots of Eq. (8) and the B,’s are given by (13). 
II. Solution of other vibrating systems. The above method can be applied to the 
case of a string with n concentrated masses M, , M, , --- , M, which separate n + 1 
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regions of different density p, , p2, --* , p, ; and different length L, , L,, --- , Z,.In 
this case we describe the system by n + 1 functions Y, , Y, , --- , which by a change of 
variable is always possible to have in the same domain, for example, 0 < zx < 1,t > O. 
We call the n + 1 regions from left to right, I, II, III, --- and, for example, in the 
regions labeled by an odd number let the space variable increase from left to right, 
while in those labeled by an even number let it increase from right to left. We are led 
to a system of n + 1 differential equations of the type (1’) with boundary conditions 
similar to conditions (2’), (3), (4’). In fact, we are led to a matrix differential system 
which is a particular case of the one discussed in the next section 

Since the solution of the mathematical problem corresponding to a vibrating string 
also solves analogous problems of torsional vibrations of shafts of circular section and 
longitudinal vibrations of bars, the above method yields exact solutions (without neglect- 
ing the mass of the bars, as is usual in Theory of Vibrations) of such systems when 
they have several concentrated masses, which separate bars of different elasticity, 
section and length 

III. General matrix differential system. The problem (5), (6), (7) as well as those that 
arise in the examples of Sec. II, are particular cases of the system formed by the general 
self-adjoint matrix differential equation [6] 


uo) = 9 (PE) + OF + (58 + Rp = Wy ay 
and boundary conditions of the type 
(A + AB)y(0) + C(dy/dz),-5 = ¥e (5) 
(D + AE)y(1) + F(dy/dz),., = 0) 


where P, R and W are n X n symmetric matrices, Q n X n antisymmetric matrix; P, 
R, W, Q are real functions of z; the unknown y is a column matrix of n components; 
\ is a scalar parameter; A, B, C, D, E, F, are n X n constant real matrices. We assume 
that P, R, W, Q, dP/dz, dQ/dz are continuous and that the determinant of P is different. 
from zero in0 < xz < 1. 

From (15) we obtain 


| ly‘L(y,) — L'(y.y,] dx 
. (16) 


= [y/P(dy,/dx) — (dy!/dx)Py; + y/Qy;]o = (A; — Ad) [ y:Wy, dz, 


where y,; and y, are two different solutions (eigen column matrices) corresponding to 
the eigenvalues A, and , respectively; and y; and L’(y,) denote the transposed of y, 
and L(y,) respectively (i.e., they are row matrices of n elements). 

The problem that we are going to deal with is such that by substituting the boundary 
conditions (15) into (16) we obtain a relation of the type: 


ly: ,yJ] = [ yi(xz) W(x)y,(x) dx + y{(0)My;(0) + yi(1)Ny,(1) = 0, (17) 


where M and N are constant n X n matrices. Relation (17) can be considered as an 
orthogonality relation that the eigen column matrices of the problem (14), (15) satisfy. 
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We are more familiar with orthogonality relations in which the last two terms of the 
left member of Eq. (17) are zero. Such is the case in the self-adjoint problem of the 
type (14), (15) but in which the parameter \ does not appear in the boundary con- 
ditions (i.e., B = EH = 0). This problem has been extensively studied by G. D. Birkhoff 
and R. E. Langer [7]; and also by Moshinsky and the author in connection with self- 
adjointness [6], and examples of it have been given in [2, 3]. As far as the present writer 
is aware the problem (14), (15) has only been studied in the scalar case (i.e., when 
nm = 1) and some papers by Bauer and Morgan have appeared recently dealing with 
it [4, 5]. 

Now that we have formulated the mathematical problem (14), (15) the next step 
would be to develop a theory dealing with properties of eigenvalues and eigen column 
matrices and with the possibility of expanding an arbitrary column matrix in terms 
of the eigen column matrices, as is needed for the applications, as, for example, in the 
problems mentioned in the above sections. 

Proceeding as W F. Bauer does in the scalar case [4], we shall give the following 
definitions: 

D1. A column matrix u(z) shall be called V-column matrix if it is real, not identically 
zero, of class C’, and satisfies the boundary conditions (15). 

D.2. The inner product corresponding to the eigenvalue problem (14), (15) is defined 


as 


al 
(h, g] = | h’(x)W(x)g(x) dx + h’(0)Mg(0) + h’(1)Ng(1), (18) 
70 
where h(x) and g(x) are column matrices whose components are bounded integrable 
functions 
D.3. The eigenvalue problem (14), (15) is normal if for every V-column matrix 
u(x) we have 
fu, uj] > 0. 


The following results can be established: 

T.1. A sufficient condition for the problem (14), (15) to be normal is that W be a 
non-zero diagonal matrix whose components are non-negative in the fundamental 
interval and N and M diagonal matrices whose components are non-negative. 

T.2. If W, N and M satisfy the same condition as in T.1., the eigenvalue problem 
(14), (15) has only real eigenvalues. The proof of T.2., is analogous to that employed 
in the scalar case [4]. 

T.3. If we define as Rayleigh’s quotient R(y) for any V-column matrix y(z) 


; (L(y), y] 

*O) = ‘Twy, 9] ’ 
where L(y) stands for the left member of (14), from (14) it follows that R(y;) = i, , 
where y; is an eigen column matrix corresponding to the eigenvalue ), . 

T.4. If W, N and M satisfy the same conditions as in T.1., and [L(y), y] > 0 for 
every V-column matrix then the eigenvalue problem (14), (15) has only positive eigen- 
values. This result follows from T.3. 

T.5. Corresponding to each real eigenvalue there is one or at the most, a finite number 
h, of linearly independent real eigen column matrices y{*’ (x), a =1, 2, +--+ , hy. (a; < 2n). 
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T.6. The eigen column matrices corresponding to different eigenvalues are orthogonal 
with respect to the inner product defined by (18). The eigen column matrices corre- 
sponding to the same eigenvalue can be normalized (we assume that the norm of all 
eigen column matrices is different from zero) and made mutually orthogonal, by the 
usual procedure [8] so that our set of eigen column matrices satisfies 


[yi yi?) = Sap bi; , (19) 
where 
a=}, 225° B= 1,2,--- ,h;; 1,j7=1,2,- 


and 5, is equal to 1 for a = 8 and equal to zero for a ¥ 8. 
T.7. If we assume the validity of an expansion of the type 


o hn 


f(z) = 2) An” yn”), (20) 
n=1 a=l 
where f(x) is a column matrix whose components are bounded integrable functions of 
x [that satisfy certain conditions that ensure the possibility of expansion (20)], and 
if the series (20) can be integrated term by term after being multiplied by y,/(r)W(z), 
it is possible, using (19), to obtain the coefficients A” , which are given by 


A® = [f(x), yo’}. 


The conditions under which the expansion (20) is possible as well as the existence 
of an infinite denumerable set of eigenvalues, and some properties of eigenvalues and 
eigen column matrices of the problem (14), (15) will be left unsettled. 


REFERENCES 


. Marcos Moshinsky, Phys. Rev. 81, 347 (1951) 

. Julidén Adem and Marcos Moshinsky, Quart. Appl. Math. 9, 424 (1951) 

Marcos Moshinsky, Bol. Soc. Mat. Mexicana 4, No. 1, 4, p. 1 (1947) 

W. F. Bauer, Quart. Appl. Math. 11, 274, (1953); see also the references mentioned in this paper 

G. W. Morgan, Quart. Appl. Math. 11, 157 (1953) 

. Juli4n Adem and Marcos Moshinsky, Bol. Soc. Mat. Mexicana 7, No. 3, 4, p. 1 (1950) 

. G. D. Birkhoff and R, E, Langer, Proc. Am. Acad. Arts Sci. 58, 51 (1923) 

R. Courant and D. Hilbert, Methods of Mathematical Physics, 1st. English ed., Interscience Pub- 
lishers, Inc., 1953, p. 50 


ore ON 


OsI1 & 








172 BOOK REVIEWS [Vol. XVII, No. 2 


BOOK REVIEWS 


(Continued from p. 164) 


Logical design of digital computers. By Montgomery Phister, Jr. John Wiley & Sons, 
Inc., New York, and Chapman & Hall, Ltd., London, 1958. xvi + 408 pp. $10.50. 


This book contains a development of a methodology of functional design of digital circuits and 
gives examples of application of the techniques to design of digital computer components. Discussion 
is limited to synchronous circuit design. The book begins with a brief discussion of the system design 
problem and the logical designers role, and of circuit components and the binary number system. Next, 
the fundamentals of Boolean Algebra are presented, followed by a description of various methods of 
simplifying Boolean functions. The general methods of application of Boolean Algebra in sequential 
circuit design are developed and utilized in a detailed discussion of the design of the major digital com- 
puter components; memory unit, input-output unit, arithmetic unit and control unit. Examples of 
complete computer design are included. 

The book provides a compilation of the techniques available in the field of logical design and, be- 
cause of both content and organization, should be useful to the logical design student and to the prac- 
ticing designer. The references to literature in the field and the supply of exercises for the student are 
comprehensive. 

The book serves as an excellent compilation of techniques, particularly in its exposition of the meth- 
ods of reduction of Boolean functions and the general methods of application. However, motivation 
for the specification of characteristics of the digital computer components is somewhat weak. This fact 
is implicitly recognized by the author in the discussion of the logical designer’s role as one who is re- 
sponsible for developing the implementation of a set of requirements provided by the system designer: 
development of characteristics is not fundamentally in the logical design province. The result of the 
lack of motivation and discussion of the prescription of system parameters limits the usefulness of the 
book to those who have some experience in the computer field. The book would not be useful, for ex- 
ample, as a text in a first course on digital computer organization and use, nor in a course on computer 
system design, but should be reserved for use in an advanced course where the objective is solely logical 


design. 
DEAN GILLETTE 


Introduction to difference equations. By Samuel Goldberg. John Wiley & Sons, Inc., New 
York, and Chapman & Hall, Ltd., London, 1958. xii + 260 pp. $6.75. 


This volume was motivated by a monograph written at the invitation of the Social Science Re- 
search Council Committee on the Mathematical Training of Social Scientists. It is intended as an intro- 
duction to the theory and application of recurrence relations of the form tuna yi: = f(u,) for students and 
research workers in the social sciences who have only a rudimentary knowledge of mathematics. 

It is a pleasure to say that the author has been eminently successful! in his task. The volume is not 
only easily intelligible, but full of important applications which illustrate the significance of this field. 
The book should certainly be useful as a basis of a course following after one patterned along the lines 
of the recent book by Kemeny, Snell and Tompson. 

Furthermore, the book would be equally useful as an introductory junior or senior course for physi- 
cal science majors. In these days of narrow specialization, it is essential that research workers in one 
domain get a glimpse of how scientists in other fields construct their mathematical models. 


RICHARD BELLMAN 


(Continued on p. 184) 
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ALGEBRAIC TOPOLOGICAL METHODS FOR CONTACT NETWORK 
ANALYSIS AND SYNTHESIS* 


BY 
CHARLES SALTZER 


Case Institute of Technology 
and 
Electronics Laboratory, General Electric Company 


1. Introduction. The use of Boolean algebra in the analysis and design of contact 
networks was introduced by C. Shannon [1], and for series-parallel circuits the problem 
was completely solved. However, for circuits of a more general nature, e.g., bridge 
networks, the design problem was not solved and the rules given for analysis depended 
on special geometric constructions carried out on the circuit diagram. The treatment 
of general circuits was approached from a topological point of view by the author [2] 
and S. Okada [3]. The method of Okada is analogous to the “loop” or “mesh” method 
of ordinary circuit analysis while the method of this paper is analogous to the ‘‘nodal’”’ 
method [4] of ordinary circuit theory. It will be shown that contact network theory can 
be based on topological considerations and topological methods for analysis and synthesis 
will be given. 

A contact network between two terminals is closed for a given state of the contacts 
if and only if there is a path between the two terminals such that each branch of the 
path is itself closed. With no loss of generality each branch may be regarded as con- 
taining a single contact and a branch is closed or open according to whether the contact 
is closed or open. The analysis problem consists of finding a formula for the state of the 
network in terms of the states of the contacts. The synthesis problem is: given such a 
formula, find a contact network which realizes it, subject usually to some sort of re- 
striction such as minimality of the number of contacts. 

The topological preliminaries will be given in Sec. 2. The analysis problem will be 
treated in Sec. 3 and the conversion of the topological formulas to formulas of Boolean 
algebra will be carried out in Sec. 4. The design problem will be treated in Sec. 5 and an 
example of the synthesis method is given in the last section. 

2. Topological preliminaries. A graph is defined as a collection of segments or 


branches denoted by 2, , %2 , -:* , 2 , and nodes which include the endpoints of the 
branches. The nodes will be denoted by y, , Ye , -** , Ya - We consider linear forms in 
these variables choosing the coefficients from a modulo two arithmetic**, i.e., 

0+0=0 

Oo+1=1+0=1 

1+1=0 

(1) 

0-0=0 

0-1l1l=1-0=0 

l1-l=1] 


*Received May 15, 1958. 
**The linear forms in the x; z2**+2» and the linear forms in y; , y2**-y, each form a group. 
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The linear forms will be called 0-chains [5, 6], C, if the terms of the sum are nodes, 
and 1-chains, ¢c, if the terms are branches. The structure of a graph is known—and the 
graph can be drawn if the relation of the nodes to the branches is given. This information 
can be described by the boundary operator 0 which is defined for segments by 


Ot, = Ye +Y (2) 
if y, and y, are the endpoints of x, and for chains by linearity, i.e., 
Oy, + y,) = OY, + OY, . (3) 


These relations may also be expressed in matrix form by the incidence matrix 
(4) A = (A,,) where 
\l if y, is an endpoint of z, , 
4. - | (4) 


(0 if not 


\ 
The rows of this matrix correspond to branches and the columns to nodes. There are 
exactly two non-zero elements in each row in the matrix e.g., for Eq. (2), the kth row 
will have ones in the rth and sth columns only. If the boundary of a chain is zero the 


chain will be said to be a cycle. In the bridge circuit of Fig. 1 the 1-chain ¢ = x, + 2, + a 
is a cycle since do = 02, + 02%, + 0%, = Yr + Yrotytystysty = 0. 








Y3 


Fia. 1. 


The boundary relations may be represented by 
Of, = 2. ArsYe (Y= 1,2, --+ , dB). 
s=] 


The incidence matrix for this bridge is: 


.. 2 a 0] 
}O 1 1 0) 
{1=/0 1 0 1| (5) 
}1 010 
10 01 1 


An alternative description is also given by the coboundary operator 6 defined for nodes by 
6(node) = sum of the branches having this node as an endpoint, (6) 


and for chains by linearity, i.e., 
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d(y, + y.) = dy, + dy, . (6’) 
Thus, in Fig. 1 
Yeo=1+tw+%, 
and 
iy +ty)=U1tuatat+w+4%3=%+%,4+97;. 
The incidence matrix A’, for this operator is defined by: 


(1 if z, has y, as an endpoint 
Ai, = 4 (7) 
| 0 if not 


A comparison with (4) shows that A’ is the transpose of A. 
In terms of the matrix elements the coboundary operator may be represented by: 


b 

by, = Z. , ee (rf = 1,2, --- ,n). (8) 
= 

A chain is a cocycle if its coboundary is zero. For example, C = y, + y2 + y¥3 + % 
is a cocycle. A pair of nodes, y, and y, are said to be connected if there is a l-chain ¢ 
such that do = y, + y, . In Fig. 1, y, and y2 are connected since dz, = y, + y, ; and 
y, and y, are connected since 0(x, + x3) = y, + yy. One such chain consists of the sum 
of the segments of any path connecting the nodes in question. Also, two nodes are said 
to be homologous if they are connected and when a graph is said to be connected it 
means that every pair of its nodes is connected. The lemma on which our considerations 


are based is: 


Lemma. A zero chain is a cocycle if and only tf, for every node in the zero chain, every 
node connected (homologous) to the given node is also present [6]. 


For completeness the proof is included. With no loss of generality we may consider 
only a connected graph. If the 0-chain consists of the sum of all the nodes then since 
every branch is a term of each of the coboundaries of exactly two nodes it is a con- 
sequence of the modulo-two addition that this chain is a cocycle. Conversely, if a node 
y, is in a cocycle then every node to which it is connected by a branch must also be in 
the cocycle; for if not, then the branch of which such a node is an endpoint would be 
present in the coboundary only once. In general, if y, is connected to y, then there is a 
chain which is a path z, + x + --- +2, such thatd(z,+2%,+-:-+"%) =YyY%t+Y,, 
and every node which is an endpoint of a branch of this chain must be in the cocycle 
since, if this were not the case, some branch would appear only once in the coboundary. 
Hence, for any given node, the cocycle must contain every node connected to the given 
node. 

3. Topological analysis of contact networks. An immediate consequence of the 


above lemma is 
Theorem I. Two nodes of a graph are connected if and only if there is no cocycle which 


contains one of these nodes but not the other. 


Thus, if we have a graph with n nodes and say we wish to determine whether y, is 
connected with y, , we form all the 0-chains which contain y, but not y,—there are 2"~” 
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such chains—and find the coboundaries of these chains. If no coboundary vanishes, i.e., 
is a cocycle, then the two nodes are connected. We may indicate this symbolically by 
writing the formal product of these coboundaries with the convention that the product 
is to be zero if and only if at least one factor is zero, and our result may be stated as: 
y, and y, are connected if and only if the product of the given coboundaries is not zero. 

An explicit formula can be written by listing the integers from zero to (2""? — 1) in 
binary form and prefixing zeros when necessary to make the total number of digits equal 
to (n — 2). We define the numbers 


a,, = (s — 1)** digit (counting from the left) of r in binary form 
[s = 2,8, --+ ,n—1;r=0,1, --+ ,(2”’ — DJ. (9) 
a, = 1 [f= 0,1, --+ ,(2”* — 1)]. 
All the 0-chains containing y, but not y, are 


n-1 
C,= Daw. , (rf = 0,1, --- , (2%? — 1). " (10) 
s=1 


Hence the function { which we shall call the state-function of the graph with respect 
to nodes y, and y, is 


(2"-%=1) 


f(a, »%p,°%** > h%) = I] 5C, P (11) 


r=0 


and is zero if and only if y, and y, are not connected. 
By formula (8) 


n—-1 n—-1l b 
0, = Sa.ty, = a,( Aut), 
s=1 


s=1 k=l 


Therefore, 


b n—1 
a. = >» (S arene , 
k s=1 


=1 


or 
iC, = >» AriLe » 
where 


n-l 
a= > a,,4,, xf =0,1,-° , (2% — 1k =1,2, ---, dj. (12) 


Thus in terms of the incidence matrix the state function f can be written 


(2%~8—1) b 
f(x, » Xo ? oe? , Xp) = I] (> rats). (13) 
k=1 


r=0 


For the circuit in Fig. 1 considered as a circuit between y, and y, , the numbers in 
binary form from zero to 2”? — 1 = 2? — 1 = 8 are 00, 01, 10, and 11. Hence the 


numbers a,, are given by 
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1 0 0| 
fa as 101 
110 
lyao4 


and the zero chains including y, but not y, are 
Co=% 5 
Cr=yHt+y, 
C,="H7t+h, 
C,=HA+t+ytys.- 

The coboundaries are 
a, = 1 +%, 
iC, = 4% + 224+ 23, 
C. = %+23,t+%; 
6C; = %3 + 75, 

and the state-function f is 

f(a: , Xo, e+ 5 Hs) = (4, + X4)(X, + Lo + Xs)(Xo + Xs + 1s)(Z3 + 25). 


Precisely the same state-function except for the order of the factors is obtained by 
excluding y, and including y, in all the zero chains. This can be seen from the fact that 
the desired 0-chains can be obtained by adding y, + y, to the given chains, since by 
the modulo-two addition this replaces y, by y, . The coboundaries of these chains are, 
therefore, obtained from the coboundaries of the given chains by adding 6(y, + y,). 
Since the sum of all the nodes is a cocycle we have 6(y, + y,) = rb ew y.). Hence we 
ean get the desired coboundaries from the given ones by adding b 5 y to the given 
0-chains instead of y, + y, . The numbers a,, which we get by this computation are 
obtained from the binary representation of r by adding one to each of the (n — 2) right- 
most digits of r modulo-two. The resulting sums must be distinct since if they were not, 
the original numbers would not have been distinct. Since there are 2"? — 1 of these, 
this list must be a permutation of the original list and hence the new (A,,) matrix is a 
row-permutation of the original and we get the same state-function. 

We define a subgraph of a given graph as the graph obtained by deleting any branches 
(but not the nodes of these branches) from the given graph; and we have 


Theorem II. The state-function f of a subgraph with respect to a given pair of nodes is 
obtained from the state-function f of the given graph with respect to the same nodes by setting 
all the branches not in the subgraph equal to zero (or deleting them) in the state-function 
of the given graph. 


The incidence matrix of the subgraph is obtained from the incidence matrix of the 
graph for the boundary operator by replacing the ones in the rows corresponding to the 
deleted branches by zeros, and the coboundaries of the nodes with respect to the sub- 
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graph are obtained from the coboundaries of the nodes with respect to the graph by 
setting them equal to zero (or deleting them) in the latter. Hence the coboundaries of 
the 0-chains C, with respect to the subgraph are obtained by setting the deleted branches 
of the graph equal to zero in the coboundaries of the 0-chains of the original graph. 
Consequently by the definition of the state-function, the theorem is proved. 

In the graph of Fig. 1, for example, if we delete branches x, and 2, the incidence 


matrix a, of the subgraph becomes 


the coboundaries become 


6C, = %2 + 343, 
6C; = 13 + 2s ; 
and the state-function becomes 
f, (Ze » Ze » Zs) = JO, Ze , Xs , 0, 25) = O, 
and we see that nodes y, and y, are not connected. 
This interpretation of the formula can be stated as 


Theorem III. The connection between two nodes is zero when a given set of branches is 
opened (deleted) if and only if this set contains all the branches in at least one of the 0-chain 
coboundaries. Conversely, the given nodes are connected if and only tf at least one branch in 
each of the 0-chain coboundaries is not opened provided the state-function ts not identically 


zero. 

This theorem shows how to write a state-function from the conditions of operation. 

We define equivalent state-functions as functions which are zero for precisely the 
same sets of branch’states and Theorem III has as a consequence a minimization method, 
namely, 

Theorem IV. If the coboundary of any 0-chain, t.e., factor of f, ts contained in another 
factor then the state function obtained from the given state-function by deleting the first factor 
ts equivalent to the given state-function. 

It is clear that the results of combinatorial switching theory can be recovered from 
this point of view. 

As the last two examples we take branches in series and parallel. 


For the series network of Fig. 2(a) the above method yields f(x, , x2.) = x,x, and for 


the parallel network of Fig. 2(b), f(z, , #2) = x, + 2,. As in the method of Boolean 
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x 
y x Yo I Y3 
l nm j Wi Cc» Y2 


te 





(a) (5) 


Fia, 2. 


algebra, if we consider the states of the branches as being determined by individual 
state-functions, we can write the state functions of series-parallel networks directly by 
inspection. 

Thus, in addition to the explicit formulas for the states of general networks the 
preceding theorems show that the theory of contact networks can be founded on an 
algebraic topological basis. It is not suggested that this approach be a substitute for the 
methods of Shannon but rather a supplement, and in the next section we consider the 
conversion of these formulas into the formulas of Boolean algebra. 

4. The reduction of state formulas of Boolean algebra. Following Shannon (1.c.) 
we introduce the Boolean elements 0* and 1* where 


0* @ 0* = 0*, 
0* @ 1* = 1* @O* = 1*@ 1* = 1*, 
0* - OF = OF - 1* = 1* - OF = OF, 


oF a ar oe RM, 


The admittance between two nodes is 1* if the circuit is closed and 0* if the circuit is 
open. Also if 8 is a Boolean variable which assumes the value 1* when a given relay is 
operated and 0* otherwise, then the states of the normally open contacts of the relay are 
denoted by 8 and the states of the normally closed contacts are denoted by the comple- 
ment of 8, i.e., 8. If we denote the state of the branch 2, of the graph of the contact 
network by a, then the state-function is transformed into a Boolean function and we 
note that the value of this Boolean function is 0*, if and only if at least one factor is 0*. 
A factor is 0* if and only if each of its terms is 0*. By comparison with the results of 
the preceding section we see that the value of this Boolean function is 0* if and only if the 
circuit between the given nodes is open. Consequently, if the circuit is closed each factor 
has the value 1*, and the value of the Boolean function is 1*. This completes the proof of 


Theorem V. The Boolean admittance function for a pair of terminals of a network is 
obtained from the state-function with respect to those terminals by replacing the branches 
by their admittances. 


The calculations are simplified if the columns of the incidence matrix A are used as 
column vectors representing the coboundaries. 

It may also be noted that the admittances between different pairs of nodes of a 
network are related in a simple way so that if we have one state-function we can get the 
others immediately. This can be seen from the following considerations. Assume that 
the factors of the state-function with respect to terminals y, and y, are known. Assume 
that these have been obtained by including y, in all the 0-chains since we get the same 
results by using either y, or y, . To compute the state-function between y, and y, we 
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need the 0-chains which include y, but not y, , i.e., we must replace y, in the first set of 
0-chains by y, . Since y, occurs in C, multiplied by a,, , the replacement can be effected 
by adding a,, (yz + y,) to C, . Hence the factors of the desired state-function are obtained 
by adding a,, (6y, + éy,) to 6C, . A similar argument holds for replacing node y, . 


< 


Theorem VI. The state-function for the nodes y; and y,; can be obtained from the state- 
function for the nodes y; and y, by adding a,, (dy, + 6y;) to the factors 56C, of the given 
state-function. 

The matrix of these state-functions when the branches are replaced by their ad- 
mittances is the admittance matrix of Lunts [7] if the diagonal elements are taken as 1*. 

5. The synthesis problem for a two-terminal network. A method for finding general 
networks, i.e., not necessarily series-parallel networks, which realizes the graph of a 
contact network for a given two-terminal Boolean admittance will be given in this 
section. We proceed by reversing the analysis in Secs. 3 and 4. The first step is to write 
the given function as a product of sums of admittances. The second step is to assign 
branches to each admittance. At this stage there may be a considerable number of 
possibilities. One extreme is to assign a different branch to each appearance of a contact 
admittance or its complement; the other extreme is to assign the same branch to each 
appearance of a variable. The third step consists of writing a state-function by replacing 
the admittance variables by the corresponding branch variables. The fourth step is to 
select an integer n such that 6 + 1 is greater than n and 2” is not less than the number 
of factors of the state-function. 

In the fifth step if 2""* is greater than the number of factors of the state-function, 
then factors subject only to the restriction that the terms of such a factor contain at 
least the set of nodes in one of the given factors are added so the total number of factors 
is 2”-*. The sixth step is to assign a number r [r = 0, 1, 2, --- (2""* —1)] to each of the 
factors. In the seventh step the equations for the elements of the incidence matrix are 
written using Eq. (12). The last step is to solve these equations. If the solution has the 
property that there are no more than two non-zero elements in any row of the solution 
displayed in matrix form, i.e., is a partial incidence matrix, a column is added so that an 
element of this column is 0 if there are two ones or no ones in this row, and | if there is 
only one 1 in this row. If the partial incidence matrix does not satisfy the above con- 
dition then the computations can be repeated for possibly different assignments of r, 
b, n, and correspondences of the admittances and branches. An example will be given 
below. That the above problem has at least one solution is obvious since a series-parallel 
network can always be drawn. 

To clarify the method we retrieve the bridge network of Sec. 1 where f(a, , 
Qa, °** , @s) = (a, + ay) (a; + a2 + a5) (a2 + a3 + a4) (a3 + a;). We choose n = 4, 
b = 5 and replace a, by x, . The numbers r = 0, 1, 2, 3 are assigned so that 


6Cyo = 11 + % 
6C, = 2%, + a+ 72; 
6C, = x7. +4,;+ 2%, 


5C; = 2%; +25. 
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Thus the matrices A and \ are given by 


[1 0 0] 1001 0 
| 
| 
a es | 0 an h=(A,,) = l l 0 0 1) 
co | lo 1 1 1 0| 
| 
11] loo101 


The solution of this system is 
An = Au» 
Ans = Yow + Axx ; 
Ais = Aor + Anz > 


which in matrix form is the partial incidence matrix 


[1 1 0] 
01 1) 
010 \ 
101 
00 1) 


We get the incidence matrix by adding a column which brings the number of ones in 


each row to two, 


j1 1 0 0] 
0 11 01 
lo 101 
1010 
lo 0 1 1) 


which is the original matrix. The contact network can be drawn immediately. 
The equations for the general case may be restated as 
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n—-1 
> a,.dy, = 6C, ir =0,1,--- ,(2"? — Dj]. (14) 
s=1 
From the definition of a,, we see that if the system has a solution it must be 
bY; = 6U, , 
642 = 5C'o + 5Co0-: , 


We call these chains a reduced set of factors. Since the ordering of the factors in the 
state-function is at our disposal we may determine solutions by choosing one of the 


factors as 6C, and adding this factor to each of the others. If from these (2"* — 1) 
chains we can select (x — 2) chains such that no branch occurs more than twice as a 
term in the set of chains consisting of these (n — 2) chains together with the chain 


designated as 6C,, then we have a partial incidence matrix. The consistency of the 
set of equations (14) must of course be verified by substitution of the trial solution. 

Thus for a given assignment of branches, nodes, and order of factors the synthesis 
problem for a given state-function can be solved if and only if the system of equations 
(14) is consistent and the set of coboundaries (15) satisfies the condition that no branch 
appears more than twice. 

The consistency equations (14) may be regarded from another point of view. If we 
add 6C, (mod 2) to each of the equations of (14), the resulting system is generated from 


the (n — 2) sets dy. , dy, , «++ 6Y,-, by taking all possible sums, i.e., the set (6C, + 6C,) 
[k = 0, 1, 2, --- (2"* — 1)] is isomorphic to a Boolean algebra with (n — 2) atoms. 


The determination of whether or not this is the case can be shown to be independent of 
the choice of 6C, , and since the order of the factors of the state-function can be assigned 
arbitrarily, the consistency can also be shown to be independent of the ordering of the 
factors. Thus the test of consistency requires only that an arbitrarily selected factor be 
added to all the factors and any (n — 1) of these sums (modulo 2) be chosen as atoms. 
If each of these sums is a linear combination of the atoms then the set of equations (14) 
is consistent, otherwise the system is inconsistent and the assignment of branches to the 
admittances must be made in another way. 

If the number of factors of the state-function is less than 2"-* then the additional 
factors have the admittance one when the circuit is closed. Thus these factors must 
arise from branches whose admittances are complementary or must contain as a subset 
the branches which are present in another factor. This can be verified by generating the 
complete reduced set (6C, + 6C,) (r = 1, 2, --- , 2"? — 1) adding 6C, and verifying 
the above condition. It follows that for a given assignment of branches and choice of n 
the verification of the consistency equations (14) can be made independently of the 
choice of 6C, , i.e., holds for all choices, and when the consistency equations are satisfied 
the problem is reduced to determining whether an assignment of the factors of the 
state-function can be made so that each branch occurs no more than twice in these 
factors. This may not be the case even for a consistent set. 

It may be remarked that the complexity of this method is due in part to the fact 
that all circuits realizing the given function can be obtained this way. 
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BOOK REVIEWS 


(Continued from p. 172) 


gas dynamics. Edited by H. W. Emmons (Volume III of high speed 


Fundamentals of 
and jet propulsion). Princeton University Press, 1958. xiii + 749 pp. 


aerodynamics 
$20.00. 

a pleasure to record the appearance of another volume of this monumental series. The title is 
not too informative, and the first duty of the reviewer is to give a short list of contents and 
The equations of gas dynamics (Tsien); One-dimensional treatment of steady gas dynamics 


It is 
perhaps 


authors: 
(Crocco); One-dimensional treatment of nonsteady gas dynamics (Kantrowitz); The basic theory of 


gasdynamic discontinuities (Hayes); Shock wave interactions (Polachek and Seeger); Condensation 
phenomena in high speed flows (Stever); Introduction to combustion (v. Kaérmaén); Gas dynamics of 
flame fronts (Emmons); Gas dynamics of detonations (Taylor and Tankin); Flow of rarified gases 
(Schaaf and Chambré). 

The author of a “handbook”’ article always has to face the choice between the restricted aim of pro- 
viding an introduction to his subject and guide to the literature, and the ambitious aim of giving a defi- 
nitive and exhaustive account. Tsien has chosen the more modest alternative and contributed a lucid 
description of the equations underlying the ‘classical’ part of Gas Dynamics, with indications of where 
to look for the ‘“‘modern”’ departures. The stress is less on the derivation of the Navier-Stokes equations 
than on the discussion of the derived systems of general equations, such as the circulation and entropy 
theorems, the stream-function approach and the variational formulation for potential flow. 

Crocco has not only chosen the other alternative, but indeed interpreted the title of his article in 
its full, literal generality. There are no steady compressible flow processes in ducted machinery so 
complicated that useful one-dimensional models could not be inferred from a physical understanding 
of the essential mechanisms—and here the models are displayed in a sweep never attempted before, 
over the whole gamut from the ideal nozzle to the ducted rocket motor. The article is a book in its own 
right, with a number of original papers thrown in for good measure, and one surmises that only illness 
prevented Professor Crocco from adding a handbook of chemical engineering! The approach is broad 
physical discussion of complex viscous, thermodynamical and chemical interactions 


and versatile 
and the number of topics caught in the Crocco net defies enumer- 


alternate with engineering charts 
ation here. No doubt, this article by itself will ensure a heavy demand for the volume in the aeronautical, 
mechanical and chemical engineering industries. 

A particular trait of the ‘‘one-dimensional’’ theory has always been that its conceptual basis is 
vague, and even this problem is attacked in the article. The starting point is the recognition that uni- 
formity of physical quantities, or even physical mechanisms, over any cross-section of the duct is not 
a necessary premise. Large parts of the theory deal with the interaction between quite different streams 
flowing side by side. The outline of a general mathematical theory of the application of the conserva- 
tion principles to mean values representing such flows is sketched and elaborated in the direction of a 
classification of processes. 

Kantrowitz’ choice, on the other hand, is a very personal introduction to the propagation of plane- 
wave pulses and the stability of transonic duct flow. The treatment is physically illuminating, even if 
it lacks the unambiguous lucidity of Tsien. The article also contains a brief review of the procedures of 
the numerical method of characteristics and a sketch of a few results relating to strong shock propagation, 

After the leisurely style of Kantrowitz, Hayes’ telegram on the theory of shock, detonation and 
deflagration fronts comes indeed as a shock. The first part is devoted to a thorough discussion of the 
system of equations representing the conservation principles and the additional relations which make 
the system determinate—one might call it the general theory of the Hugoniot curve. In the second 
part, the ordinary differential equation governing the structure of the front on the basis of the Navier- 
Stokes equations for one-dimensional steady motion is discussed, both for the shock and the front with 
one-parameter exothermic reaction. For weak shocks, the description is extended to the case of slow 
variation in time. As with all telegrams, the onus of deciphering is on the reader—the reviewer passes 
on with the feeling that his appreciation of this article is rather incomplete. 

The discussion of the systems of equations governing the reflection and refraction of shocks follows 


(Continued on p. 219) 
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FOLDING OF A LAYERED VISCOELASTIC MEDIUM DERIVED FROM 
AN EXACT STABILITY THEORY OF A CONTINUUM UNDER 
INITIAL STRESS* 


By 
M. A. BIOT 
Shell Development Company, New York 


Abstract. A layer of viscoelastic material embedded in an infinite medium also 
viscoelastic tends to fold through instability when the system is compressed in a direction 
parallel with the layer. This phenomenon which was treated previously by using an 
approximate plate theory, [l, 2], is analyzed here by using an exact theory for the 
deformation of a continuum under prestress [3, 4, 5, 6]. The effect of the compressive 
prestress in the embedding medium is taken into account, and it is found that although it 
is not always negligible, it tends to be very small under conditions where the instability 
of the layer is strong. The same conclusion holds for the error involved in the use of a 
plate theory for the layer instead of the exact equations for a prestressed continuum. 
In the course of the analysis we have also treated the problem of a semi-infinite visco- 
elastic half space subject to a uniform internal compression parallel with the boundary 
and a surface load normal to this boundary. The compressive load produces an increase 
of the surface deflection under the normal load. This effect appears through an ampli- 
fication factor which is evaluated numerically for the particular example of an elastic 
body. It is shown that under certain conditions the free surface of the compressed semi- 
infinite medium may become unstable and will tend to wrinkle. This is suggested as a 
probable explanation for the wrinkles which appear on the surface of a body subjected 
to a plastic compression. 

1. Introduction. We consider a layer of viscoelastic material surrounded by an un- 
bounded medium also viscoelastic, the whole system being subject to compressive stresses 
parallel with the direction of the layer. Such stresses may be set up for instance if the 
whole system is subject to a uniform compressive strain parallel with the layer. 

As another example we have the case where a compressive stress in the layer alone 
will be produced if a swelling of the layer occurs while longitudinal expansion is prevented. 
The question which we are concerned with here is the stability of the layer under such 
a system of prestress and the prediction of the deformation of this sytem. We have already 
analyzed this problem by introducing certain approximations [1, 2]. In particular we 
have assumed that the layer behaves like a plate, i.e., that it obeys equations of flexure 
of the type used in strength of materials theory and generalized to apply to viscoelastic 
media. Another assumption introduced in this previous work is that the prestress in the 
surrounding medium has a negligible effect. 

Our purpose is to check the accuracy of the previous simplified theories by solving 
the problem from the exact equations for small deformations of a continuum under 
prestress. General equations for incremental deformations of a body under prestress 
were derived by this writer some years ago in a series of publications [3, 4, 5, 6]. While 
the theory was developed in the particular context of elastic continua, the equations 


*Received June 26, 1958. 
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are applicable to any type of incremental deformations of elastic or anelastic bodies. 
This is due to the fact that in the development of the theory the physics is separated 
from the geometry of the problem, and the incremental stress strain relations are left 
arbitrary. For the reader’s convenience, we have rederived the equations in Sec. 2 for 
the particular case of plane strain. In Sec. 3 these equations are first applied to the 
evaluation of the surface deflection under a normal load for a semi-infinite viscoelastic 
medium, which is at the same time subject to a compression in a direction parallel with 
the boundary. 

The relations between the incremental stresses and strains are assumed to be the 
same as those derived by the present writer from thermodynamics [7]. The nature of 
the approximation which may be involved in this assumption is discussed in a previous 
publication [2]. It is found that the lateral compression produces an amplification of 
the surface deflection and its magnitude is evaluated. The possibility of an instability 
of the free surface under a lateral compression is analyzed in Sec. 4. It is pointed out that 
although the magnitude of this instability is in general quite small, the effect may be 
large for incremental deformations in the plastic range where the so-called tangent 
modulus is the significant parameter. The instability should manifest itself in the form 
of a wrinkling of the surface, a phenomenon which is actually observed on the surface 
of a compressed solid in the plastic range. 

In Sec. 5 a complete analysis is given for the folding of a layer of uniform thickness 
embedded in a medium of infinite extent. We have assumed that there is no friction at 
the interface of the layer and the medium. This affords a direct comparison with the 
results obtained in [1] and [2] by an approximate theory and with the same assumption 
of perfect interfacial slip.* Results of the previous section are used in order to evaluate 
the influence of the compressive stress in the surrounding medium. It is found that in 
general this compression may be neglected, as it does not affect substantially the wave 
length of the folding. The case where the layer and the surrounding medium are both 
incompressible viscous fluids is investigated in detail. 

The rate of growth of the amplitude of folding as a function of the wave length has 
been derived. The dominant wave length—i.e., that which exhibits the maximum rate 
of amplitude growth—is evaluated as a function of the viscosity ratio of the two media. 
This dominant wave length is independent of the load. Results are compared with the 
results of the approximate theory of [1]. It is found that this approximate theory agrees 
very well with the present results when the !ayer viscosity is larger than about seventy 
times the viscosity of the surrounding medium. This is also the condition for the magni- 
tude of the instability to become significant so that the folding amplitude increases at a 
rate fast enough to be observable in practice. 

2. The theory of deformation of a prestressed solid. We shall first recall briefly 
the theory of incremental stress and deformation of a solid under initial stress. We 
assume the incremental deformation to be small so that the theory is linear. Such a 
theory was developed by the writer [3, 4, 5, 6] for the case of an elastic body. However, 
it is independent of the physical nature of the body and is therefore applicable to visco- 
elastic media. In order to simplify the presentation, we shall limit ourselves here to a 
two dimensional deformation. 

We are dealing with the problem of a continuum undergoing small deformation 


*An evaluation of the effect of interfacial adherence has been given in [2]. 


to 
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from an initial state which is not one of zero stress. The initial stress field is defined 
by components S,; , Sis = Sa; , So , referred to x, y axes. This being a two dimensional 
case, we assume that the other components of the initial stress field vanish or do not 
appear in the equations. The condition of equilibrium of this stress field is 





Sa. dSre ; 
oon + 908 + pX = 0, 

al (2.1) 
AS2, | 9S22 a 


with a mass density p and a body force X, Y per unit mass. The problem is to formulate 
the laws of deformation of this medium where small incremental stresses and strains 
are superimposed on the initial state. 

In treating this problem we have departed from the traditional viewpoint of the 
mathematician and defined the incremental tensors in a different way. Let us first 
consider the strain. A point P of initial coordinates x, y in the prestressed state is 


displaced to P’ of coordinates 


Pere (2.2) 
n=yrt. 
In this transformation the material around P undergoes a rotation 
1 (# ou) 
= se 2.3 
oa 2 \ox oy (2.3) 
and a deformation 
au ov 1 (# au) 
-~ = = » = =i — —j}. 2.4 
wo Ox “ow Oy fev = 9 \ar + oy G4 


In order to deal with the problem of relating stress to strain, it is important to consider 
the phenomenon from the view point of an observer rotating with the material through 
an angle w. If we denote by E;; the strain components referred to rotated axes, then 
it is clear that to the first order these are the same as e,; defined above, 


Ei; = €,; « (2.5) 


"69 a7 


This is not true, however, for the incremental stress. If we refer the stress to axes rotating 
with the material, the total stress acting on the element is 


: Y ) 
eget. (2.6) 
So, + 82; vs — Soo 


The incremental stress s,; is referred to rotated axes and depends only on the strain 
e,; . The nature of this relation is a physical problem which we shall consider in the 
next section. We shall deal here only with the geometrical aspects of the problem and 
established equations for the incremental stress field s;; which express the condition 
of equilibrium of the field. This is done by first referring the stress field to the unrotated 
initial directions x, y. These stress components are to the first order 
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ca, = Si + $i1 ieee 28,20, 
Oy, = Soo + S22 + 28,20, (2.7) 
Or, = Sie + Si2 + (Si ia Soo)w. 


Note that in the deformation this stress field is the one at point = x + u andy = y+ a. 
Also that the mass density has now become p’ by the incremental volume change.* The 
equations of equilibrium of the field o;; are 


Oo. Oc-, eee 
— + — + eX = 0, 
0& On 
(2.8) 
0o., 00,1, i 
——— oe og Y = 0. 
0& On 


We now express these equations by means of the original coordinates xz, y. Using the 


relations 
4] 4] 
dt = (1 + au) dx + dy 
O21 oy (2.9) 
av Ov 
1 = a +(1 2) 
Teer + oy 
and solving for dz, dy, we find the partial derivatives of x, y with respect to &,n, e.g. 
Ox ] ( 2) 
—= l —~I ete. (2.10 
a" BY ay * ) 


D is the Jacobian, i.e., the determinant of the transformation (2.9). With these ex- 
pressions we may write 

0022 _ 90,, OX Oo., OY _ 1 Bees ( a) 1 9022 Ov (2.11) 

0é Ox dé Oy ot D dx oy D dy dx 
Conservation of mass requires 

Dp’ = p. (2.12) 

The equilibrium equations (2.8) may then be expressed in terms of x and y. Substituting 
the values (2.7) for o,;; retaining only first order terms, gives 


iy ¥ ~ 3 = (Sia) > - (Si, — S22)o] 
, (2.13) 
OSi; Ov = OS,, dv _ OS. Ou OSi2 Ou -_ 
Ox Oy Oy Ox Ox Oy Oy Ox 
Another equation is obtained by permutation of the symbols and changing w in — w. 
These equations are further simplified by introducing the following identities 
; (2.14) 
Ou 
ay = €,,— Ww 


*For simplicity we assume here that the body force is independent of the coordinates. The more 


general case is considered in [4]. 





1959] FOLDING OF A LAYERED VISCOELASTIC MEDIUM 189 


and by taking into account the equilibrium condition (2.1) for the initial stress field. 
Equations (2.13) are then written 


AS 1 AS 9 - ’ a) ’ Y 
= + — + p¥w — 28y, 5" + (Si — Su) oe 
03 oy c y (2.15) 
aS, (28 dSi> dS,> 
ue — (2 4+ = <—2,.=0. 
Ton ou ay + ax + ay 


The other equation is also obtained by permutation of the variables and changing 
w in — w. If the initial field is uniform, the terms written on the second line disappear 
and the equations become 
Os . 0 eo 
ax + —_ + pYw = 2812 = + (S,, — Soo) a = 0, 
” . (2.16) 


08), 





OS; O8o0 ” Ow , , OW 
ean —— = y. S y we S = Ben — = e 
ox + ay pXw + 282 ay + ( 11 22) ar 0 


It is interesting to see that in the absence of a body force the only additional terms 
due to the pre-stress depend on the total initial shear and disappear if the initial stress 


is hydrostatic. 
The boundary conditions are found from the stresses o;; . The force acting on an 


arc element dé, dn, after deformation is 
dF, = Cop dn —~ Oxy dé, (2.17) 
dF, = a,,dn — oy, dé. 


The force F is acting on matter lying to the left side of the arc element dé, dy in a counter- 
clockwise coordinate system and to the right in a clockwise system. Substituting the 
values (2.7) for the stresses, expressions (2.9) for dt, dy and retaining only first order 


terms, we find 
dF, _ — [812 + Sie a Sow + S202: ‘ini Sii€zy] dz 
+ [S11 + 8u — Siw + Sii€vy seit Sj 2€2y] dy, (2.18) 


dF, = —[822 + Soo + Siow + Soe, — Sizesy] dx 
+ [8:2 + Si2 + S,0 + Sye, — So2€z2y] dy. 

This gives the force acting on an element of arc initially of components dz, dy, at point 
2, Yy. 
3. The viscoelastic half space under combined surface loading and lateral com- 
pression. We shall apply the previous equations to the problem of determining the de- 
flection of the surface of a viscoelastic half space under the action of a load normal to 
the surface and a condition of prestress which is a compression in a direction parallel 
to the surface (Fig. 1). 

The surface is along the z-axis and the y-axis is directed inward. The initial stress 
system is assumed to be 
So = Siz = 0, Su = —P, (3.1) 
ie., a uniform compression along the z-direction. We assume the body force to be zero. 


Equations (2.16) become 
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O81, 08,2 » Ow 0 
Ox On ay Sti‘ ’ 
, . (3.2) 
OS; O85 ‘ 
ed a P — a ® 
Ox oy Ox 


The next step is now to introduce a relation between the incremental stresses s;; and the 
incremental strain e;; . Such relations have not yet been developed in the most general 
case for a pre-stressed material, and some difficulty arises in this connection. We have 
shown [4] that for an elastic body under prestress, the elastic constants cannot be the 
same as in the unstressed condition unless the prestress system is hydrostatic. We there- 
fore point out that, in general, it must be considered an approximation to use—as we 
shall do here—the same stress-strain relations for the prestressed and the initially 
stress-free medium. We have discussed this point more extensively in [2]. 


! G =4 cos Ix 





+ 
Vo = cos Ix 
= o= Vv —_—— 
P — = a— P2z2-Sy 
y 
_—_ = 


Fic. 1. Viscoelastic half space under combined surface loads and lateral compression. 


Another point arises if the body is viscous for slow deformations. Then the initial 
state is one of steady strain rate, and there are additional incremental strains which are 
not due to the incremental stresses s;; alone. 

In the applications below we shall, therefore, consider the total displacement 
field as u; + u* where u* is the displacement associated with the steady strain rate due 
to the initial stress S;; . The term u,; then represents the departure of the displacement 
from the initial uniform strain rate. The incremental strain and rotation appearing in 
equations (2.15), (2.16) and (2.18) should then be written as e;; + e* and w + w* where 
e* and w* correspond to the initial steady state. 

In the applications below, the components e* and w* drop out of the equations 
and do not appear in the formulation of the problem. We may, therefore, disregard the 
term u* in the present treatment. In general, of course, since e* and w* are linear func- 
tions of time, they will disappear from the equations by introducing second time de- 
rivatives of all variables. 

The operational stress-strain relations for an initially stress-free medium were 
established in [7] and found to be 
8,; = 2Qe,;; + 6,,Re (3.3) 


with the unit matrix 6,; and e = e,, + @,,. 
The operators are 


Q =p [ - a) y(r) dr + Q + Q’p, (3.4) 
70 Pp ae 

- r*“ Rr ; oe 

Kk = p | > + : y(r)dr+R+R P; (3.5) 








oe 


Sar PT 
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where p designates the time differential operator 
d 
,=—: 
F dt 
By a well-known property of such operators, the relations are also valid if p is a real or 
a complex quantity, provided all variables contain the exponential factor e”’. 
Substituting expressions (3.3) for the stress into the equations (3.2), we find 


(3.6) 


OVin+ 16+ RIX - PS =0, 

Ee , (3.7) 
Vv 5 4 py) 2C _ poe _ 
QV0+ G+ RIS - PS =O. 


It can be verified that solutions of these equations with a sinusoidal factor along z are 
u = —sin lr[Ae~'” + Ck(1 — B)e™"*”), 
| (3.8) 
v = —cos lr[Ae'” + C(1 — Bk’)e™"*”]. 
These solutions are analogous to those developed in [3, 6] for the purely elastic case. In 
these expressions A and C are unknown operators, and we have put 
2) + P 1— ¢\'” P 
g=—t, ga=(—S), gar. (3.9) 
20+R 1+ 20 
If k is chosen positive the solutions (3.8) represent deformations confined to a region 
near the surface. The same is true if k is complex and chosen such that its real part is 
positive. 

Consider now the case where a normal load of sinusoidal distribution is applied to 
the surface (y = 0). We wish to determine the deflection at the surface due to the simul- 
taneous action of the surface load and the compressive preload P in the material. The 
surface load will be represented as 


an = q = q cos lx. (3.10) 


It is acting downward in the y direction as indicated in Fig. 1. The tangential component 
of the boundary force vanishes 
dF, = 0. (3.11) 


Introducing expressions (2.18) into the relations (3.10) (3.11), we obtain the boundary 
conditions 
¢., = @, 


8.2 = —@ cos I, = 20e,, + Re. (3.12) 
Substituting in these conditions the values u and v from expressions (3.8) and putting 
y = 0 yields 
Y . 


by : - 
A+-5(1 +k — 26k’) = 0, 


CR . q ’ 
— kal — k*) = -—-—ti- 3.13) 
20 B ( 


A + Ck(1 — Bk’) + 501 
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The conditions determine the operators A and C. The deflection of the surface is also 


sinusoidal and represented by 


Vo = V cos iz. (3.14) 

The ratio of the load q to the deflection v) at the surface y = 0 is then given by 
 @€. te kQ + ¢)° — | - er 
oe B| ¢ — af (3.15) 


with the notation 
10(Q + R) 
20+R ’ 
+h (3.16) 
Q 
Y= a — © 
Q+R 


The effect of the compressive load P arises through the bracketed factor in Eq. 
(3.15). This factor is 


= : [AC -f. ¢)° 7 1] — at. (3.17) 
$ 


The deflection under a load q is 


V = —®. 3.18) 
iB 7 no 
If a constant sinusoidally distributed load is suddenly applied at time t = 0 with a 
maximum value unity, 
0 ft 0 , 
% = l(t) = . (3.19) 
ll t>0 
it may be represented by the complex integral 
a= b= = f< ‘ (3.20) 
o = \ = te ap. 0.4 
2at Jc-io P I 
Hence the deflection under this load is 
4 Oe [ co (3.21) 
== a 1, 3.2 
} tlt Jc-i2 pBey P 


In this expression B and ¢ are functions of p. 
The effect of the lateral compression P appears through the factor ¢. If this lateral 


compression is zero, then 


¢ = 0, g= 1 (3.22) 
and the deflection is 
1 C+ia e”! 
i) == — dp. (3.2% 
V(d) 7 [ . pb 3.23) 


In order to give an idea of the magnitudes involved, let us consider first the case of a 
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purely elastic incompressible material. Incompressibility corresponds to 


R=o0 or a=0. (3.24) 

With a shear modulus G, we put 
Q=G. (3.25) 

Hence 

.¢ Ste 2 

- c (3.26) 
with 

¢ = P/. 


The deflection is amplified by the lateral compression P, and the amplification 
factor is 1/g. Values of this amplification factor as a function of ¢ are given in Table 1. 


TABLE 1. 

§ l/¢g 
0.50 1.67 
0.60 2.14 
0.70 3.29 
0.80 10.00 
0.85 a 


We see that when ¢ = 0.85 i.e. for P = 1.7G the amplification factor is infinite and 
the surface becomes unstable. We also notice that the amplification factor becomes 
appreciable only for values ¢ > 0.5; hence, for compressive loads P larger than G these 
loads are of the order of the elastic modulus. In a material of constant modulus, they 
ean only occur for loads corresponding to a very large compressive strain. However, 
we must remember that the present theory is for incremental deformations. Hence, the 
modulus G appearing in the above expression is actually the so-called tangent modulus. 
In many cases where plastic deformations are involved, this modulus decreases with the 
load, and the above values may well fall within observable ranges. 

In regard to the significance of the amplification factor, it should be remarked that 
it is the same for all wave lengths. Hence, the deflection of the surface under an arbitrary 
load distribution may be derived from the deflection under the same surface load and 
zero compression (P = 0) by simply multiplying the latter by the amplification factor. 
The compressive load introduces no distortion of the deflection but simply a change of 
magnitude. 

It can easily be seen that the amplification effect of the compressive load on the 
surface deflection of a half space also applies when the initial deflection is not due to 
a surface loading but due to irregularities and waviness. This departure from perfect 
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smoothness is amplified without distortion by the compressive loads. This suggests 
that under high compressive loads, surfaces which appear to be smooth initially will 
tend to wrinkle. Since the amplification factor is dependent only on the tangent modulus, 
this mechanism suggests itself particularly as the explanation for the wrinkling which 
appears on surfaces of materials compressed in the plastic range. 

4. Instability of the free surface of a half space under compression. Let us now 
examine the condition of instability of the free surface of a half space under compression. 
This corresponds to the existence of solutions different from zero for the deflection V 
while the surface load gq is zero. Equation (3.18) may be written 

do = $BeVI. (4.1) 


The condition gq = 0 with V ¥ 0 amounts to 


] 2 
g=-(k1l+9* -—1] -—eat =0 (4.2) 
$ 
which may also be written 
kd + 0° 1 _ 
§ § 
or 
2— 2° — & = Wt t+a’e’. (4.4) 


Let us first examine the elastic case. The operators Q and FR are replaced by the 
Lamé constants 
Qg=G kh = X. (4.5) 


Poisson’s ratio being designated by »v we find 


G 
— 2». (4.6) 


Ao 


Also 
¢ = P/Q. (4.7) 
The critical value of ¢ is the real root of Eq. (4.4). This root depends on Poisson’s 


9 


ratio only and its value is shown in Table 2. 


TABLE 2 
v ¢ 
0 .55 
1 .68 
i 84 





The instability occurs at values of the compressive load P comparable with the 
shear modulus G. However, as already pointed out above, since the latter is the tangent 
modulus, its value may be low enough in the plastic range for the surface instability 


to occur within the practical range of P. 
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As an example of a viscous material, let us now consider the case of a viscous in- 
compressible fluid. Some clarification is needed here as to what is meant by instability. 
We assume that the half space is initially under the influence of a compressive load P 
parallel with the free surface and that the fluid is initially in a steady state of uniform 
rate of deformation under this load. All velocities and strain rates are initially constant 
in time and space. The stability of the free surface under those conditions is determined 
by considering the possible exponential growth of a disturbance starting from a certain 
instantaneous configuration considered as the initial state.* The possibility of such 
instability is easily determined from the above analysis. For an incompressible fluid, 


the operators become 


Q=up, R=, (4.8) 
with » the viscosity coefficient. Since the condition R = © implies a = 0, Eq. (4.4) 
becomes 
2-—2°°-¢ =0 (4.9) 
with 
ne eee 
$20 ~ Qup 


Equation (4.9) is a cubie in the unknown p. For values of p which are roots of this 
equation, there exist disturbances proportional to the factor e”’. Instability corresponds 
to roots which are real and positive or complex with a positive real part. The three roots 
p of (4.9) are 


P 
pi = 1.192 Ou , 
_— me 
Pp. = (—0.596 + 0.2552) 2, ? (4.10) 
n 


> 
ps = (—0.596 — 0.255%) 2. 
Qu 


When substituted in Eq. (4.2), i.e. ¢ = 0, the complex roots do not furnish a solution 
if we choose k to have a positive real part. Hence they do not correspond to the physical 
problem. 

The root p, corresponds to instability. In order to appreciate the magnitude of this 
instability, let us consider a time ¢, during which the undisturbed uniform flow produces 
a 25% compressive strain. The rate of strain in the z direction under the load P is 


hg EE es (4.11) 
t; 4u 
Hence, 
t == pit, = 0.596. (4.12) 


,” 
*In dealing with a viscous fluid whose initial state of stress gives rise to a steady strain rate, we must 
refer to the remark Sec. 3 by which it is justified to consider not the actual incremental deformation, 
but the departure of the deformation from the steady state itself. 
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This means that in the time interval ¢, during which the material is squeezed by 25%, 
any disturbance of the free surface from a perfect plane is amplified by a factor e”*** 
1.81. This instability is of such a mild form that it will be very difficult to observe in 


practice.* 

A similar procedure may be used for any kind of viscoelastic half space by solving 
Eq. (4.4) as an equation for p with the operators Q and R as appropriate functions of 
the unknown p. 

5. Characteristic equation for the stability of the embedded layer. We consider a 
layer of thickness h embedded in an infinite medium. The layer and the surrounding 
medium have different viscoelastic properties. The layer is subject to an initial com- 
pression P and the surrounding medium to a compression P, (Fig. 2). This initial state 
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Fic. 2. Layer embedded in infinite medium under lateral compression. 


of stress will depend on the type of initial deformation imposed on the system and on 
the viscoelastic parameters of the two media. For example, if the two media are in- 
compressible viscous fluids and if the system is squeezed uniformly in a direction parallel 


with the layer, the compressive stresses are in the ratio: 


P 
c (5.1) 


where uy is the viscosity coefficient of the layer and y, that of the surrounding medium. 
The stress could also be due to other causes, as, for instance, an increase of volume of 
the layer while longitudinal elongation is prevented. In this case the stress P, in the 
surrounding medium can be zero. 

In order to simplify the formulation, we shall assume that the interface between 
the layer and the medium is perfectly lubricated so that there are no tangential stresses 
at the interface. A previously developed approximate theory for the case of perfect and 
imperfect adherence [2] indicates that this assumption does not modify the result sub- 
stantially. The properties of the layer are represented by the operators Q and R and 
that of the surrounding medium by Q, and R, . We consider the layer separately and 
locate the z-axis in the plane equidistant from the two faces (Fig. 3). The y-axis is normal 
to the layer. Flexure of the layer is represented by antisymmetric solutions of Eqs. 


*Since the theory is only valid for small deformations, its application to the case of a large com- 
pressive strain of 25 per cent will, of course, give only an indication of the orders of magnitude of the 


amplification effect. 
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Fic. 3. Forces acting on the layer. 


(3.7). They are obtained by adding two solutions of the type (3.8) by using both positive 
and negative exponents. We find 

u = —sin lz[A sinh ly + Ck(1 — 8) sinh lky], 

v = cos lx[A cosh ly + C(1 — Bk’) cosh lky], (5.2) 
where 6 and k are defined by (3.9). 

This also coincides with the solution obtained previously for the purely elastic case 
[3]. The normal force applied to the layer at the boundary y = h/2 and due to the 
surrounding medium is 

q = Qo cos lx. (5.3) 

This load is taken positive in the negative y direction, and the tangential force vanishes 

owing to the condition of perfect slip at the interface. Because of symmetry, boundary 


conditions imposed at one interface are automatically satisfied at the other. We shall 
consider the interface at y = h/2. Applying (2.18) we find the boundary conditions 


dF, h 
== —s,, — Pe,, = 0 y= - 
dx . 2 ” 
(5.4) 
dF, i 
= —S8,. = —p cos lx. 
dx = d 
From expressions (3.3) for s,. and s.. in terms of the strains we find at y = h/2, 
e., = 0 ‘ 
r . (5.5) 
Qe,, + Res + Cy) = —o cos Ix. 


Substituting the expressions (5.2) for uw and v yields two equations for the operators A 
and C in terms of g . We must also consider the boundary value of the normal displace- 
ment v. This is expressed by substituting y = h/2 in the second equation (5.2). We find 


v= Veoslx = E cosh . + C(1 — Bk’) cosh a cos lx. (5.6) 


The normal deflection at the boundary y = h/2 is V cos lx. By substituting in (5.6) the 

expressions for A and C determined from the two equations (5.5), we find a relation 

between the luad g and the interface deflection V. This relation may be written 
(1+ ¢) 


9 | - 2 
So = +2) 4, an Seal. b tanks b x 
VB ¢ + af’) tanh P k tanh ky, (5.7) 
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where 


2 


20 ’ 


te 


y = 3lh, f= 


and B and a are operators depending solely on the viscoelastic properties of the layer 
__9 
sSihed (5.8) 
_ 40@Q + R)_ 
20+R 


Q@ 


At this point we still have to match boundary conditions for the layer and the half 
space representing the embedding medium. This is readily done if we consider the ratio 
of the load gq to deflection V for the surface of the half space which was obtained in 
Sec. 3 and given by expression 


== A zs 
“i ams (5.9) 
where 
r ) / J = 
B. == a Ee 
me (5.10) 
2 — 261 wat o a 
ik = = “EDs 


1 84 EY 4 1 


ee 3 _ (= 4)" _~_%_. 
a ao eae a 2 ee 

The quantities B, and a, depend only on the viscoelastic properties of the surrounding 
medium, while the others include the value of the prestress P; in this medium. This 
prestress is contained in the expression g, and, as was pointed out above, we may con- 
sider 1/y, as an amplification factor for the deflection at the surface of a half space. 
The value of the amplification factor is unity when the prestress is zero (P,; = ¢, = 0). 
Since the load g, and the normal deflection V are the same for both the layer and the 
medium at the interface matching of the boundary, conditions amount to equating the 
values of g,/l1V appearing on the left side of the two equations (5.7) and (5.9). This 
yields the characteristic equation 


l 2 ( ¢)? 
— g, = —(1 + af’) tanhy — = k tanh ky. (5.11) 
é 


It is true that the values of g, and V for the layer and the medium are not those of 
the same location of the interface because of the relative slip of the two media. However, 
the error is of the second order and is irrelevant in a first order theory. 

The parameters of the surrounding medium are contained solely in B,g, on the left 
hand side. The significance of this characteristic equation lies in the fact that for given 
materials and a given prestress field P and P, , it constitutes a relation between the non- 
dimensional wave number y = 3lh of the folding and the value of p. This value of p is 
the coefficient in the time factor e” which multiplies the amplitude of the folding for 
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any particular wave number y initially present as a disturbance in the layer. In a system 
with random initial disturbance, those wave lengths will appear for which p has the 
largest positive value. We shall disregard the possibility that there are complex values 
for p as encountered above in the case of the viscoelastic half space and assume that 
such complex values are not physically significant. 

6. Numerical results and discussion. An approximate expression for the charac- 
teristic equation (5.11) is obtained if we expand the hyperbolic tangents in power series 
to the third power in y. We find 


B . , 
R ¢, = fy(1 +a) — Fy [2 — (1 — ade]. (6.1) 

If the wave length of the folding is sufficiently large, y and ¢ are small. We may 
further neglect (1 — a)t in the bracket. Also we may assume ¢, to be small, which 
amounts to neglecting the prestress P, in the surrounding medium. Hence we may 
put g, = 1. Equation (6.1) then reduces to 


3 
St = gl +a) — fy’. (6.2) 
This may also be written 
l so 
P= 12 B(th)? + Th (6.3) 


This last equation coincides with Eq. (3.1) of [1] obtained by neglecting the prestress 
P, in the surrounding medium and using a “thin plate” type theory for the layer. A 
complete discussion of the approximate equation (6.3) was given in [1] for the general 
case of two media of arbitrary viscoelastic properties and for a certain number of par- 
ticular cases. 

The exact equation (5.11) provides a means to test the accuracy of the approximate 
theory. For this purpose it is sufficient to compare the two theories on a simple example. 
We shall choose the case of two incompressible fluids. Incompressible fluids are repre- 


sented by 


R=R, = @, 

om m= 9, (6.4) 
Q = up, 
Q, = MiP, 


where u is the viscosity coefficient of the layer and y, that of the surrounding medium. 
Introducing these expressions into the exact characteristic equation (5.11), we derive 





Mi = 1 tanh y— (i+ 9° k tanh ky, (6.5) 
u t 
where 
. £- i= i. 
MO” kl +o)? #1? i 
. P % P, 
sats Qup —_ Quip 
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There are two types of approximations involved in Eq. (6.2): that due to the use of 
thin plate theory and the other due to neglecting the prestress P, of the surrounding 
medium. Let us discuss the effect of each separately. Therefore we put g, = 1 in Eq. 
(6.5), i.e., we consider first the error due to the use of plate theory alone. The character- 


istic equation is then 

(1 y-\ 2 a= 
——_$/ k tanh ky. (6.7) 
§ 


I 
Fi - tanh y — 
$ 


This is a relation between these variables representing a function ¢ of y with a param- 
eter u/u, . The variable ¢ versus y is plotted in Figs. (4) and (5) for six values of n/p, . 
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Fic. 4. Plot of ¢ = P/2up versus y and the parameter u/y; from the exact relation (6.7) and the ap- 
proximate solution (6.9). 


p/m, = 6, 12, 36, 72, 144, 288. 


Consider now the approximate equation (6.2) which corresponds to the thin plate 
theory and, also, to neglect of the prestress P, . For incompressible fluids, this equation 
becomes 


(6.8) 
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or 


pa thee (6.9) 

Yo 
This also represents a family of functions ¢ of y which are plotted as dotted lines in Figs. 
(4) and (5) for the same values of u/u, . Examination of the two plots affords an im- 
mediate comparison between the exact and approximate theories. Of particular interest 
is the minimum value of ¢ denoted by ¢,,;, and the value y at which this minimum occurs. 
It is easily seen that the minimum of ¢ for a given value of the compressive load P 
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Fic. 5. Same as Fig. 4. 


corresponds to maximum values of p. The significance of a root p of Eq. (6.7) is given 
by the fact that for such roots a disturbance of wave number / of the layer increases 
proportionally to the factor e”’. Hence, maximum values of p correspond to wave lengths 
of maximum rate of growth. We have referred to this wave length as the dominant wave 
length. The wave number 1, of this dominant wave length is therefore determined by 


h 


6.10) 
a - 


ve = lh=e 


(Lz = dominant wave length). 
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From formula (6.9) we find the approximate values 


3 by 
Teo es , (6.11) 
4p 
Snin = (3 mJ" ; (6.12) 
4p 


The approximate value (6.11) of y, as a function of (u,/u)'”* is shown in Fig. 6 by the 
dotted straight line, while the exact value of y, [determined from Eq. (6.7)] is represented 
by the full line. Comparison between the two curves indicates that a discrepancy in the 
value of y, begins to be noticeable for values (u,/u)'’* > 1/4 or u/u, < 64. The error in 
the dominant wave length may rise to 10-15% when u/p, becomes of the order of ten. 
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Fic. 6. Value y, = zh/La (Le = dominant wavelength) as a function of u/y: by exact and approximate 
theory. 


We have also plotted the exact value of ¢,,;, versus (u,/u)** as determined from Eq. 
(6.7). It is represented by the full line in Fig. 7. The approximate value (6.12) is repre- 
sented by a dotted straight line. The error becomes noticeable for (u,/u)** > 1/16; 
hence, again for u/u, < 64 For smaller values of u/y, , the error may reach 10-15%. 

Let us now look upon the other approximation involved in Eq. (6.9) and which is 
due to neglect of the prestress P, in the surrounding medium. In order to evaluate the 
effect of this prestress, we must reintroduce the factor yg, which we had put equal to 
unity. This amounts to using the same graphs as in Figs. 4 to 7, but with a corrected 
value y,¢, instead of yu, for the viscosity of the surrounding medium. Let us evaluate 
the value of the correction factor ¢g, in the vicinity of the dominant wave number (near 
Smin)- Assuming the layered material to be squeezed uniformly so that the strains are 
the same in the layer and the surrounding medium the compressive prestresses P, and 
P satisfy relation (5.1) and we derive 


t= ¢, (6.13) 


We may, therefore, evaluate ¢g, as a function of ¢. Values of yg, are given in Table 3. 
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TABLE 3. 


c ¥1 | (¢,)3 gis 

1 946 | .980 | .962 
3 804 930 | .865 
5 


598 841 | 710 


The values of ¢, show a decrease in the “effective value”’ u,¢, of the viscosity coefficient 
of the surrounding medium due to the prestress. The decrease is about 5% or less for 
values of f.:, < 0.1, ie., for u/u, > 64 and may be as high as 40% for the lower values 
of u/u, . The effect on the dominant wave length is much smaller and is proportional to 
(v,)'’*. From the values of (g,)'”* in Table 3 and from Fig. 6, it can be seen that the domi- 
nant wave length will be increased from 2 to 16 per cent in the same range. Reference to 
Fig. 7 and values of (y,)*”* show that ¢,,;, Will be decreased by an amount somewhat larger 
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Fic. 7. The minimum value of ¢ as a function of u/m by exact and approximate theory. 


due to the prestress. In general, we may conclude that the effect of the prestress P, in 
the embedding medium becomes negligible if the viscosity ratio u/u, of the layer to the 
medium is larger than about 64. It is interesting to note that this was also found to be 
the condition for the error, due to thin plate theory, to be negligible. 

In order to complete the present discussion, we must also pay attention to the degree 
of instability involved in the phenomenon of folding. The dominant wave length de- 
termined by y, is that for which an initial waviness in the layer increases at the fastest 
rate. After a time ¢, the amplitude is multiplied by the factor e”’ where p is given by 





P 
= . 6.14) 
P 2uS min ( 


We may ask, for instance, what is the value of the factor e” after a time ¢, such that 
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the layered material is squeezed by 25 per cent under the compressive stresses P and P, . 
Referring to the same argument as in Sec. 4 in the case of the half space, we find 


t, = ‘ (6.15) 
and 
pl, a l 2S min . 


From Fig. 7 giving {mi, aS a function of u/u, , we may derive the amplification factor 
3 ry e rn 
e”’* in terms of u/u, . Values are shown in Table 4. 


TABLE 4. 


Ms M1 exp (ph) 
: 3.78 
27 14.8 
64 127 
125 1940 


The amplification factor increases sharply and becomes quite large for values 
u/u, > 64. Below these values the magnitude of the amplification factor becomes 
relatively small, indicating that the instability in that range is not too significant and 
will not exhibit very sharp features. This lower range of values is also the one in which 
the errors become appreciable if we apply thin plate theory and neglect the prestress in 
the surrounding medium. We may conclude that in the region of significant instability 
Eq. (6.2) of the approximate theory is applicable. 

The present discussion is limited to the case of incompressible viscous fluids. A similar 
discussion can be carried out for any two materials of arbitrary viscoelastic properties 
represented by four suitable operators. If the two materials are incompressible, the 
plots in Figs. 4 to 7 represent master graphs by which the folding may be analyzed 
quantitatively by means of very simple calculations. 
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THERMAL SHOCK IN AN ELASTIC BODY WITH A SPHERICAL CAVITY* 


BY 
ELI STERNBERG anv J. G. CHAKRAVORTY 


Brown University 


Summary. This investigation aims at the dynamic thermoelastic response of an 
infinite medium with a spherical cavity to a sudden uniform change in the temperature 
of its internal boundary. By means of the Laplace transform, a closed solution to this 
problem—exact within classical elastokinetics—is obtained in terms of error functions 
of real and complex arguments. The ensuing temperature stresses are compared with 
the corresponding quasi-static results. 

Introduction. The conventional quasi-static treatment of transient thermoelastic 
problems rests on the assumption, apparently introduced first by Duhamel [1]', that the 
inertia terms may be neglected in the governing field equations. The quality of the 
approximate solutions thus obtained evidently depends both upon the magnitude of 
the time-gradients inherent in the temperature field and upon the size of the relevant 
intrinsic inertia parameters. 

Duhamel’s hypothesis has been the object of several recent studies. The earliest 
investigation of this kind is due to Danilovskaya [2] (1950), who determined the dynamic 
thermal stresses induced in an elastic half-space by a sudden uniform heating (or cooling) 
of the entire boundary, if the body is constrained to a uniaxial motion perpendicular 
to the bounding plane*. Later on, Danilovskaya [4] generalized her solution to accommo- 
date convective boundary conditions. Still more recently, the present authors [5] estab- 
lished explicitly the thermal displacements associated with the stresses given in [2], [3] 
and extended the completed solution to the case of a gradual (ramp-type) change of 
the surface-temperature. Related studies in the theory of beams and plates were carried 
out by Boley [6] and by Boley and Barber [7]. Finally, we mention thermoelastic in- 
vestigations by Nowacki [8] and Ignaczak [9] in which the inertia terms are taken into 
account, although neither of these papers is concerned with a quantitative assessment 
of the dynamic effects thus arising. 

As was pointed out in [5], the problem treated in [2], [3], [5] suffers from two deficiencies 
as a vehicle for examining the adequacy of the traditional quasi-static methods in 
transient thermoelasticity. First, the artificial constraint of the medium imposed in 
[2] results in a severe degeneracy of the corresponding quasi-static stress distribution, 
the normal stress perpendicular to the boundary being identically zero. Second, the 
absence of a characteristic length in this particular problem causes the pertinent inertia 
parameter to be hidden in the dimensionless time and distance scales whence its influence 
upon the dynamic effects sought, is obscured. 

lor the preceding reasons it appears desirable to examine the inertia effects en- 


*Received July 1, 1958. The results communicated in this paper were obtained in the course of an 
investigation conducted under Contract Nonr 562(25) with the Office of Naval Research, Department 
of the Navy, Washington, D. C. 

‘Numbers in brackets refer to the list of references at the end of the paper. 

*The results contained in [2] were independently reached by Mura [3] (1952). 
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countered in a less degenerate thermoelastic space problem. This leads us to consider 
an infinite, homogeneous and isotropic, elastic medium with a spherical cavity, the 
boundary of which is exposed to an instantaneous change of temperature. We assume 
that the body, which is free from loading, is initially at rest, and seek the ensuing 
thermal stresses and displacements if the temperature field obeys the (uncoupled) 
heat-conduction equation. The corresponding quasi-static solution was given and dis- 
cussed in [10]. It should be mentioned that the present dynamic thermoelastic problem 
is related to an ordinary problem of elastokinetics, i.e. to the determination of the stresses 
produced by a sudden uniform pressure applied to the surface of a spherical cavity in 
an infinite elastic medium. An exact closed solution to this more elementary problem 
was obtained first by Jeffreys [11] (1931) and has since been rediscovered by various 
subsequent authors. 

Formulation of problem. Dimensionless variables. Let (r, 6, ¢) respectively denote 
the radial coordinate, the co-latitude, and the longitude of a spherical coordinate system. 
Suppose that the medium under consideration occupies the region a < r < ©, where 
a is the radius of the spherical cavity. The thermoelastic problem described in the 
Introduction is characterized by polar symmetry about the origin, whence 


2 = 20. 2 u, = ulr, 0), ue =u, = 0, (1) 


in which ¢ is the time, 7 the temperature, while (u, , us , u,) designate the spherical 


components of the displacement vector. 
If 7, is the constant temperature suddenly assumed by the internal boundary r = a, 


and x stands for the thermal diffusivity of the material, we may introduce a dimensionless 
radial coordinate, time and temperature, by means of 


. i (2) 


With reference to these dimensionless variables, the independent temperature problem 
is governed by the (uncoupled) heat-conduction equation 


ay 2 dy = de (3) 
dp pop Or’ ‘ 


subject to the initial condition 
g(p, 0) = 0, (4) 
together with the boundary and regularity conditions 
g(1, 7) = h(n), g(p,7) 70 as poo, (5) 
where h(7) is the Heaviside step-function, given by 


h(r) = 0 for —2e7 <7r<0, h(r) = 1 for O<r< o. (6) 


We turn next to the characterization of the thermal displacements and stresses 


associated with the temperature field g(p, r). To this end, let a dimensionless radial 


displacement u*, as well as dimensionless spherical components of normal stress o* 


and o% , be defined through 
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1959] 
l—p } 
es —————— si, 
, (1 + v)aTya (7) 
l-—vp l-—vp 


O we ction cee sie 
= ’ 0, 7 = mn a) 
21+ vaTymu "’ 21 + vial yp ~’) 
in which o, , ¢» are the corresponding physical components of normal stress, whereas 
u, v, and @ denote the shear modulus, Poisson’s ratio, and the coefficient of thermal 
expansion, respectively. In view of (1), the spherical components of shear stress vanish 


identically and c, = o» . Finally, it is expedient to introduce a dimensionless inertia 
parameter y, defined by 


~ (8) 


where c is the velocity of irrotational waves, 


2(1 as v) 





c 


8 being the mass density of the material. 
The displacement equation of motion now takes the form 
2, % * * 2 Fs) 
ou 2 du 2u 20u* , dp (10) 


ap" p ap Z alt ar ap’ 
while the stress-displacement relations become 


ou* 


* 
(1 ~ Met = (1-9) = 4+ H—= 0 — de, 
Op p 


(11) 


, ju* ut 
(1 — 2Qv)o% =v = + a (1 — vg. 
p p 


Once ¢(p, r) is known, the solution of the thermoelastic problem at hand reduces to 
the determination of a function u*(p,r) which satisfies (10), meets the initial conditions 


* 
u*(p, 0) = 0, | = @, (12) 
OT JSip,0) 
and is such that the stresses (11) conform to the boundary condition 
o*(1, r) = 0, (13) 


as well as to the regularity conditions 


*(p, r) 3 0, at(p, 7) 70 as poo. (14) 


We observe that the quasi-static formulation of the present problem®* is obtained by 
setting y equal to zero in (10) and by deleting (12). 

Introduction of the Laplace transform. Solution in the transform domain. The solution 
of the elementary heat-conduction problem governed by (3), (4), (5), is well known‘ 


and admits the closed representation 





‘See [10]. 
‘See, for example, [12], p. 209. 
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‘ l — |] ote 
oe,7) =-¥®, §=o oz (15) 
p 2(T) 
where 
= ° 
y(t) = erfe (§) = —p3 | exp (— 7°) dn (16) 
T ve 

is the complementary error function. For future convenience we recall* here that 

y¥(0) = 1, vii) -0 as E> om, ¥(—) = 2 — vy), (17) 


and cite the power-series development 
2 = (-1)"¢"" 
y(t) = 1 — a =——— (18) 
= mw” Sa0 (2n + 1) |n 


as well as the semi-convergent expansion, valid as Re(é) — + @, 


) N n 12n 
» exp (—£) a 1) F -2N-3 r . 
(¢ - - - - + < = - se0od ( 
s/ nw & Pp n (2¢)"" Og (A 1,2, , (19) 


With a view toward the determination of the thermal displacements and stresses 
arising from the temperature distribution (15), let 


U(p, s) = L{u*(p, r)} = | exp (—sr)u*(p, 7) dr, (20) 


/0 


so that U(p, s) is the Laplace transform with respect to r of u*(p, 7), s being the transform 
parameter. From (15), 


] 2 6 
®(p,s) = Lig(p, r)} = exp [—(p — 1)s'”’], (21) 
ps 


whence applying the transform to (10), and bearing in mind (12), we reach the ordinary 
differential equation 
@U , 2dU 2 = 1 ; 
“ie ie - ( a + 78°) == alk 35 + =) exp [—(p — 1)s'’*]. (22) 
up p dp p / ps ps 


Furthermore, by virtue of (11) to (14), (20), (21), we arrive at the transformed boundary 


and regularity conditions 


ca ee ae 0 “ef 
") dp ides p ' je! (23) 


U(p,s) — 0 as pO. 


Equation (22) is reducible to an inhomogeneous modified Bessel equation. Its general 


solution is given by 


U = V(p, 8) fe W(p, 8), (24) 


5See 13], p. 147. 


6See [14], p. 246. 
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where W(p», s) is a particular integral of the complete equation, and 


‘ As r2(5 ss ” 
(9,8) = Min Tals) + Me K, (80). (25) 


Here, \,(s), \2(s) are arbitrary functions of s, while J;,. and K3,2 designate modified 


Bessel functions of the first and second kind, respectively. As is readily confirmed, W 
may be taken in the form 


W(p,s) = mo sell exp [—(p — 1)s'”’]. (26) 


On the other hand, (23) to (26), after an elementary computation, yield 


Ho. a) = Lt rely’ + Iii + 1, ‘97 
V(p, 8) = ps(1 as v's) [(ys + p) + q] exp [ y(p 1)s], (27) 


provided 


9 


— 


- i. = 2Qv)"’ 


p<, go (28) 





Inverse transforms. Solution of problem. We now proceed to establish the desired 
inverse transform of U(p, s). For this purpose we note from (26) that 


W = W,(p, 3)W2{(p, 8), _ (29) 
if 
W, = p(y’s— 1)", W. = (8? + ps”) exp [—(o — 1)s'”]. (30) 
Moreover, according to known results’, 
1 ) 
wt(p, r) = L7'{W,(p, 8)} = a3 exp (7/7), 
YP 


(31) 
w*(p, 7) = L'{W.(p, 8)} = (po + 1)(r/m)'” exp (-—#) + L “= P wt) + rv), 
& J 


where £ and y/(£) are defined in (15), (16). From (29), (31), and the convolution theorem, 


follows 


w*(p, 7) = L™*{W(p,s)} = [ w*t(p, 7 — n)w%(p, n) dn. (32) 


70 


Substitution of (31) into (32), and repeated use of integration by parts, ultimately 


lead to 


] , ae 
w*(p, 7) = 2p [(p — 27 — 1 — 2y')Vv) (33) 


— 2(p + 1)(r/)'” exp (-—#) + v'dilp, 2) — veve(p, 7)]- 


The auxiliary functions y, , ¥2 appearing in (33) are accounted for by 


7See [14], No. 1, p. 229; No. 4, p. 245; No. 7, p. 246. 
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7? — | 
vi(o, 2) = o(¢ + X-) exp (4 + 2—+) 


y (34)° 





. ue m r) oe (+ a e—t) 
Y Y ef 


Our next task consists in finding the inverse Laplace transform of V(p, s) in (27). 
Evidently, 
v*(p, r) = L'{V(p, 8)} = h@L'{G(p, 8)} (35) 
with 
w= 7r— y(p— 1), (36) 
if h again denotes the Heaviside step-function given by (6), and G is the coefficient of 
the exponential function in (27). Moreover, a partial fraction decomposition of G(p, s) 


results in 


vqe(1 + 2yp + 2y’p) 


where ‘“‘Re’’ stands for “the real part of”, k is the complex parameter defined by 


G(p,s) = Se eS |  — Fo |, (37) 
8s—y¥ s-—k 


k _ mae at + 71q (38) 
¥ , 
and 

F(p) = —q(y + op) + 2[1 + ply — 2ye — p)]. (39) 
The inverse of G(p, s) is readily found with the aid of the available inverse transform 
( 3/2 cand : 9 1/2 : 
L~*s p= ro? — YrA'?7'”)] exp (Ar) — 2 (2) . (40)” 

s— Aj \ \r 


in which } is a constant. Thus, (35) and (37) imply 


; h(w) ) 
v*(p, 7) = _—+ Ec + p) exp (w 7) 
p 
, 347 1/2 2 3/2 1/2 (41)"° 
+A Re {F(p) ly ¥@ ‘*/y) exp @/y) — k-"Wlkw) “~ exp (kw) ; 


+ B exp (kw) + Clo + 2/x)'”?] + Di |. 


The letters A, B, C, D, appearing in (41), refer to the following auxiliary functions of 


vy and y: 


’Note that the arguments of the complementary error function y here are — + r3/y or — — r/y, as 


indicated. 
See [14], No. 2, p. 233. 
Here, as well as subsequently, k* stands for the principal value of this root. 
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{  ——— ae 
vq(l + 2yp + 2y'p) ’ (42) 





B an : + E + a. C = y a b, D _ y' e EF. 


From (24), (32), (35), we have 
y* = v*(p, 7) + w*(p, 7), (43) 


where v* and w* are given by (41) and (33). 
The thermal stresses associated with the displacement (43), are obtained by sub- 


stitution of (43) into (11). If 


the results of these lengthy computations become: 


_ h(w) 


o*% = ro G + 2yp(y + p)] exp w/y’) 
+ A Re {[F(p)(p + 2vp) — Evply’*Wo'”’/y) exp @/y’) 
+ [Ep — F(p)(2p + ykp)\k* 7 y(k'“w'”) exp (kw) 
+ B[F(p)(2p + ykp) — Ep] exp (kw) 


+ C[2pF(p) — Ep]lw + 2@/x)'””] (45) 


+ F(p)(2pD + ypC) — Dep} | 


~ ; [(2r + 277 +1 — pV) + Ap + 1)(r/x)'” exp (—#) 





a (£ + 7) le, 7) + vev(p, 7), J 


and, setting 


(46) 


2 
Se ae 
™ 1 — 2p p ; 
o%4 = ~ | me — yy + p)] exp @/y’) 
+ A Re {[F(p)(mp — y) — Eymphy*¥@"'”’/y) exp @/7’) 
+ [Emp — F(p)(ymkp — Ik??? W(k'w'””) exp (kw) 
+ BIF(p)(ymkp — 1) — Emp] exp (kw) 


— C[F(p) + Emp]lw + 2(w/x)'””] (47) 


+ F(p)(Cymp — D) — DEme| | 


| 


~ [(2r + 27 +1 + p)vO + Ap + 1)(r/n)”? exp (-#) 





— (y’ — mp’)¥,(p, 7) + veve(p, 7)]. 
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to 
to 


This completes the formal solution to the thermoelastic problem at hand. One con- 
firms by direct substitution, with the aid of (17), (18), (19), that the displacement (43) 
and the stresses (45), (47), indeed satisfy the field equations (10), (11) if w ¥ 0, and con- 
form to conditions (12), (13), (14). 

Properties of solution. Quasi-static and steady-state limits. Examining the structure 
of the solution obtained in the preceding section, we observe that the displacement 
(43) and the stresses (45), (47), each contain a term which has the step-function h(w) 
as a multiplier. If y > 0, these terms reflect the presence of a spherical shock-wave 
issuing from the boundary p = 1 and propagating outward with the velocity c 
of irrotational waves, given by (9). The remaining terms are diffusive in character and 
correspond to a signal which is instantaneously received throughout the entire medium. 

In accordance with (6), (36), we have 


tr = y(p — 1), (48) 
provided p and 7 momentarily designate the dimensionless radius of the wave-front 
and the dimensionless time at which this radius is attained. While the time and space 
derivatives of the radial displacement suffer discontinuities at r = y(p — 1), u* itself 
depends continuously on p and rfor1 < p< »,0 <7 < ™, as required by the physical 
continuity of the material. On the other hand, the stresses exhibit finite jump-discon- 
tinuities at the wave-front. Specifically, , 





, ' l—y 
o*[p, y(p — 1) +] — o*[p, y¥(o — 1) —-)] = = ; 
(1 — 2y)p 
(49) 
r v 
sip, (9 — 1) +] — oflp, (0 — 12) —-] = —— 
a%lp, v(p +] [p, v(p ] ie 
Thus, o% remains continuous only if y = 0; the jumps in both stresses are independent 


of the inertia parameter y and become unbounded as v — 1/2, if uw remains fixed in this 
limit. 

It is interesting to compare (49) with the stress-discontinuities generated in Jeffreys’ 
problem [11] by a radial pressure of magnitude o, , which is suddenly applied to the 
spherical boundary and is steadily maintained thereafter. Taking, in this instance, 


r = ct/a, o* = o,/00, 0% = oo/o,,, one has here 
*[ | * , ] | 
o*le, (eo — 1) +] — o%[e,(o- 1) -] = —-, 
. \ (50) 
| 
*/[ *f / = v | 
ot[p,(p — 1) +] — o%[p, (op — 1) —] = —__——. | 
(1 — v)p 


In contrast, the stress-discontinuities arising from the thermal shock in Danilovskaya’s 
problem [2] of the half-space are independent of the position of the plane wave-front 
and persist undiminished for all time. 

As the inertia parameter y approaches zero, while p and 7 are held fixed, the dynamic 
solution to the present problem should tend to the corresponding quasi-static solution. 
Carrying out this cumbersome limit process, which necessitates repeated use of the 
properties of y cited in (17), (18), (19), we obtain: 
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u* _ 


l 
a V(p, T), 
p \ 


‘ 1 
= 2 V(p, 7), o% = —3 V(p, r) — o(p, 7), 
p p ) 


where 


ap 


W(p, 7) = n(n, 7) dn 
1 


< 


(52) 
1/2 


(p> — 27 — 1)W(t) — (op + L)(r/m)'” exp (—&) + + + 2(7/z) 


Equations (51), (52) are identical with the quasi-static results deduced directly in [10]. 
As pointed out in [10], the displacement u* in (51) vanishes identically on the boundary 
p = 1, which remains fixed for all time. This curious conclusion fails to apply once inertia 
effects are taken into account. It should be mentioned that a quasi-static solution (in 
infinite series form) to the companion problem of a solid sphere whose surface temperature 
is suddenly altered, was given by Griinberg [15], as early as 1925; the same problem 
was reconsidered by Melan [16], [17]. Trostel [18] presented a formal quasi-static series 
solution to the general axisymmetric transient thermal-stress problem for a spherical 
shell of arbitrary thickness. 

Next, consider the limit of (15), (43), (45), (47), as r — ©, for fixed p and y. Here 
we find, after a tedious computation, 


‘ 
g=-, u* = (1 — p’), 


o*=p'—p', of = —-Hp + 2').J 


These values of displacement and stresses also characterize the limit of (51) as r — 0. 
Consequently, the dynamic and the quasi-static solution approach the same steady 
state. Formulas (53) are in agreement with the appropriate steady-state results for 
a spherical shell."’ 

Numerical results. Discussion. We turn, finally, to the discussion of certain numerical 
results based on the dynamic solution to the transient thermoelastic problem treated in 
this paper. Our main objective, in this connection, is to illustrate the character, and 
assess the magnitude, of the departures from the analogous quasi-static results. The 
space and time dependence of the temperature field (15), as well as of the quasi-static 
displacement and stresses (51), were discussed extensively in [10]; the corresponding 
diagrams need not be reproduced here. 

The dynamic solution, in contrast to (51), depends not only on p and 7, but involves 
also the inertia parameter y and Poisson’s ratio” v. An inspection of (43), (45), (47), 
together with the equations defining the various auxiliary functions, reveals that these 
formulas are expressible exclusively in terms of elementary functions and error functions 
of real argument, provided w = +r — y(p — 1) < 0. Forr > y(p — 1), however, that 
is, for points lying behind the wave-front, the displacement and stresses contain, in 


See, for example, [19], p. 420. See also [20]. 
12N ote that the dimensionless quantities u*, o* , and c¥ in (51) are independent of v. See (7). 
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addition, the complementary error function of complex argument ¥(k'*w'’’), in which 
k is the complex parameter (38). According to (16) and (38), 


oo 


ae 2 r - 9 ° 
V(kw)'? = = | exp (y' — 9 — 2iyn) dn, 


1/2 1/2 
~ — | & (91/2,1/72 =, | & (91/2,,1/2 
a E (2°"p »)| » Y [2 (2°) + »| , 


where p and gq are given by (28). 

While dense tabulations of ¥(z), for real z, are available in [21], the existing tables 
for complex z appear to be inadequate for our purposes.’* This fact, along with the 
general unwieldiness of (43), (45) and (47), complicates the numerical evaluation of 
the present results. The magnitude of y is bound to be exceedingly small compared to 
unity in applications which are of physical interest.’* It is natural, therefore, to seek 


oX 


(54) 











4 
-0.9- 
dynamic 
-1.0 1 l a — 
re) i 2 3 ' 
quasi-static steady - state 
-tif 
12+ 
“1.35 





Fic. 1. Time-dependence of o% at p = 1 for y = 1/5, vy = 1/4. 


13Approximate methods for the computation of the error function of complex argument are dis- 


cussed by Rosser [22]. 
M41Tf a = 1 in. and the material is steel, y is of the order of 10~’. 
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Fic. 2. Time-dependence of o% at p = 2 for y = 1/5, » 
asymptotic expansions for u*, o* , and o% , in the neighborhood of y = 0, with the aid of 
the semi-convergent development (19). Unfortunately, the asymptotic expansions thus 
reached break down at p = 1 and at r = 0; for this reason, they fail to supply useful 
approximations to the solution in the most significant range of position and time. 

Figures 1, 2, 3 show the time-dependence of o% at p = 1, 2, 3, respectively, for » = 1/4 
and y = 1/5. This unrealistically large value of the inertia parameter was chosen in 
order to bring out clearly the qualitative nature of the dynamic effects here involved. 
The underlying computations necessitated appropriate tabulations of the complex 
error function (54), which were carried out on a 650 I.B.M. electronic computer. The 
corresponding quasi-static curve (y = 0), has been included on each of the diagrams 
under discussion. 

In Fig. 1, the dynamic stress values are seen to undergo a pronounced oscillation, 
after which they approach the constant quasi-static value o$(1, 7) = — 1, which coincides 
with the steady-state boundary-value of o% . It is interesting to recall that the boundary- 
values of the transverse normal stress in Danilovskaya’s problem [2], [5], do not exhibit 
any dynamic effects. In Figs. 2, 3, the dynamic curves display the jump discontinuity 
(49) at the wave-front and thereafter rapidly approach the corresponding quasi-static 
curves which, in turn, tend to the common steady-state asymptote. As the shock-wave 
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Fic. 3. Time dependence of o% at p = 3 for y = 1/5, » = 1/4. 


progresses, the difference between the values of o% immediately behind and ahead of 
the wave-front decreases in accordance with (49); it should be kept in mind that this 
difference is independent of y. 

Figure 4 illustrates the quantitative influence of inertia upon the transient thermal 
stresses. Here, the dynamic peak-values of o% on either side of the wave-front, and the 
quasi-static values at the wave-front, are plotted as a function of r. The dynamic curves 
are based on vy = 1/4 and y = 7.3 X 10°°. This choice of the parameters is descriptive 
of steel, if a = 1 ft. The current value of y is obtained from (8), (9), by use of the approxi- 
mate values » = 5 X 10°° poundals/ft?, 8 = 484 lbs/ft®, x = 1.28 X 10™ ft?/sec. In 
this instance c = 1.75 X 10*ft/sec. Furthermore, for the preceding choice of a and x, 
according to (2), ¢ = 7.8 X 10°r. 

The curves appropriate to the dynamic o% behind the wave-front (w = 0+) and to 
the quasi-static o% at the wave-front (w = 0), are indistinguishable from the coordinate 
axes for the time-scale underlying the main diagram in Fig. 4. The inset diagram—drawn 
to a magnified time-scale—exhibits what might be called a boundary-layer effect with 
respect to time, which is characteristic of the behavior of o% during the initial stage 
of the motion. Thus, the dynamic peak of o% ahead of the wave-front is reduced to about 
one third of its initial limiting value of —3/2 during the first 4 X 10°” second of the 
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motion. In the same time-interval, the dynamic peak of c% behind the wave-front as 
well as the quasi-static o% at the front of the wave, decrease to less than five per cent 


of their joint initial value of —1. 

In conclusion we emphasize that the present results, which presuppose a step- 
function dependence upon time of the surface temperature, are bound to require severe 
modifications once the fiction of an instantaneous heating of the boundary is abandoned, 
even if the rate of heating, though finite, is extremely rapid by physically realistic 
standards.” 

Acknowledgment. The authors are indebted to Mr. J. J. Kalker who conducted that 
portion of the numerical work which was carried out on the I.B.M. 650 electronic 


computer. 
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toy. Astron. 


15Compare the conclusions reached in [5]. 
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(Continued from p. 184) 


naturally on the theory of the Hugoniot curve. Polachek and Seeger’s article is notable for the balance 
of its coverage and for the felicitous presentation of experimental records illustrating the theory. 

Stever’s discussion of condensation phenomena from the point of view of the wind-tunnel man in- 
cludes both a guide through the kinetic theory of condensation and a description of Oswatitsch and 
Wegener's theories of shocks and gradual condensation. It is rounded off by a discussion of the experi- 
mental evidence and its relation to the theory. 

After an illuminating introduction to one-parameter combustion by v. Kérmaén, Emmons discusses 
present knowledge on flames. No scope here for a definitive account, so Emmons concentrates instead 
on inviting the reader to join in the fray! The few topics singled out—flame structure, normal and 
oblique discontinuous fronts, vorticity production, steady flame geometry and beginnings of a sta- 
bility theory—are treated with elegance, and the gaps in present knowledge are vividly underlined. 

The chapter on detonations by Taylor and Tankin is distinguished from the others by the extensive 
use of numerical examples to support the physical argument. After a disucssion of the detailed compu- 
tation of a Hugoniot curve, the viewpoint foreshadowed in Hayes’ chapter is taken up that the Hugoniot 
relation for the detonation front needs to be complemented by a consideration of the unsteady motion 
of the burned gas behind the front, the possible types of stable front being determined in many cases by 
boundary conditions far behind the front. The same idea is applied to the case of a shock followed at 
some distance by a flame front, to explain one of the possible mechanisms of transition from deflagration 
to detonation, and the discussion of both cases is extended from the plane front to the spherical front 
moving at a constant speed. The chapter closes with a brief description of the theory of spinning de- 
tonation and its experimental support, and of the tentative extension of the gas-dynamical theory to 
explosions in liquids and solids. 

The chapter on rarefied gas dynamics describes the situation at the time of writing and some of the 
reasons why the situation will be very different very soon. A brief survey of the method and main 
results of free molecule flow calculations is followed by one of the derivation of the Thirteen Moment 
and Burnett equations, but Schaaf and Chambré make no secret of their feeling—so well confirmed 
since—that their study should be postponed in favor of that of the good old Navier-Stokes equations. 
To give this chapter at least maximum short-term usefulness, they add an extensive survey of present 
experimental results, and a bibliography which surpasses even those given by the other authors. 

With two exceptions, this collection of ten articles meets the specification—that it should not only 
be, but also remain for many years, the foremost reference work in its field. 


R. E. MEYER 


Jets, wakes, and cavities. By G. Birkhoff and E. H. Zarantonello. Academic Press, Inc., 
New York, 1957. xii + 353 pp. $10.00. 


This book is a noteworthy essay in bringing together the known theories and facts concerning the 
flows described in the title, all of which are characterized by the presence of free streamlines or streamlines 
of discontinuity. For this, if for.no other reason, the authors are to be congratulated in filling a real 
gap in the existing literature of fluid motion. 

To quote from the preface “our book draws on the resources of pure and applied mathematics and 
on experimental physics, and it sheds light on numerous problems of hydraulics and aeronautics. There- 
fore, it will perhaps have the greatest interest for readers whose scientific curiosity spans all the fields 
just mentioned. However, we hope that others will also find it a useful and stimulating reference in 
connection with many special questions”’. 

Owing to the development of high-speed computing methods which enable many calculations to be 
performed which hitherto would have taken a prohibitive amount of time, approximations are now 
the order of the day. It may therefore be permissible to say that the present work falls, using the authors’ 
delightful phrase (page 220), “into two roughly equal halves’; the first concerned with exact mathe- 
matical methods applied to the two-dimensional motion of an inviscid liquid, the second with the sig- 











220 BOOK REVIEWS [Vol. XVII, No. 2 


nificant theory of axisymmetric jets and cavities which has developed in the last two decades, and the 
reinterpretation of the fundamental facts about vortex trails, and turbulent jets and wakes. 

The exact two-dimensional theory of Chapters II-IX proceeds naturally by use of the complex 
variable and is well-classified to permit the application of uniform methods, but one would like to have 
seen at least one problem involving gravity flow pursued to a conclusion instead of stopping (page 204) 
at the point where it becomes interesting. The most valuable contributions of these chapters are contained 
in Chapter IV where there are brought together accounts of recent work on comparison of flows originated 
by Lavrentieff and simplified in a series of papers by Gilbarg and Serrin with applications to uniqueness 
and minimum cavity drag, and in Chapter VII on existence and uniqueness. 

The remaining chapters concern (X) axially symmetric flows, (XI) unsteady potential flows, (XII) 
steady viscous wakes and jets, (XIII) periodic wakes, (XIV) turbulent wakes and jets, (XV) miscel- 
laneous experimental facts. These chapters treat generally increasingly complicated situations. At each 
stage of complication the authors try to give a picture of what is known, organized as rationally as 
possible. Here rigorous mathematical theory is not always forthcoming and hypothesis and empirical 
methods play a large part. Nevertheless, these chapters constitute an admirable and absorbing part 


of the book. 
L. M. Mitne-TuHomson 


Elements of gas dynamics. By H. W. Liepmann and A. Roshko. John Wiley & Sons, Inc., 
New York, and Chapman & Hall, Ltd., London, 1957. xv + 439 pp. $11.00. 


Many of us engaged in teaching high speed aerodynamics have eagerly awaited the publication of 
this book for some time now. It is the opinion of this reviewer that a book as excellent as this one was well 
worth waiting for. The author’s aim was to provide a comprehensive text for senior-first-year graduate 
study in modern gas dynamics, and in this respect they could not have achieved greater success. The 
presentation is almost always novel; the material is discussed in the clearest physical and mathematical 
terms; there is never any loss in rigor; and yet, there is a clarity which makes it ideal for student use. 
If anything, the student may be led to believe that many of the topics covered are too elementary, 
however this feeling should be dispelled by the fine set of problems which require much more than just 
“putting in numbers” to obtain a solution. 

Perhaps one of the best features of the book is that it is as modern as the latest research in the 
field of gas dynamics. Thus, in Chapter 1 the concepts of thermodynamics are introduced including the 
thermodynamics of reacting and dissociating gas mixtures, a knowledge of which is so essential for the 
study of hypersonic flow. Chapter 2 considers steady one-dimensional gas flow, and Chapter 3 unsteady 
one-dimensional wave motion along with a discussion of shock tubes. As throughout, experimental results 
and photographs are included to give the reader a vivid illustration of the “comparison of theory with 
experiment.” Chapter 4 is primarily concerned with wave and wave interaction phenomena in steady 
two-dimensional supersonic flow, although more general problems such as the detached shock wave 
are also considered. Flows in ducts and wind tunnels are considered in Chapter 5, while Chapter 6 out- 
lines methods of measuring supersonic flow phenomena. In addition to discussing the now conventional 
techniques such as schlieren photography, pitot measurements, etc., the more unconventional and 
modern X-ray absorption method, direct skin friction technique, and hot-wire methods, are also pre- 
sented. In Chapter 7 the equations of frictionless flow are derived, along with general flow theorems, 
and a discussion of the character of the equations of motion. The small perturbation equations and 
boundary conditions are derived in Chapter 8 for subsonic, transonic, supersonic, and hypersonic flow, 
and Ackeret’s two-dimensional thin airfoil theory is presented. Chapter 9 considers supersonic linear- 
ized potential flow over bodies of revolution and a clear distinction, so usually lacking, is made between 
first-order and slender body theory. All of the similarity laws for the high speed flow regimes are derived 
in Chapter 10. Chapter 11 is devoted to the elements of transonic flow and, as in the rest of the book, 
a very clear physical picture is given of the problems involved in this difficult field. In Chapter 12 the 
method of characteristics is developed by using intrinsic or natural coordinates, so that the discussion 
is unencumbered by the usual mass of algebra associated with a presentation of the method in cartesian 
coordinates. An extremely concise but nonetheless clear chapter is 13 which considers the effects of 
viscosity and conductivity both for laminar and turbulent flow. Here the authors introduce the problems 
associated with boundary layers, including the effects of dissociation, by the use of the simple and 
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instructive Couette flow. The most recent work in boundary layer-inviscid flow interaction phenomena 
are also considered. In the closing Chapter 14 the concepts of gas kinetics of importance to aerodynamics 
are considered, starting with the fundamental idea of probability concepts, and including such topics 
as the flow of highly rarefied gases, relaxation phenomenon, and continuum theory limits. 

As is always true with a book as comprehensive as this one, there are some minor criticisms. For 
example, the book would probably have been enhanced by more adequate referencing. In addition, 
this reviewer would take exception with some of the views expressed in the topics which are now actively 
being studied, such as boundary layer-induced interactions and continuum theory limits. However, 
such comments can indeed be considered trivial in the light of the excellence of the book, including, 
by the way, the format and clarity of printing. Apart from a “must” text in any theoretical aerodynamics 
course, this reviewer would also recommend the book for anyone involved in high speed aerodynamics 
research, as well as for those persons with a suitable background who wish to obtain an overall picture 
of modern gas dynamics. The authors say that they hope this book will be followed by one of a more 
advanced and specialized nature—we can only look forward to its publication with the greatest interest. 


RONALD F. PROBSTEIN 


The calculation of atomic structures. By Douglas R. Hartree. John Wiley & | ons, Inc., 
New York, and Chapman & Hall, Ltd., London, 1957. xiii + 181 pp. $5.00. 


This book is based on a series of lectures given by the author at Haverford College in 1955. The 
practical methods of calculating atomic structures are presented in a way to be of immediate help to 
anyone concerned with such calculations. 

Unfortunately, the results of the calculations on atomic structures are not given but a list of refer- 
ences to such results is given in an appendix. Recent work in calculations of this type has been appreciable 
with thirty or more papers appearing during the past seven or eight years. The fact that no really ade- 
quate discussion of these calculation methods has previously appeared in one place should make this 


book most welcome. 
Roun TRUELL 


Rheology, theory and applications. Edited by Frederick R. Eirich. Academic Press, Inc., 
New York, 1956. xiii + 761 pp. $20.00. 


This volume, the first of three intended to provide “integrated surveys’’ of ‘‘well-demarcated areas 
of rheology”, described as ‘‘the science of deformation and flow,” contains seventeen articles on an 
impressive range of topics. Most of the work is of an experimental, physical, or chemical nature and need 
not be described here. Furthermore, most if not all of the articles concerning mechanics or applied math- 
ematics consist in material easily available in other surveys, often by the same authors, and it is possible 
that the main virtue of the volume is to put between two covers a number of different conceptions of 
this fluid science. 

“Phenomenological macrorheology” by M. Reiner and ‘Dynamics of viscoelastic behavior’ by 
Turner Alfrey, Jr., and E. F. Gurnee deal with what seems to be the same topic from the mechanical 
point of view and differ in that the former emphasizes terminology, the latter, the consequences of 
linearity and the associated operator techniques. While both articles contain passages filled with tensor 
equations, the methods and problems are essentially one-dimensional, and prominence is given to the 
devising of complicated materials by imagined couplings of simpler materials. 

Only three articles approach the mathematical theories of genuine large deformation and flow from 
a standpoint that seems tenable in view of recent researches in general mechanics. J. G. Oldroyd’s 
“Non-newtonian flow of liquids and solids’ deals mainly with rectilinear flows. As is now well known, 
specialization of correct three-dimensional equations to such flows does not yield the same results as 
were assumed in the purely one-dimensional treatments of the older studies on rheology. The article of 
Oldroyd, however, concentrates upon its author’s researches of some years ago and cannot be regarded 
as giving a just view of the field today. While the general equations for an incompressible fluid are the 


author’s equations (46), viz. 


Pik = Dix + sy z. D;;D jx - Pox ’ 
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virtually all of the article is restricted to the case y = 0; that is, to a fluid which might be described as 
devoid of cross-viscosity though possessed of variable shear viscosity. Rivlin’s general solutions for 
y =~ Oare alluded to but not described. Moreover, more recent researches of Ericksen! and Stone? indicate 
that rectilinear flow generally fails to exist if ¥/n * const.; in a special case, Green and Rivlin® have 
obtained an approximate solution in which a circulating cross flow is superimposed upon the flow down 
the tube. 

William Prager’s “Finite plastic deformation’’ is an interesting exception to the usual literature 
of plasticity in that it attempts to select those theoretical problems which may be interpreted as in- 
volving large deformation. Except for a few pages at the end, only the case of a body rigid up to yield 
and perfectly plastic thereafter is considered; even here, ‘‘finite’’ must be qualified, since the inertial 
terms in the equations of motion are neglected. This article is very clearly written and is illuminating 
for its well chosen discussion, examples, and figures. Most of it concerns plane flows, and the complicated 
calculations are omitted in favor of explanation of the results in a form comprehensible to those who are 
not specialists in plasticity. 

“Large elastic deformations” by R. S. Rivlin attempts a task less ambitious than those set them- 
selves by the other authors, since it is devoted to a definite and well poised problem of mechanics. Its 
value lies principally in explaining to those who might wish to apply it the general theory of elasticity 
in the form and status a considerable body of specialists now see for it. The article presents only the parts 
of the theory where the author has made notable contributions; this results in a certain lack of apprecia- 
tion of other aspects, and in particular it seems to me that the work of Signorini does not deserve to be 
dismissed along with that of Seth on p. 353, and that the general approximation process initiated on 
pp. 382-383 was given, substantially, twenty years ago by Signorini. The reader of Rivlin’s article will 
come away with an accurate summary of certain definite predictions of finite elasticity, but if he has no 
more mathematical apparatus than the author requires of him, it is unlikely that he can read further 
in the literature of the subject. 

Another exceptional article, of a very different kind, is ““The statistical mechanical theory of irrevers- 
ible processes in solutions of macromolecules,” by J. Riseman and J. G. Kirkwood. This is a condensed 
extract from one aspect of the group of statistical theories proposed by Kirkwood and his students. A 
much higher level of mathematical and physical preparation is needed for the reader of this article, 
which concerns statistical mechanics of chains and employs a Riemannian geometry for the subspace 
obtained by imposing constraints on the motions in the molecular configuration space. Approximation 
procedures are introduced to get definite results for special models. 

Looking back at the volume as a whole, while we find much overlap of qualitative discussion, es- 
pecially regarding names and symbols, at the same time this heavy tome seems to have found no space 
to present the recent exact studies of visco-elastic and time rate processes. The fundamental memoir 
of Oldroyd‘ is cited thrice (pp. 51, 657, 676), but in such a way that one not already familiar with it 
could scarcely divine its subject; the memoirs of Noll5 and of Rivlin and Ericksen*® perhaps appeared 
too late for mention. However, it seems unlikely that fundamental problems of large deformation can 
ever be described successfully, even a posteriori, at the mathematical level most of the surveys in this 


volume strive to maintain. 


1Q. Appl. Math. 14, 318-321 (1956). 

7Q. Appl. Math. 15, 257-262 (1957). 

*Q. Appl. Math. 14, 299-308 (1956). 

‘Proc. Roy. Soc. London A 200, 523-547 (1950). 
5J. Rational Mech. Anal. 4, 3-81 (1955). 


6J. Rational Mech. Anal. 4, 323-425 (1955). 
C. TRUESDELL 


Experimental designs in industry. Edited by Victor Chew. John Wiley & Sons, Inc., 
New York, and Chapman «& Hall, Ltd., London, 1958. xi + 268 pp. $6.00. 


During recent years experimental designs have been employed at an increasing rate in industrial 
research and development work. The various classical designs as Latin squares, balanced incomplete 
block designs, etc., although they have their place in industrial research are quite inadequate for many 
problems that arise in industry where often the number of factors and the number of levels of each 
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factor are far too large to work with complete replications. The fractional replications introduced by 
Finney in 1945 proved particularly useful in this situation. The methods of fitting response surfaces 
first. introduced by Box and Wilson in 1951 seem to have proven their value in situations where the 
factors may be applied in continuously varying levels. 

It was thought that these designs had been used for a long enough time to exchange information 
on experiences with their use. To this purpose a symposium was held at North Carolina State College 
on November 5-9, 1956. The present volume contains a selection of papers presented at this symposium. 
The first paper by V. Chew gives a general survey of the basic ideas of the analysis of variance. The 
next paper by R. L. Anderson deals with complete and fractional factorials and the technique of con- 
founding. This is followed by a paper on multiple regression analysis by R. J. Hader and A. H. E. Gran- 
dage. An exposition by J. E. P. Box and J. 8. Hunter of the methods of fitting response surfaces intro- 
duced by Box concludes the theoretical part of the volume. The second part reports on experiences 
in the application of designs in industry in papers by W. S. Connor, W. H. Horton, M. B. Carroll and 
O. Dykstra Jr., De Baun and A. M. Schneider and concludes with a report by C. A. Bicking on experi- 
ences with designs and needs for designs in ordnance experimentation. 

Although all papers presented will be of interest to anyone in the field of mathematical statistics, 
the exposition of the methods of fitting response surfaces seems of particular significance. It seems 
to the reviewer that these methods introduce ideas that are not only basically new but have also arisen 
from a very urgent practical need. 

It is of course not the purpose of the papers presented in this volume to give a complete account 
of the mathematical foundations of the methods presented. The authors therefore restrict themselves 
to a presentation of the calculation procedures and a description of the situations in which the methods 
apply. Numerous applications to problems that have actually arisen in industrial research are presented. 
For those who are interested in a more thorough study a large index of references is given at the end 
of each paper. 


H. B. Mann 


Principles of noise.By J. J. Freeman. John Wiley & Sons, Inc., New York, and Chapman 
& Hall, Ltd., London, 1958. x + 299 pp. $9.25. 


This book is designed to introduce graduate students, principally in electrical engineering, to the 
fundamentals of noise analysis. Thus, the mathematical tools—Fourier series and integrals, probability 
theory, stochastic processes—are developed in sufficient detail for the purpose and the physical sources 
of noise are described before proceeding to their analysis. The work may also serve as an elementary 
introduction to several recent more advanced treatises on random signals and noise. There are many 


well chosen examples and problems for the reader. 
W. F. FrerperGeR 


The preparation of programs for an electronic digital computer. By Maurice V. Wilkes, 
David J. Wheeler, and Stanley Gill. Addison-Wesley Publishing Co., Reading, Mass., 
1957. xiv + 238 pp. $7.50. 

When the first edition of this work was published in 1951 it represented the first effort to make the 
art of programming accessible in bookform and formed an invaluable aid in spreading knowledge of 
a relatively new field. Since then, developments have been rapid and books on computers are rolling 
off the presses. Thus, this second edition of a pioneering treatise inevitably seems somewhat outdated. 
The language in which the programs are presented is still the order code of the EDSAC and although 
material on other machines is included, it is this code the reader will have to master to understand the 
book. One of the most striking developments of recent years, at least in America, has been the increasing 
use of automatic coding systems, and there is no discussion of this interesting and important subject 
in the book. 

This is, thus, strictly a book for the specialist who will make the effort of mastering it for the sake 
of the innumerable hints on methods and techniques strewn throughout its pages; the authors’ wide 


experience vouchsafes their effectiveness. 
W. F. FRerBeRGER 





TENTH INTERNATIONAL CONGRESS OF APPLIED MECHANICS 


The Tenth International Congress of Applied Mechanics will be held in the Congress 
Building at Stresa (Italy) from Wednesday, 31 August through Wednesday, 7 September 
1960. 

Apart from a number of invited general lectures the technical sessions of the Con- 
gress will be held in two sections, viz.: 


Section 1: Fluid dynamics (hydrodynamics and aerodynamics). 
Section 2: Mechanics of solids (rigid body dynamics, vibrations, elasticity, plas- 
ticity and theory of structures). 


It should be noted that thermodynamics and computational methods as such are not 
included, although specific applications of computational methods to pertinent problems 
of one of the two sections mentioned above are acceptable subjects for papers to be read 
at the Tenth Congress. 

Previous Congresses have demonstrated the desirability of an adequate period of 
time for the presentation and discussion of individual papers. In order to allow a period 
of 45 minutes for each paper (30 mins. for presentation and 15 mins. for discussion) a 
Program Committee will make a selection from papers submitted for presentation. 
Abstracts of papers should be submitted in four copies to the Secretary of the Interna- 
tional Committee (Prof. Mekelweg 2, Delft, Netherlands) before 1 January 1960. Prefer- 
ably they should not exceed two typewritten pages (double-spaced) and in no case 
should they exceed four pages. In order to facilitate the work of the Program Committee 
it is recommended that abstracts be in two of the official Congress languages (English, 
French, German and Italian). Authors are urged to make their abstracts as clear as 
possible, since selection of papers has to be based upon them. Decisions of the Program 
Committee are final, and it will be understood that it is impossible to enter into correspond- 
ence about them with authors of papers. They will be informed promptly of the deci- 
sion on each paper. 

Day-to-day organization of the Congress is effected by the Italian Organizing Com- 
mittee (President: Prof. G. Colonnetti; Secretary: Dr. F. Rolla, Consiglio Nazionale 
delle Ricerche, Ufficio relazioni internazionali, Piazza delle Scienza 7, Roma). All corre- 
spondence (apart from submission of papers) should be addressed to the Italian Organ- 
izing Committee. Information on accommodation, also registration forms, will be obtain- 
able from Dr. Rolla on and after 1 September 1959. 


The Executive Committee of the 
International Committee for the 
Congresses of Applied Mechanics 
C. B. Brezeno, President; 
RICHARD V. SOUTHWELL; 
W. T. Korrer, Secretary; 
(Prof. Mekelweg 2, Delft). 




















